MINLP World [ MINLP World Home | Board | Solvers | MINLPLib | Links | GamsWorld Group | Search | Contact ]

ghg_1veh.gms:

Reference:

Point: p1 
Best known point (p1): Solution value 7.78 (global optimum, BARON and LINDOGLOBAL certificate)


$offlisting * MINLP written by GAMS Convert at 08/05/10 06:59:15 * * Equation counts * Total E G L N X C * 38 18 9 11 0 0 0 * * Variable counts * x b i s1s s2s sc si * Total cont binary integer sos1 sos2 scont sint * 30 18 12 0 0 0 0 0 * FX 1 1 0 0 0 0 0 0 * * Nonzero counts * Total const NL DLL * 131 40 91 0 * * Solve m using MINLP minimizing objvar; * Variables b1,b2,b3,b4,b5,b6,b7,b8,b9,b10,b11,b12,x13,x14,x15,x16,x17,x18,x19 ,x20,x21,x22,x23,x24,x25,x26,x27,x28,x29,objvar; Positive Variables x24,x25,x26,x27,x28,x29; Binary Variables b1,b2,b3,b4,b5,b6,b7,b8,b9,b10,b11,b12; Equations e1,e2,e3,e4,e5,e6,e7,e8,e9,e10,e11,e12,e13,e14,e15,e16,e17,e18,e19 ,e20,e21,e22,e23,e24,e25,e26,e27,e28,e29,e30,e31,e32,e33,e34,e35,e36 ,e37,e38; e1.. - 21.6*x15*x16*x20 + x24 =E= 0; e2.. - (0.007852585706*x13**3 + 0.154288922601*x14**3 + 0.352933730854*x15**3 - 0.004816150342*x13**2*x14 - 0.00547943134*x14**2*x13 - 0.02533808214* x13**2*x15 + 0.00021201136*x15**2*x13 - 0.057118497613*x14**2*x15 - 0.042739509965*x15**2*x14 - 0.01583097252*x13*x14*x15 - 0.001028174658*x13 **2 - 0.805369774847*x14**2 - 0.655580550098*x15**2 + 0.057270405947*x13* x14 + 0.07973036236*x13*x15 + 0.342091579946*x14*x15 - 0.191345333621*x13 + 1.188971392024*x14 - 0.346682012779*x15) + x20 =E= 4.960068215723; e3.. - (2.21406746341*x13**3 + 1.086659693552*x14**3 + 5.577874978662*x15**3 - 0.815241697738*x13**2*x14 + 0.509578110533*x14**2*x13 + 1.561758113326* x13**2*x15 + 2.212321055022*x15**2*x13 - 0.612567680918*x14**2*x15 + 0.254008083604*x15**2*x14 - 0.159429747244*x13*x14*x15 - 8.905599398536* x13**2 - 6.095001164559*x14**2 - 15.207539664993*x15**2 + 0.089172114876* x13*x14 - 3.273526677614*x13*x15 + 2.498376358946*x14*x15 + 2.621894664006 *x13 + 9.284846067558*x14 + 5.837143728557*x15) + x21 =E= 57.679680208231; e4.. - (1.456640469666*x13**3 - 5.495718264905*x14**3 - 28.456261951645*x15**3 + 0.912917970314*x13**2*x14 - 0.88119920631*x14**2*x13 - 1.049763024383* x13**2*x15 - 0.308107344863*x15**2*x13 + 2.043536297441*x14**2*x15 + 15.609611231641*x15**2*x14 + 0.336486837518*x13*x14*x15 - 4.634160849469* x13**2 + 31.478262635483*x14**2 + 34.016843490037*x15**2 + 1.153148892739* x13*x14 + 1.168601192983*x13*x15 - 32.056936006397*x14*x15 + 3.405095041238*x13 - 54.472915571467*x14 + 9.56987912824*x15) + x17 =E= 44.230616625681; e5.. - (3.334445194766*x13**3 - 2.265666208775*x14**3 - 20.256566414583*x15**3 + 0.413782262402*x13**2*x14 - 3.523622273943*x14**2*x13 - 0.285910055687* x13**2*x15 - 10.110726634622*x15**2*x13 + 1.95072196814*x14**2*x15 + 10.308512463418*x15**2*x14 + 5.808426325827*x13*x14*x15 - 6.932398033967* x13**2 + 15.80019426934*x14**2 + 39.197963873266*x15**2 + 7.900706395772* x13*x14 + 6.58186092156*x13*x15 - 30.119438887106*x14*x15 - 6.733798415788 *x13 - 26.385308892431*x14 - 4.098268423019*x15) + x18 =E= 32.102172356117; e6.. - (-0.194075741585*x13**2 - 0.004843420334*x14**2 + 0.04736686635*x15**2 + 9.4029979e-5*x13*x14 + 0.011329651793*x13*x15 - 0.001017352942*x14*x15 + 0.382275988592*x13 + 0.019484588535*x14 - 0.077357069039*x15) + x19 =E= 0.140278656706; e7.. x17 =L= 11; e8.. x18 =L= 11; e9.. x19 =G= 0.32; e10.. exp(-0.029595*x24)*(33.7894914681534 + x24) + x26 =E= 33.7894914681534; e11.. exp(-0.029595*x24) + x27 =E= 1; e12.. - 0.134723681728774*(0.010073140669*x13**2 + 0.011394190823*x14**2 + 0.052910213683*x15**2 + 0.000159410872*x13*x14 + 0.008036404292*x13*x15 - 0.003423392047*x14*x15 + 0.097124049148*x13 + 0.03829180344*x14 + 0.370440556286*x15) + x22 =E= 0.29587368369345; e13.. - 0.134723681728774*(0.46598008632*x13**2 - 0.00797004615*x14**2 - 0.01779288613*x15**2 - 0.01429434551*x13*x14 - 0.03832188467*x13*x15 + 0.00970510229*x14*x15 - 0.88981702163*x13 + 0.07730602595*x14 + 0.39988032723*x15) + x23 =E= 0.194162178290626; e14.. - 150000*(0.0181052631578947/x20 + 0.03458*x22 - 0.03458*x23)*x24/(4320* x15 - 5187*x23) + x25 =E= 0; e15.. exp(-0.029595*x25)*(33.7894914681534 + x25) + x28 =E= 33.7894914681534; e16.. exp(-0.029595*x25) + x29 =E= 1; e17.. b1 + b2 + b3 =E= 1; e18.. b1*x24 =L= 0; e19.. b2*x24 =G= 0; e20.. b2*(-200 + x24) =L= 0; e21.. b3*(-200 + x24) =G= 0; e22.. b4 + b5 + b6 =E= 1; e23.. b8*b4*x25 =L= 0; e24.. b8*b5*x25 =G= 0; e25.. b8*b5*(-200 + x25) =L= 0; e26.. b8*b6*(-200 + x25) =G= 0; e27.. b7 + b8 + b9 =E= 1; e28.. (-150000 + 124927.703875072*x15/x23)*b7 =L= 0; e29.. (-150000 + 124927.703875072*x15/x23)*b8 =G= 0; e30.. (150000 - 4320*x15/(0.0172/x20 + 0.03458*x22))*b8 =G= 0; e31.. (150000 - 4320*x15/(0.0172/x20 + 0.03458*x22))*b9 =L= 0; e32.. b7*(-1 + b4) =G= 0; e33.. b9*(-1 + b4) =G= 0; e34.. b2 + b4 + b8 =L= 2; e35.. b3 + b4 + b8 =L= 2; e36.. b3 + b5 + b8 =L= 2; e37.. - (1.87912853526074 + (376.046780997472/x21 + 0.997312113279821*( 0.854659090909091/x20 - 11.34/x21)*x24)*b1 + (0.854659090909091*x26/x20 + 11.34*(33.1610917987189 - x26)/x21 + (0.854659090909091/x20 - 11.34/ x21)*x24*(0.997312113279821 - x27))*b2 + 28.341428570246*b3/x20 + 0.573023666281862*(b4*b8 + b9)*x15 + b1*b5*b8*(0.6*(0.03458*x23*x28 + ( 0.0181052631578947/x20 + 0.03458*x22 - 0.03458*x23)*x24*x29) + 0.01728* x15*(33.1610917987189 - x28)) + b2*b5*b8*(0.6*((0.0181052631578947/x20 + 0.03458*x22)*x26 + 0.03458*x23*(x28 - x26) + (0.0181052631578947/x20 + 0.03458*x22 - 0.03458*x23)*x24*(x29 - x27)) + 0.01728*x15*( 33.1610917987189 - x28)) + 0.6*(b1*b6*b8 + b1*b7)*(0.997312113279821*( 0.0181052631578947/x20 + 0.03458*x22 - 0.03458*x23)*x24 + 1.1467105543997 *x23) + 0.6*(b2*b6*b8 + b2*b7)*((0.0181052631578947/x20 + 0.03458*x22)* x26 + 0.03458*x23*(33.1610917987189 - x26) + (0.0181052631578947/x20 + 0.03458*x22 - 0.03458*x23)*x24*(0.997312113279821 - x27)) + 19.8966550792313*(b3*b6*b8 + b3*b7)*(0.0181052631578947/x20 + 0.03458*x22 ))*b10 - 8.20275610163388*b11 - 14.6264770436496*b12 + objvar =E= 0; e38.. b10 + b11 + b12 =E= 1; * set non default bounds x13.lo = 0.526315789473684; x13.up = 1.05263157894737; x14.lo = 0.961538461538462; x14.up = 2.11538461538462; x15.lo = 0.2; x15.up = 1; x16.fx = 0.8; x17.lo = 6; x17.up = 13; x18.lo = 6; x18.up = 13; x19.lo = 0.26; x19.up = 0.35; x20.lo = 4.9; x20.up = 5.5; x21.lo = 55; x21.up = 63; x22.lo = 0.296392099803303; x22.up = 0.404171045186323; x23.lo = 0.134723681728774; x23.up = 0.229030258938916; x24.up = 90; x25.up = 200; x26.up = 26; x27.up = 1; x28.up = 34.1; x29.up = 1; $if set nostart $goto modeldef * set non default levels b2.l = 1; b5.l = 1; b8.l = 1; b10.l = 1; x13.l = 1; x14.l = 1; x15.l = 0.5; x17.l = 13; x18.l = 10.9460692020431; x19.l = 0.3215334333865; x20.l = 5.218428550001; x21.l = 58.1078648496005; x22.l = 0.344077403769737; x23.l = 0.16888643257787; x24.l = 45.0872226720086; x25.l = 49.2705196815703; x26.l = 13.01907481523; x27.l = 0.736672389572227; x28.l = 14.4644383631733; x29.l = 0.767336256792154; * set non default marginals $label modeldef Model m / all /; m.limrow=1; m.limcol=0; $if NOT '%gams.u1%' == '' $include '%gams.u1%' $if not set MINLP $set MINLP MINLP Solve m using %MINLP% minimizing objvar;