def main(instance, model_file):
print("Server: " + str(platform.node()) + "\n")
# Obtem os valores do arquivo e cria os retângulos
filename = instance
boxes, Lu, Wu, Hu = get_data(filename)
filename = filename.replace('.data', '', 1)
filename = filename.replace('instances/', '', 1)
print(filename)
m = read(model_file)
print('Problema: ' + m.getAttr('modelName'))
m.setParam('TimeLimit', 1 * 3600.0)
vars = m.getVars()
z_obj, L, W, H = vars[0], vars[1], vars[2], vars[3]
X, Y, Z = [], [], []
for i, box in enumerate(boxes, start=0):
li = box.get('dim')[0]
wi = box.get('dim')[1]
hi = box.get('dim')[2]
X_i = [p for p in range(0, Lu - li)]
Y_i = [p for p in range(0, Wu - wi)]
Z_i = [p for p in range(0, Hu - hi)]
X.append(X_i)
Y.append(Y_i)
Z.append(Z_i)
m.optimize()
# Recuperar as variáveis binárias
x_bin = []
for i, box in enumerate(boxes):
rotations = list(permutations(box.get('dim')))
k_list = []
for k, b in enumerate(rotations):
p_list = []
for p in X[i]:
q_list = []
for q in Y[i]:
r_list = []
for r in Z[i]:
index = 'X_' + '{}_{}_{}_{}_{}'.format(i, k, p, q, r)
x = m.getVarByName(index)
r_list.append(x)
q_list.append(r_list)
p_list.append(q_list)
k_list.append(p_list)
x_bin.append(k_list)
print(' = Variables Added')
if True:
m.write('results/resultado_' + filename + '.sol')
print('Container: ', Lu, Wu, Hu, sep='\t')
container = [Lu, Wu, z_obj.X]
print('\nFunção objetivo: ',
container[0],
container[1],
container[2],
sep='\t')
# Recupera as posições das caixas
solutions = []
for i, box in enumerate(boxes, start=0):
sol_box = []
for k in range(6):
for p in X[i]:
for q in Y[i]:
for r in Z[i]:
if x_bin[i][k][p][q][r].X == 1:
sol_box.append([i, (p, q, r), k])
solutions.append(sol_box)
ret = []
for i, box in enumerate(boxes):
color_b = np.random.rand(1, 3)
rotated_box = list(permutations(box.get('dim')))
for j in range(box.get('qtd')):
sol = solutions[i][j]
r = rotated_box[sol[2]]
r = Retangulo(vertex=r, color=color_b)
r.change_vertex(list(sol[1]))
ret.append(r)
ret.sort(key=lambda x: max(x.vertex[2]), reverse=True)
file_write_list = open('results/retangulos_' + filename + '.bin', 'wb')
pickle.dump(ret, file_write_list)
visualiser = PolyVisualiser(
array_polygons=ret,
x_lim=container[0],
y_lim=container[1],
z_lim=container[2],
obj=container[2],
alpha=.8,
title=filename,
)
visualiser.animate(no_animation=False)
visualiser.scrol()
visualiser.show()
def main(instance):
print("Server: " + str(platform.node()) + "\n")
# Obtem os valores do arquivo e cria os retângulos
# filename = 'problema1.data'
filename = instance
boxes, Lu, Wu, Hu = get_data(filename)
filename = filename.replace('.data', '', 1)
filename = filename.replace('instances/', '', 1)
print(filename)
# Cria os conjuntos X,Y,Z
BigX, BigY, BigZ = [], [], []
for i in boxes: # Percorrer cada tipo de caixa
p, q, r = 0, 0, 0
for e in range(1, i.get('qtd') + 1):
l = i.get('dim')[0]
p += e * l
for e in range(1, i.get('qtd') + 1):
l = i.get('dim')[1]
q += e * l
for e in range(1, i.get('qtd') + 1):
l = i.get('dim')[2]
r += e * l
BigX.append(p)
BigY.append(q)
BigZ.append(r)
container = [Lu, Wu, Hu]
BigX = [i for i in range(sum(BigX))]
BigY = [i for i in range(sum(BigY))]
BigZ = [i for i in range(sum(BigZ))]
# print(BigX, BigY, BigZ, sep='\n')
m = Model("3D-ODRPP")
m.setParam('TimeLimit', 1 * 3600.0)
L = m.addVar(vtype=GRB.CONTINUOUS, name="L", lb=0, ub=GRB.INFINITY)
W = m.addVar(vtype=GRB.CONTINUOUS, name="W", lb=0, ub=GRB.INFINITY)
H = m.addVar(vtype=GRB.CONTINUOUS, name="H", lb=0, ub=GRB.INFINITY)
z_obj = m.addVar(vtype=GRB.CONTINUOUS, name="Z", lb=0, ub=GRB.INFINITY)
m.addConstr(z_obj == L + W + H, "z_objective")
m.setObjective(z_obj, GRB.MINIMIZE)
X, Y, Z = [], [], []
for i, box in enumerate(boxes, start=0):
li = box.get('dim')[0]
wi = box.get('dim')[1]
hi = box.get('dim')[2]
X_i = [p for p in range(0, Lu - li)]
Y_i = [p for p in range(0, Wu - wi)]
Z_i = [p for p in range(0, Hu - hi)]
X.append(X_i)
Y.append(Y_i)
Z.append(Z_i)
# print(X, Y, Z, sep='\n')
# Criação das variáveis binárias
x_bin = []
for i, box in enumerate(boxes):
rotations = list(permutations(box.get('dim')))
# print(box.get('dim'), rotations, sep='\t')
k_list = []
for k, b in enumerate(rotations):
p_list = []
for p in X[i]:
q_list = []
for q in Y[i]:
r_list = []
for r in Z[i]:
index = '[' + str(i) + '][' + str(k) + '][' + str(
p) + '][' + str(q) + '][' + str(r) + ']'
x = m.addVar(vtype=GRB.BINARY, name='X_' + index)
r_list.append(x)
q_list.append(r_list)
p_list.append(q_list)
k_list.append(p_list)
x_bin.append(k_list)
print(' = Variables Added')
# Constraint 2
for s in BigX:
for t in BigY:
for u in BigZ:
soma_r2 = LinExpr()
for i, box in enumerate(boxes, start=0):
rotations = list(permutations(box.get('dim')))
for k, b in enumerate(rotations):
li, wi, hi = b[0], b[1], b[2]
for p in X[i]:
for q in Y[i]:
for r in Z[i]:
if s - li + 1 <= p <= s:
if t - wi + 1 <= q <= t:
if u - hi + 1 <= r <= u:
soma_r2 += x_bin[i][k][p][q][r]
# index = str(i)
index = '[' + str(s) + '][' + str(t) + '][' + str(u) + ']'
m.addConstr(soma_r2 <= 1, name='Restricao2_' + index)
print(' = Constraint 2 adicionada')
# Constraint 3
for j, box in enumerate(boxes, start=0):
bj = box.get('qtd')
soma_r3 = LinExpr()
rotations = list(permutations(box.get('dim')))
for k, b in enumerate(rotations):
for p in X[j]:
for q in Y[j]:
for r in Z[j]:
soma_r3 += x_bin[j][k][p][q][r]
m.addConstr(soma_r3 == bj, name='Restricao3_' + str(j))
print(' = Constraint 3 adicionada')
# Restrições 4, 5 e 6
for i, box in enumerate(boxes, start=0):
rotations = list(permutations(box.get('dim')))
for k, b in enumerate(rotations):
for p in X[i]:
for q in Y[i]:
for r in Z[i]:
index = '[' + str(i) + '][' + str(k) + '][' + str(
p) + '][' + str(q) + '][' + str(r) + ']'
li, wi, hi = b[0], b[1], b[2]
algo_p = (p + li) * x_bin[i][k][p][q][r]
algo_q = (q + wi) * x_bin[i][k][p][q][r]
algo_r = (r + hi) * x_bin[i][k][p][q][r]
m.addConstr(algo_p <= L, name='Restricao4_' + index)
m.addConstr(algo_q <= W, name='Restricao5_' + index)
m.addConstr(algo_r <= H, name='Restricao6_' + index)
print(' = Restrições 4,5,6 adicionadas')
m.write('models/modelo_' + filename + '.lp')
print(' = Model written!')
m.optimize()
if m.status == GRB.Status.OPTIMAL:
m.write('results/resultado_' + filename + '.sol')
print('Container: ', Lu, Wu, Hu, sep='\t')
print('\nFunção objetivo: ',
str(round(L.X)),
str(round(W.X)),
str(round(W.X)),
sep='\t')
# Recupera as posições das caixas
solutions = []
for i, box in enumerate(boxes, start=0):
sol_box = []
for k in range(6):
for p in X[i]:
for q in Y[i]:
for r in Z[i]:
if x_bin[i][k][p][q][r].X == 1:
sol_box.append([i, (p, q, r), k])
solutions.append(sol_box)
# solutions = [[[0, (3, 0, 0), 0], [0, (4, 0, 0), 0]], [[1, (0, 0, 0), 1], [1, (0, 3, 0), 1]],
# [[2, (1, 0, 4), 2], [2, (2, 0, 4), 4]], [[3, (1, 1, 5), 1], [3, (2, 2, 5), 4], [3, (1, 3, 5), 5]]]
ret = []
for i, box in enumerate(boxes):
color_b = np.random.rand(1, 3)
rotated_box = list(permutations(box.get('dim')))
for j in range(box.get('qtd')):
sol = solutions[i][j]
r = rotated_box[sol[2]]
r = Retangulo(vertex=r, color=color_b)
r.change_vertex(list(sol[1]))
ret.append(r)
ret.sort(key=lambda x: max(x.vertex[2]), reverse=True)
file_write_list = open('results/retangulos_' + filename + '.bin', 'wb')
pickle.dump(ret, file_write_list)
