* Cute AMPL model (translation to GAMS) * * AMPL Model by Hande Y. Benson * * Copyright (C) 2001 Princeton University * All Rights Reserved * * Permission to use, copy, modify, and distribute this software and * its documentation for any purpose and without fee is hereby * granted, provided that the above copyright notice appear in all * copies and that the copyright notice and this * permission notice appear in all supporting documentation. * Source: * A partial specification of problem SPANHYD. * SIF input: Ph. Toint, Sept 1990. * classification LNR2-MN-97-33 $Set N 97 Set i /i1*i%N%/; parameter l(i) / i1 77.0 , i2 1124.52 , i3 158.0 , i4 16.0 , i5 0.0 , i6 783.65 , i7 11.0 , i8 49.0 , i9 2155.17 , i10 252.0 , i11 0.0 , i12 0.0 , i13 0.0 , i14 0.0 , i15 0.0 , i16 0.0 , i17 0.0 , i18 0.0 , i19 0.0 , i20 0.0 , i21 0.0 , i22 0.0 , i23 0.0 , i24 0.0 , i25 0.0 , i26 0.0 , i27 0.0 , i28 0.0 , i29 0.0 , i30 77.0 , i31 403.4 , i32 158.0 , i33 16.0 , i34 0.0 , i35 502.0 , i36 11.0 , i37 49.0 , i38 915.3 , i39 252.0 , i40 0.0 , i41 0.0 , i42 0.0 , i43 0.0 , i44 0.0 , i45 0.0 , i46 0.0 , i47 0.0 , i48 0.0 , i49 0.0 , i50 0.0 , i51 0.0 , i52 0.0 , i53 0.0 , i54 0.0 , i55 0.0 , i56 0.0 , i57 0.0 , i58 0.0 , i59 77.0 , i60 403.4 , i61 158.0 , i62 16.0 , i63 0.0 , i64 505.64 , i65 11.0 , i66 49.0 , i67 915.3 , i68 252.0 , i69 0.0 , i70 0.0 , i71 0.0 , i72 0.0 , i73 0.0 , i74 0.0 , i75 0.0 , i76 0.0 , i77 0.0 , i78 0.0 , i79 0.0 , i80 0.0 , i81 0.0 , i82 0.0 , i83 0.0 , i84 0.0 , i85 0.0 , i86 0.0 , i87 0.0 , i88 77.0 , i89 1100.0 , i90 158.0 , i91 16.0 , i92 0.0 , i93 700.0 , i94 11.0 , i95 49.0 , i96 2000.0 , i97 252.0 / ; parameter u(i) / i1 77.01 , i2 1124.53 , i3 158.01 , i4 16.01 , i5 0.0 , i6 783.66 , i7 11.01 , i8 49.01 , i9 2155.18 , i10 252.01 , i11 397.84 , i12 222.32 , i13 205.63 , i14 205.63 , i15 205.63 , i16 124.83 , i17 127.01 , i18 61.08 , i19 614.84 , i20 778.08 , i21 3024.0 , i22 3024.0 , i23 3024.0 , i24 3024.0 , i25 7257.6 , i26 1209.6 , i27 907.2 , i28 7257.6 , i29 7257.6 , i30 77.0 , i31 1312.0 , i32 158.0 , i33 16.0 , i34 0.0 , i35 928.46 , i36 11.0 , i37 49.0 , i38 2611.6 , i39 252.0 , i40 397.84 , i41 222.32 , i42 205.63 , i43 205.63 , i44 205.63 , i45 124.83 , i46 127.01 , i47 61.08 , i48 614.84 , i49 778.08 , i50 3024.0 , i51 3024.0 , i52 3024.0 , i53 3024.0 , i54 7257.6 , i55 1209.6 , i56 907.2 , i57 7257.6 , i58 7257.6 , i59 77.0 , i60 1312.0 , i61 158.0 , i62 16.0 , i63 0.0 , i64 928.46 , i65 11.0 , i66 49.0 , i67 2611.6 , i68 252.0 , i69 397.84 , i70 222.32 , i71 205.63 , i72 205.63 , i73 205.63 , i74 124.83 , i75 127.01 , i76 61.08 , i77 614.84 , i78 778.08 , i79 3024.0 , i80 3024.0 , i81 3024.0 , i82 3024.0 , i83 7257.6 , i84 1209.6 , i85 907.2 , i86 7257.6 , i87 7257.6 , i88 77.01 , i89 1100.01 , i90 158.01 , i91 16.01 , i92 0.0 , i93 700.01 , i94 11.01 , i95 49.01 , i96 2000.01 , i97 252.01 / ; parameter x0(i) / i1 77.0 , i2 1124.52 , i3 158.0 , i4 16.0 , i5 0.0 , i6 783.65 , i7 11.0 , i8 49.0 , i9 2155.17 , i10 252.0 , i11 51.38 , i12 140.21 , i13 142.79 , i14 21.91 , i15 164.7 , i16 58.19 , i17 58.19 , i18 61.08 , i19 566.43 , i20 583.57 , i21 0.0 , i22 0.0 , i23 0.0 , i24 0.0 , i25 0.0 , i26 0.0 , i27 0.0 , i28 0.0 , i29 0.0 , i30 77.0 , i31 1049.53 , i32 158.0 , i33 16.0 , i34 0.0 , i35 738.43 , i36 11.0 , i37 49.0 , i38 1835.36 , i39 252.0 , i40 32.06 , i41 0.0 , i42 4.2 , i43 48.37 , i44 52.57 , i45 59.85 , i46 59.85 , i47 58.24 , i48 0.0 , i49 18.76 , i50 0.0 , i51 0.0 , i52 0.0 , i53 0.0 , i54 0.0 , i55 0.0 , i56 0.0 , i57 0.0 , i58 0.0 , i59 77.0 , i60 1081.87 , i61 158.0 , i62 16.0 , i63 0.0 , i64 696.71 , i65 11.0 , i66 49.0 , i67 1972.77 , i68 252.0 , i69 18.13 , i70 0.0 , i71 0.0 , i72 18.13 , i73 18.13 , i74 5.81 , i75 5.81 , i76 0.0 , i77 0.0 , i78 6.02 , i79 0.0 , i80 0.0 , i81 0.0 , i82 0.0 , i83 0.0 , i84 0.0 , i85 0.0 , i86 0.0 , i87 0.0 , i88 77.0 , i89 1100.0 , i90 158.0 , i91 16.0 , i92 0.0 , i93 700.0 , i94 11.0 , i95 49.0 , i96 2000.0 , i97 252.0 / ; Variable x(i) , obj ; Equation n1 , n2 , n3 , n4 , n5 , n6 , n7 , n8 , n9 , n10, n11 , n12 , n13 , n14 , n15 , n16 , n17 , n18 , n19 , n20, n21 , n22 , n23 , n24 , n25 , n26 , n27 , n28 , n29 , n30, n31 , n32 , n33 , Def_obj ; n1.. x['i1'] - x['i11'] - x['i21'] - x['i30'] + 51.38 =e= 0; n2.. x['i2'] + x['i11'] - x['i12'] + x['i21'] - x['i22'] - x['i31'] + 13.84 =e= 0; n3.. x['i3'] + x['i12'] - x['i13'] + x['i22'] - x['i23'] - x['i32'] + 2.58 =e= 0; n4.. x['i4'] - x['i14'] - x['i24'] - x['i33'] + 21.91 =e= 0; n5.. x['i5'] + x['i13'] + x['i14'] - x['i15'] - x['i34'] =e= 0; n6.. x['i6'] - x['i16'] - x['i25'] - x['i35'] + 12.97 =e= 0; n7.. x['i7'] + x['i16'] - x['i17'] + x['i25'] - x['i26'] - x['i36'] =e= 0; n8.. x['i8'] + x['i17'] - x['i18'] + x['i26'] - x['i27'] - x['i37'] + 2.89 =e= 0; n9.. x['i9'] + x['i15'] + x['i18'] - x['i19'] + x['i23'] + x['i24'] + x['i27'] - x['i28'] - x['i38'] +20.84 =e= 0; n10.. x['i10'] + x['i19'] - x['i20'] + x['i28'] - x['i29'] - x['i39'] + 17.14 =e= 0; n11.. x['i30'] - x['i40'] - x['i50'] - x['i59'] + 32.06 =e= 0; n12.. x['i31'] + x['i40'] - x['i41'] + x['i50'] - x['i51'] - x['i60'] + 0.28 =e= 0; n13.. x['i32'] + x['i41'] - x['i42'] + x['i51'] - x['i52'] - x['i61'] + 4.2 =e= 0; n14.. x['i33'] - x['i43'] - x['i53'] - x['i62'] + 48.37 =e= 0; n15.. x['i34'] + x['i42'] + x['i43'] - x['i44'] - x['i63'] =e= 0; n16.. x['i35'] - x['i45'] - x['i54'] - x['i64'] + 18.13 =e= 0; n17.. x['i36'] + x['i45'] - x['i46'] + x['i54'] - x['i55'] - x['i65'] =e= 0; n18.. x['i37'] + x['i46'] - x['i47'] + x['i55'] - x['i56'] - x['i66'] - 1.61 =e= 0; n19.. x['i38'] + x['i44'] + x['i47'] - x['i48'] + x['i52'] + x['i53'] + x['i56'] - x['i57'] - x['i67'] +26.6 =e= 0; n20.. x['i39'] + x['i48'] - x['i49'] + x['i57'] - x['i58'] - x['i68'] + 18.76 =e= 0; n21.. x['i59'] - x['i69'] - x['i79'] - x['i88'] + 18.13 =e= 0; n22.. x['i60'] + x['i69'] - x['i70'] + x['i79'] - x['i80'] - x['i89'] =e= 0; n23.. x['i61'] + x['i70'] - x['i71'] + x['i80'] - x['i81'] - x['i90'] =e= 0; n24.. x['i62'] - x['i72'] - x['i82'] - x['i91'] + 18.13 =e= 0; n25.. x['i63'] + x['i71'] + x['i72'] - x['i73'] - x['i92'] =e= 0; n26.. x['i64'] - x['i74'] - x['i83'] - x['i93'] + 9.1 =e= 0; n27.. x['i65'] + x['i74'] - x['i75'] + x['i83'] - x['i84'] - x['i94'] =e= 0; n28.. x['i66'] + x['i75'] - x['i76'] + x['i84'] - x['i85'] - x['i95'] - 5.81 =e= 0; n29.. x['i67'] + x['i73'] + x['i76'] - x['i77'] + x['i81'] + x['i82'] + x['i85'] - x['i86'] - x['i96'] +9.1 =e= 0; n30.. x['i68'] + x['i77'] - x['i78'] + x['i86'] - x['i87'] - x['i97'] + 6.02 =e= 0; n31.. x['i20'] + x['i29'] + x['i49'] + x['i58'] + x['i78'] + x['i87'] - 608.35 =e= 0; n32.. -x['i1'] - x['i2'] - x['i3'] - x['i4'] - x['i5'] - x['i6'] - x['i7'] - x['i8'] - x['i9'] - x['i10'] +4626.34 =e= 0; n33.. x['i88'] + x['i89'] + x['i90'] + x['i91'] + x['i92'] + x['i93'] + x['i94'] + x['i95'] + x['i96'] +x['i97'] - 4363.0 =e= 0; Def_obj.. obj =e= - x['i1']; x.lo(i) = l(i) ; x.up(i) = u(i) ; x.l(i) = x0(i) ; Model linspanh /all/; Solve linspanh using nlp minimize obj; display x.l ; display obj.l ;