reset;

option solver cplexamp;
model MTC3_fix_y.mod;  						# load the optimization model

#data 20100311_MTC3_10j_v_lj.dat;   					# load the model data
data 20100504_MTC3_10j.dat

problem MTC3_mach: t_mach, z_mach, y_mach, s_mach, h_mach, Finish_times_and_tardiness_mach, opconstraints_mach,
					flexconstraints_mach, both_y_not_1_constraints_mach, at_least_one_y_constraints_mach,
					earliest_start_constraints_mach, start_constraints1_mach,
					start_constraints2_mach, completion_time_constraints_mach, interop_time_constraints_mach,
					tardiness_constraints_mach;
option cplex_options
	'mipdisplay=2'
	'timing=1'
	'mipinterval=2000'
#	'mipgap=0.1'
;

problem MTC3_feas: t, z, y, s, h, Process_times_and_tardiness_feas, opconstraints_feas,
					flexconstraints, both_y_not_1_constraints, at_least_one_y_constraints,
					earliest_start_constraints, same_job_constraints, start_constraints1,
					start_constraints2, completion_time_constraints, interop_time_constraints,
#					t2_fixed_constraints_feas,					
					y_fixed_constraints_feas, z2_fixed_constraints_feas, tardiness_constraints;

option cplex_options
	'mipdisplay=2'
	'timing=1'
	'mipinterval=2000';
option presolve_eps 1.0e-10;


solve MTC3_mach;          		
option omit_zero_rows 1;
print "";
print "----------Modell MTC3_mach - Only machining processing operations--------------";
print "----------z_mach_jk (which job j is scheduled to which resource k)-------------";

display z_mach;
 
option omit_zero_rows 0, display_transpose 10;
print "----------Starting time t_mach_j (hours from time 0)---------------------------";
display t_mach; 
print "----------Completion time s_j (hours)------------------------------------------";
display s_mach; 
print "----------Tardiness h_j (hours)------------------------------------------------";
display h_mach; 

let {j in JOBS, q in JOBS, k in K_RESOURCES} y_mach_solution[j,q,k] := 0;	#so y_mach_solution has a value for the non-machining resources
let {j in JOBS, q in JOBS, k in K_mach_RESOURCES} y_mach_solution[j,q,k] := y_mach[j,q,k];

let {j in JOBS, k in K_RESOURCES} z_mach_solution[j,k] := 0;	#so z_mach_solution has a value for the non-machining resources
let {j in JOBS, k in K_mach_RESOURCES} z_mach_solution[j,k] := z_mach[j,k];

solve MTC3_feas;

print "";
print "----------Modell MTC3_feas - The feasibility problem---------------------------";
print "----------z_ijk (which operation (i,j) is scheduled to which resource k)-------";

option omit_zero_rows 0, display_transpose -10;
display z;
 
option omit_zero_rows 0, display_transpose 10;
print "----------Starting time t_ij (If (tr) jobs j in rows,(hours)-------------------";
display t; 
print "----------Completion time s_j (hours)------------------------------------------";
display s; 

print "----------Tardiness h_j (hours)------------------------------------------------";
display h;

let proc_tardi_objective := sum{j in JOBS} (s[j]-t[1,j]+h[j]);

print "---Sum of processing times and tardiness for the solution of the feasibility problem---";
display proc_tardi_objective;

  
#display y_mach;
#display z_mach;