$ontext gppC3 pooling problem data. Author: Mohammed Alfaki, Wed Nov 9 15:36:09 2011 $offtext $eolcom # # Declare sets set i / 1*20 /; set s(i) / 1*8 /; set t(i) /15*20 /; set k / 1*4 /; alias (i,j); # The arc unit cost c_{ij} table c(i,j) 9 10 11 12 13 14 15 16 17 18 19 20 1 3.00 0.00 0.00 0.00 3.00 3.00 0.00 0.00 0.00 -4.00 -7.00 0.00 2 4.00 4.00 0.00 4.00 4.00 0.00 -3.00 0.00 -9.00 -3.00 0.00 0.00 3 5.00 5.00 0.00 5.00 0.00 5.00 0.00 0.00 0.00 0.00 0.00 -5.00 4 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 5 2.00 0.00 2.00 2.00 2.00 2.00 0.00 0.00 0.00 0.00 0.00 0.00 6 0.00 3.00 0.00 3.00 0.00 0.00 0.00 0.00 -10.00 0.00 -7.00 0.00 7 1.00 1.00 0.00 1.00 0.00 0.00 -6.00 0.00 -12.00 0.00 0.00 0.00 8 0.00 0.00 2.00 0.00 0.00 2.00 0.00 0.00 0.00 0.00 -8.00 -8.00 9 0.00 0.00 0.00 0.00 0.00 0.00 -7.00 -14.00 -13.00 -7.00 -10.00 0.00 10 0.00 0.00 0.00 0.00 0.00 0.00 -7.00 0.00 -13.00 -7.00 -10.00 -10.00 11 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 -13.00 -7.00 -10.00 0.00 12 0.00 0.00 0.00 0.00 0.00 0.00 0.00 -14.00 0.00 -7.00 -10.00 -10.00 13 0.00 0.00 0.00 0.00 0.00 0.00 -7.00 -14.00 -13.00 0.00 0.00 -10.00 14 0.00 0.00 0.00 0.00 0.00 0.00 -7.00 0.00 -13.00 -7.00 0.00 -10.00 ; # The adjacency matrix (the arcs set A) table a(i,j) 9 10 11 12 13 14 15 16 17 18 19 20 1 1 0 0 0 1 1 0 0 0 1 1 0 2 1 1 0 1 1 0 1 0 1 1 0 0 3 1 1 0 1 0 1 0 0 0 0 0 1 4 1 0 1 1 1 0 0 0 0 0 0 0 5 1 0 1 1 1 1 0 0 0 0 0 0 6 0 1 0 1 0 0 0 0 1 0 1 0 7 1 1 0 1 0 0 1 0 1 0 0 0 8 0 0 1 0 0 1 0 0 0 0 1 1 9 0 0 1 1 1 1 1 1 1 1 1 0 10 0 0 1 0 1 1 1 0 1 1 1 1 11 0 0 0 1 1 1 0 0 1 1 1 0 12 0 0 0 0 0 1 0 1 0 1 1 1 13 0 0 0 0 0 0 1 1 1 0 0 1 14 0 0 0 0 0 0 1 0 1 1 0 1 ; # Source qualities/terminal quality upper bounds table q(i,k) 1 2 3 4 1 4.00 7.00 0.00 4.00 2 9.00 0.00 4.00 7.00 3 2.00 2.00 5.00 7.00 4 9.00 1.00 9.00 2.00 5 8.00 1.00 0.00 4.00 6 1.00 6.00 0.00 4.00 7 4.00 5.00 0.00 2.00 8 6.00 1.00 0.00 8.00 15 4.00 6.00 4.00 6.00 16 6.00 2.00 6.00 6.00 17 5.00 5.00 2.00 2.00 18 6.00 6.00 5.00 4.00 19 3.00 6.00 4.00 5.00 20 4.00 4.00 6.00 3.00 ; # Node capacity lower bound parameter bl(i) / 1 0.00 2 0.00 3 0.00 4 0.00 5 0.00 6 0.00 7 0.00 8 0.00 15 0.00 16 0.00 17 0.00 18 0.00 19 0.00 20 0.00 / ; # Node capacity upper bound parameter bu(i) / 1 22.00 2 53.00 3 27.00 4 57.00 5 35.00 6 56.00 7 23.00 8 23.00 9 40.00 10 57.00 11 26.00 12 51.00 13 26.00 14 53.00 15 42.00 16 59.00 17 36.00 18 52.00 19 22.00 20 58.00 / ; $include xmodel.gms