$ontext gppC2 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 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 -6.00 -8.00 0.00 0.00 2 4.00 0.00 4.00 4.00 4.00 0.00 0.00 0.00 -2.00 0.00 -5.00 -5.00 3 0.00 4.00 4.00 4.00 4.00 4.00 0.00 0.00 0.00 0.00 0.00 -5.00 4 4.00 0.00 0.00 4.00 0.00 4.00 0.00 0.00 0.00 0.00 -5.00 0.00 5 2.00 2.00 0.00 2.00 2.00 0.00 0.00 0.00 0.00 -6.00 -7.00 0.00 6 4.00 4.00 0.00 4.00 4.00 0.00 0.00 -9.00 0.00 0.00 0.00 0.00 7 1.00 0.00 0.00 0.00 0.00 1.00 0.00 -12.00 0.00 -7.00 0.00 0.00 8 0.00 3.00 0.00 3.00 3.00 3.00 0.00 0.00 0.00 0.00 0.00 -6.00 9 0.00 0.00 0.00 0.00 0.00 0.00 -13.00 -13.00 -6.00 0.00 0.00 -9.00 10 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 -6.00 0.00 0.00 -9.00 11 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 -6.00 -8.00 -9.00 -9.00 12 0.00 0.00 0.00 0.00 0.00 0.00 -13.00 0.00 0.00 0.00 -9.00 0.00 13 0.00 0.00 0.00 0.00 0.00 0.00 -13.00 0.00 0.00 -8.00 -9.00 -9.00 14 0.00 0.00 0.00 0.00 0.00 0.00 -13.00 -13.00 -6.00 -8.00 0.00 0.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 1 0 0 1 0 0 0 1 1 0 0 2 1 0 1 1 1 0 0 0 1 0 1 1 3 0 1 1 1 1 1 0 0 0 0 0 1 4 1 0 0 1 0 1 0 0 0 0 1 0 5 1 1 0 1 1 0 0 0 0 1 1 0 6 1 1 0 1 1 0 0 1 0 0 0 0 7 1 0 0 0 0 1 0 1 0 1 0 0 8 0 1 0 1 1 1 0 0 0 0 0 1 9 0 0 0 1 1 1 1 1 1 0 0 1 10 0 0 1 1 0 0 0 0 1 0 0 1 11 0 0 0 1 0 1 0 0 1 1 1 1 12 0 0 0 0 0 1 1 0 0 0 1 0 13 0 0 0 0 0 1 1 0 0 1 1 1 14 0 0 0 0 0 0 1 1 1 1 0 0 ; # Source qualities/terminal quality upper bounds table q(i,k) 1 2 3 4 1 7.00 4.00 0.00 3.00 2 6.00 9.00 9.00 6.00 3 4.00 3.00 8.00 5.00 4 1.00 6.00 6.00 2.00 5 7.00 9.00 5.00 2.00 6 2.00 2.00 8.00 7.00 7 8.00 0.00 5.00 0.00 8 9.00 4.00 6.00 5.00 15 3.00 4.00 5.00 5.00 16 5.00 4.00 5.00 5.00 17 5.00 3.00 5.00 2.00 18 6.00 2.00 2.00 5.00 19 5.00 2.00 3.00 4.00 20 5.00 3.00 3.00 2.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 48.00 2 34.00 3 27.00 4 39.00 5 42.00 6 22.00 7 44.00 8 46.00 9 37.00 10 24.00 11 43.00 12 24.00 13 20.00 14 20.00 15 27.00 16 44.00 17 45.00 18 26.00 19 58.00 20 47.00 / ; $include xmodel.gms