Claudemir Woche titowoche30

🎯

Focusing

Always learning 🧐

titowoche30 / TSP.py

Created November 19, 2019 00:39

titowoche30 / TSP.py

Created November 19, 2019 00:31

	#Solves
	solver = pyEnv.SolverFactory('cplex')
	result = solver.solve(model,tee = False)

	#Prints the results
	print(result)

titowoche30 / TSP.py

Created November 19, 2019 00:27

	#Prints the entire model
	model.pprint()

titowoche30 / TSP.py

Created November 19, 2019 00:15

	def rule_const3(model,i,j):
	if i!=j:
	return model.u[i] - model.u[j] + model.x[i,j] * n <= n-1
	else:
	#Yeah, this else doesn't say anything
	return model.u[i] - model.u[i] == 0

	model.rest3 = pyEnv.Constraint(model.U,model.N,rule=rule_const3)

titowoche30 / TSP.py

Created November 19, 2019 00:14

	def rule_const2(model,N):
	return sum(model.x[N,j] for j in model.M if j!=N) == 1

	model.rest2 = pyEnv.Constraint(model.N,rule=rule_const2)

titowoche30 / TSP.py

Created November 19, 2019 00:07

	def rule_const1(model,M):
	return sum(model.x[i,M] for i in model.N if i!=M ) == 1

	model.const1 = pyEnv.Constraint(model.M,rule=rule_const1)

titowoche30 / TSP.py

Created November 19, 2019 00:02

	def obj_func(model):
	return sum(model.x[i,j] * model.c[i,j] for i in model.N for j in model.M)

	model.objective = pyEnv.Objective(rule=obj_func,sense=pyEnv.minimize)

titowoche30 / TSP.py

Created November 19, 2019 00:01

	#Cost Matrix cij
	model.c = pyEnv.Param(model.N, model.M,initialize=lambda model, i, j: cost_matrix[i-1][j-1])

titowoche30 / TSP.py

Created November 19, 2019 00:00

	#Decision variables xij
	model.x=pyEnv.Var(model.N,model.M, within=pyEnv.Binary)

	#Dummy variable ui
	model.u=pyEnv.Var(model.N, within=pyEnv.NonNegativeIntegers,bounds=(0,n-1))

titowoche30 / TSP.py

Created November 18, 2019 23:59

	#Model
	model = pyEnv.ConcreteModel()

	#Indexes for the cities
	model.M = pyEnv.RangeSet(n)
	model.N = pyEnv.RangeSet(n)

	#Index for the dummy variable u
	model.U = pyEnv.RangeSet(2,n)