def RunTest(network_graph, model_original, model_decomposer, solvers,
solve_original = False, solve_decomposed = False, max_iter = 100,
validate_feasability = False, recover_feasible = False, recover_feasible_original = False,
draw_progress = False, draw_solution = False):
#overall results connection
results = {}
#solve decomposed problem
if solve_decomposed:
#solve decomposed problem
solver_master = solvers['Master']
solvers_local = solvers['Decomposed']
start_time_decomposed = t.perf_counter()
solution = model_decomposer.Solve(solver_master, solvers_local, max_iter)
elapsed_time_decomposed = t.perf_counter() - start_time_decomposed
#show decomposition iterating data
if draw_progress:
data = model_decomposer.RecordedData
dr.PlotIterationData('Objective', data['MasterObj'])
for name, data in data['Multipliers'].items():
dr.PlotIterationData(name, data)
if solution is None:
results['Decomposed'] = None
print('###################!RESULTS!#############################')
print('Decomposed:\nSolution was not found!')
print('__________________________________________________________')
else:
objective, objective_dual, strains, routes, time = ( solution['Objective'], solution['ObjectiveDual'],
solution['Strain'], solution['Route'],
solution['Time'] )
results['Decomposed'] = { 'Objective': objective, 'ObjectiveDual': objective_dual, 'Time': time,
'Total time': elapsed_time_decomposed}
print('###################!RESULTS!#############################')
print(f'Decomposed:\nObjective: {objective}, Objective Dual: {objective_dual}, Time: {time}, Total time: {elapsed_time_decomposed}')
#validate constraints violations
if validate_feasability:
violations = mp.FindConstraintsViolation(model_original.cmodel, strains, routes)
cc_vn = violations['capacity_constraints'][1]
rc_vn = violations['route_constraints'][1]
if cc_vn == 0 and rc_vn == 0:
print('Feasible')
else:
print(f'Capacity constraint violations number: {cc_vn}')
print(f'Route constraint violations number: {rc_vn}')
print('__________________________________________________________')
#draw if needed
if draw_solution:
path_list_sol = [ [edge for edge, value in route.items() if abs(value - 1) <= 1e-6] for route in routes ]
flow_list_sol = strains
network_graph.SetPath(*path_list_sol)
network_graph.SetFlows(*flow_list_sol)
dr.PlotNetwork(network_graph, 'DECOMPOSED')
#recover feasible for given routes
if recover_feasible:
solver = solvers['Recovered']
start_time_recovered = t.perf_counter()
solution = mp.RecoverFeasibleStrain(model_original, routes, solver)
elapsed_time_recovered = t.perf_counter() - start_time_recovered
if solution is None:
results['Recovered'] = None
print('###################!RESULTS!#############################')
print('Recovered:\nSolution was not found!')
print('__________________________________________________________')
else:
objective, strains, routes, time = ( solution['Objective'], solution['Strain'],
solution['Route'], solution['Time'] )
results['Recovered'] = { 'Objective': objective, 'Time': time, 'Total time': elapsed_time_recovered }
print('###################!RESULTS!#############################')
print(f'Recovered:\nObjective: {objective}, Time: {time}, Total time: {elapsed_time_recovered}')
#validate constraints violations
if validate_feasability:
violations = mp.FindConstraintsViolation(model_original.cmodel, strains, routes)
cc_vn = violations['capacity_constraints'][1]
rc_vn = violations['route_constraints'][1]
if cc_vn == 0 and rc_vn == 0:
print('Feasible')
else:
print(f'Capacity constraint violations number: {cc_vn}')
print(f'Route constraint violations number: {rc_vn}')
print('__________________________________________________________')
#draw if needed
if draw_solution:
path_list_sol = [ [edge for edge, value in route.items() if abs(value - 1) <= 1e-6] for route in routes ]
flow_list_sol = strains
network_graph.SetPath(*path_list_sol)
network_graph.SetFlows(*flow_list_sol)
dr.PlotNetwork(network_graph, 'RECOVERED')
#solve original problem
if solve_original:
network_graph = cp.deepcopy(network_graph)
solver = solvers['Original']
start_time_original = t.perf_counter()
solution = solver.Solve(model_original.cmodel)
elapsed_time_original = t.perf_counter() - start_time_original
if solution is None:
results['Original'] = None
print('###################!RESULTS!#############################')
print(f'Original:\nSolution was not found!, Total time: {elapsed_time_original}')
print('__________________________________________________________')
else:
objective, strains, routes, time = ( solution['Objective'], solution['Strain'],
solution['Route'], solution['Time'] )
results['Original'] = { 'Objective': objective, 'Time': time, 'Total time': elapsed_time_original }
print('###################!RESULTS!#############################')
print(f'Original:\nObjective: {objective}, Time: {time}, Total time: {elapsed_time_original}')
#validate constraints violations
if validate_feasability:
violations = mp.FindConstraintsViolation(model_original.cmodel, strains, routes)
cc_vn = violations['capacity_constraints'][1]
rc_vn = violations['route_constraints'][1]
if cc_vn == 0 and rc_vn == 0:
print('Feasible')
else:
print(f'Capacity constraint violations number: {cc_vn}')
print(f'Route constraint violations number: {rc_vn}')
print('__________________________________________________________')
#recover feasible for given primal routes
if recover_feasible_original:
solver = solvers['Recovered']
start_time_recovered = t.perf_counter()
solution = mp.RecoverFeasibleStrain(model_original, routes, solver)
elapsed_time_recovered = t.perf_counter() - start_time_recovered
if solution is None:
results['Recovered_from_primal'] = None
print('###################!RESULTS!#############################')
print('Recovered_from_primal:\nSolution was not found!')
print('__________________________________________________________')
else:
objective, strains, routes, time = ( solution['Objective'], solution['Strain'],
solution['Route'], solution['Time'] )
results['Recovered_from_primal'] = { 'Objective': objective, 'Time': time, 'Total time': elapsed_time_recovered }
print('###################!RESULTS!#############################')
print(f'Recovered_from_primal:\nObjective: {objective}, Time: {time}, Total time: {elapsed_time_recovered}')
#validate constraints violations
if validate_feasability:
violations = mp.FindConstraintsViolation(model_original.cmodel, strains, routes)
cc_vn = violations['capacity_constraints'][1]
rc_vn = violations['route_constraints'][1]
if cc_vn == 0 and rc_vn == 0:
print('Feasible')
else:
print(f'Capacity constraint violations number: {cc_vn}')
print(f'Route constraint violations number: {rc_vn}')
print('__________________________________________________________')
#draw if needed
if draw_solution:
path_list_sol = [ [edge for edge, value in route.items() if abs(value - 1) <= 1e-6] for route in routes ]
flow_list_sol = strains
network_graph.SetPath(*path_list_sol)
network_graph.SetFlows(*flow_list_sol)
dr.PlotNetwork(network_graph, 'ORIGINAL')
#show figures if needed
if draw_solution or draw_progress:
dr.ShowAll()
return results
# {'NodesNumber': 4, 'EdgesNumber': 12, 'ExternalEdgesNumber': 2 }
# ]
# network = net.NetworkGraph.GenerateSmallWorld(world_init_data, *CAPACITY_BOUNDS)
#network.GenerateRandomSrcDst(12)
# #Network large
# world_init_data = [ {'NodesNumber': 3, 'EdgesNumber': 8, 'ExternalEdgesNumber': 3 },
# {'NodesNumber': 12, 'EdgesNumber': 32, 'ExternalEdgesNumber': 3 },
# {'NodesNumber': 6, 'EdgesNumber': 22, 'ExternalEdgesNumber': 2 },
# {'NodesNumber': 4, 'EdgesNumber': 12, 'ExternalEdgesNumber': 2 },
# {'NodesNumber': 9, 'EdgesNumber': 16, 'ExternalEdgesNumber': 1 },
# {'NodesNumber': 8, 'EdgesNumber': 20, 'ExternalEdgesNumber': 2 },
# {'NodesNumber': 4, 'EdgesNumber': 12, 'ExternalEdgesNumber': 1 },
# {'NodesNumber': 3, 'EdgesNumber': 6, 'ExternalEdgesNumber': 1 },
# ]
# network = net.NetworkGraph.GenerateSmallWorld(world_init_data, *CAPACITY_BOUNDS)
# network.GenerateRandomSrcDst(32)
# test graph getter funcitions
n_list = network.GetNodeList()
f_list = network.GetFlowList()
sd_dict = network.GetSrcDstDict()
a_list = network.GetArcList()
c_dict = network.GetCapacityParam()
# plot network
dr.PlotNetwork(network, 'Test graph')
#show all drawn
dr.ShowAll()
