#----------------------------------------------------------------------------- # Solve the model and display the result #----------------------------------------------------------------------------- # Solve model print("Solving model....") msol = mdl.solve(FailLimit=30000, TimeLimit=10) print("Solution: ") msol.print_solution() if msol and visu.is_visu_enabled(): load = [CpoStepFunction() for j in range(NB_RENEWABLE)] for m in MODES: itv = msol.get_var_solution(modes[m['id']]) if itv.is_present(): for j in range(NB_RENEWABLE): dem = m['demandRenewable'][j] if dem > 0: load[j].add_value(itv.get_start(), itv.get_end(), dem) visu.timeline("Solution for RCPSPMM " + filename) visu.panel("Tasks") for t in TASKS: tid = t['id'] visu.interval(msol.get_var_solution(tasks[tid]), tid, str(tid)) for j in range(NB_RENEWABLE): visu.panel("R " + str(j + 1)) visu.function(segments=[(INTERVAL_MIN, INTERVAL_MAX, CAPACITIES_RENEWABLE[j])], style='area', color='lightgrey') visu.function(segments=load[j], style='area', color=j) visu.show()
#----------------------------------------------------------------------------- print("Solving model....") msol = mdl.solve(TimeLimit=10) print("Solution: ") msol.print_solution() if msol and visu.is_visu_enabled(): visu.timeline("Solution SchedTime", origin=10, horizon=120) visu.panel("Schedule") for t in ALL_TASKS: visu.interval(msol.get_var_solution(tasks[t.id]), t.id, t.name) for t in ALL_TASKS: if t.release_date is not None: visu.panel('Earliness') itvsol = msol.get_var_solution(tasks[t.id]) cost = fearliness[t].get_value(itvsol.get_start()) visu.function(segments=[(itvsol, cost, t.name)], color=t.id, style='interval') visu.function(segments=fearliness[t], color=t.id) if t.due_date is not None: visu.panel('Tardiness') itvsol = msol.get_var_solution(tasks[t.id]) cost = ftardiness[t].get_value(itvsol.get_end()) visu.function(segments=[(itvsol, cost, t.name)], color=t.id, style='interval') visu.function(segments=ftardiness[t], color=t.id) visu.show()
# Solve model print("Solving model....") msol = mdl.solve(FailLimit=30000, TimeLimit=10) print("Solution: ") msol.print_solution() ############################################################################## # Display result ############################################################################## if msol and visu.is_visu_enabled(): load = [CpoStepFunction() for j in range(nb_renewable)] for m in modes_data: itv = msol.get_var_solution(modes[m]) if itv.is_present(): for j in range(nb_renewable): if 0 < m.demand_renewable[j]: load[j].add_value(itv.get_start(), itv.get_end(), m.demand_renewable[j]) visu.timeline("Solution for RCPSPMM " + filename) visu.panel("Tasks") for t in tasks_data: visu.interval(msol.get_var_solution(tasks[t]), int(t.name[1:]), t.name) for j in range(nb_renewable): visu.panel("R " + str(j + 1)) visu.function(segments=[(INTERVAL_MIN, INTERVAL_MAX, cap_renewables[j])], style='area', color='lightgrey') visu.function(segments=load[j], style='area', color=j) visu.show()
############################################################################## # Solve model print("Solving model....") msol = mdl.solve(FailLimit=100000, TimeLimit=10) print("Solution: ") msol.print_solution() ############################################################################## # Display result ############################################################################## if msol and visu.is_visu_enabled(): load = [CpoStepFunction() for j in range(nbResources)] for i in range(nbTasks): itv = msol.get_var_solution(tasks[i]) for j in range(nbResources): if 0 < demands[i][j]: load[j].add_value(itv.get_start(), itv.get_end(), demands[i][j]) visu.timeline("Solution for RCPSP " + filename) visu.panel("Tasks") for i in range(nbTasks): visu.interval(msol.get_var_solution(tasks[i]), i, tasks[i].get_name()) for j in range(nbResources): visu.panel("R " + str(j + 1)) visu.function(segments=[(INTERVAL_MIN, INTERVAL_MAX, capacities[j])], style='area', color='lightgrey') visu.function(segments=load[j], style='area', color=j) visu.show()
return task[0].upper() + loc # Solve model print("Solving model....") msol = mdl.solve(FailLimit=10000, TimeLimit=10) print("Solution: ") msol.print_solution() # Display result if msol and visu.is_visu_enabled(): workersF = CpoStepFunction() cashF = CpoStepFunction() for p in range(5): cashF.add_value(60 * p, INT_MAX, 30000) for task in all_tasks: itv = msol.get_var_solution(task) workersF.add_value(itv.get_start(), itv.get_end(), 1) cashF.add_value(itv.start, INT_MAX, -200 * desc[task].duration) visu.timeline('Solution SchedCumul') visu.panel(name="Schedule") for task in all_tasks: visu.interval(msol.get_var_solution(task), house[task], compact(task.get_name())) visu.panel(name="Workers") visu.function(segments=workersF, style='area') visu.panel(name="Cash") visu.function(segments=cashF, style='area', color='gold') visu.show()
############################################################################## # Solve model print("Solving model....") msol = mdl.solve(FailLimit=100000, TimeLimit=10) print("Solution: ") msol.print_solution() ############################################################################## # Display result ############################################################################## if msol and visu.is_visu_enabled(): load = [CpoStepFunction() for j in range(nbResources)] for i in range(nbTasks): itv = msol.get_var_solution(tasks[i]) for j in range(nbResources): if 0 < demands[i][j]: load[j].add_value(itv.get_start(), itv.get_end(), demands[i][j]) visu.timeline("Solution for RCPSP " + filename) visu.panel("Tasks") for i in range(nbTasks): visu.interval(msol.get_var_solution(tasks[i]), i, tasks[i].get_name()) for j in range(nbResources): visu.panel("R " + str(j + 1)) visu.function(segments=[(INTERVAL_MIN, INTERVAL_MAX, capacities[j])], style="area", color="lightgrey") visu.function(segments=load[j], style="area", color=j) visu.show()
print("Solution: ") msol.print_solution() ############################################################################## # Display result ############################################################################## def compact(name): # Example: H3-garden -> G3 # ^ ^ loc, task = name[1:].split('-', 1) return task[0].upper() + loc if msol and visu.is_visu_enabled(): workers_function = CpoStepFunction() for v in all_tasks: itv = msol.get_var_solution(v) workers_function.add_value(itv.get_start(), itv.get_end(), 1) visu.timeline('Solution SchedState') visu.panel(name="Schedule") for v in all_tasks: visu.interval(msol.get_var_solution(v), house[v], compact(v.get_name())) visu.panel(name="Houses state") for f in all_state_functions: visu.sequence(name=f.get_name(), segments=msol.get_var_solution(f)) visu.panel(name="Nb of workers") visu.function(segments=workers_function, style='line') visu.show()
val = msol.get_value(interval) if val != (): tasks.append( (msol.get_var_solution(interval), 1, interval.get_name())) cost_sum += energy_intervals_array[val[2] - 1].get_value( val[0]) * task['power_consumption'] # Add segments to cost function for i in range(val[0], val[1]): cost_i = energy_prices[i] * task['power_consumption'] energy_costs.add_value(i, i + 1, cost_i) # Do not show this machine if no task if assigned to it if tasks and ons: visu.timeline("Machine " + str(m_id), 0, int(TIMESLOTS)) visu.panel("Tasks") visu.sequence(name='Machine', intervals=ons) visu.sequence(name='Tasks', intervals=tasks) visu.function(name='Cost={}'.format(cost_sum), segments=energy_costs) for j in range(NUM_RESOURCES): visu.panel('resources_{}'.format(j)) res = CpoStepFunction() for task, interval in task_intervals_on_machines[m_id]: val = msol.get_value(interval) if val != (): res.add_value(val[0], val[1], task['resource_usage'][j]) visu.function(segments=res, color=j) visu.show()
############################################################################## # Display result ############################################################################## def compact(name): # Example: H3-garden -> G3 # ^ ^ loc, task = name[1:].split('-', 1) return task[0].upper() + loc if msol and visu.is_visu_enabled(): workersF = CpoStepFunction() cashF = CpoStepFunction() for p in range(5): cashF.add_value(60 * p, INT_MAX, 30000) for task in all_tasks: itv = msol.get_var_solution(task) workersF.add_value(itv.get_start(), itv.get_end(), 1) cashF.add_value(itv.start, INT_MAX, -200 * desc[task].duration) visu.timeline('Solution SchedCumul') visu.panel(name="Schedule") for task in all_tasks: visu.interval(msol.get_var_solution(task), house[task], compact(task.get_name())) visu.panel(name="Workers") visu.function(segments=workersF, style='area') visu.panel(name="Cash") visu.function(segments=cashF, style='area', color='gold') visu.show()
# 1. Calling the solve print("\nSolving model....") msol = mdl.solve(url=None, key=None, FailLimit=30000) print("done") # 2. Displaying the objective and solution print("Cost will be " + str(msol.get_objective_values()[0])) # 3. Viewing the results of sequencing problems in a Gantt chart rcParams['figure.figsize'] = 15, 3 workers_function = CpoStepFunction() for h in Houses: for t in TaskNames: itv = msol.get_var_solution(task[h, t]) workers_function.add_value(itv.get_start(), itv.get_end(), 1) visu.timeline('Solution SchedState') visu.panel(name="Schedule") for h in Houses: for t in TaskNames: visu.interval(msol.get_var_solution(task[h, t]), h, t) visu.panel(name="Houses state") for h in Houses: f = state[h] visu.sequence(name=f.get_name(), segments=msol.get_var_solution(f)) visu.panel(name="Nb of workers") visu.function(segments=workers_function, style='line') visu.show()
# 전체 인력 사용 수준 그래프 print("workforce level for" + str(team)) rcParams['figure.figsize'] = 15, 3 workforce=CpoStepFunction() workforce_standard=CpoStepFunction() workforce_whole=CpoStepFunction() for key, value in acts.items(): wf=sol.get_var_solution(act_var[key]) workforce.add_value(wf.get_start(),wf.get_end(), acts[key].worker) workforce_whole.add_value(0, last_end, sum(worker)) workforce_standard.add_value(0,last_end ,standard) # 공종 별 인력 사용 수준 그래프 visu.panel(name="Workers") visu.function(segments=workforce, style='area', origin=0, horizon=last_end, color='lightgreen') visu.function(segments=workforce_standard, style='line', origin=0, horizon=last_end) visu.function(segments=workforce_whole, style='line', origin=0, horizon=last_end, color='blue') visu.show(pngfile="overall_workforce") print(worker) for i in range(len(worker)): workforce_bywork = CpoStepFunction() workforce_limit = CpoStepFunction() visu.panel(name="Work " + str(i + 1)) for key, value in acts.items(): if acts[key].process == "W" + str(i + 1): itv = sol.get_var_solution(act_var[key]) workforce_bywork.add_value(itv.get_start(), itv.get_end(), acts[key].worker) workforce_limit.add_value(0, last_end, worker[i])
############################################################################## print("Solving model....") msol = mdl.solve(TimeLimit=10) print("Solution: ") msol.print_solution() ############################################################################## # Display result ############################################################################## if msol and visu.is_visu_enabled(): visu.timeline("Solution SchedTime", origin=10, horizon=120) visu.panel("Schedule") for t in ALL_TASKS: visu.interval(msol.get_var_solution(tasks[t.id]), t.id, t.name) for t in ALL_TASKS: if t.release_date is not None: visu.panel("Earliness") itvsol = msol.get_var_solution(tasks[t.id]) cost = fearliness[t].get_value(itvsol.get_start()) visu.function(segments=[(itvsol, cost, t.name)], color=t.id, style="interval") visu.function(segments=fearliness[t], color=t.id) if t.due_date is not None: visu.panel("Tardiness") itvsol = msol.get_var_solution(tasks[t.id]) cost = ftardiness[t].get_value(itvsol.get_end()) visu.function(segments=[(itvsol, cost, t.name)], color=t.id, style="interval") visu.function(segments=ftardiness[t], color=t.id) visu.show()