ljcluster.gms
<![CDATA[
* NLP written by GAMS Convert at 06/26/06 16:29:37
*
* Equation counts
* Total E G L N X C
* 46 46 0 0 0 0 0
*
* Variable counts
* x b i s1s s2s sc si
* Total cont binary integer sos1 sos2 scont sint
* 66 66 0 0 0 0 0 0
* FX 0 0 0 0 0 0 0 0
*
* Nonzero counts
* Total const NL DLL
* 271 46 225 0
*
* Solve m using NLP minimizing objvar;
Variables x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,x13,x14,x15,x16,x17,x18,x19
,x20,x21,x22,x23,x24,x25,x26,x27,x28,x29,x30,x31,x32,x33,x34,x35,x36
,x37,x38,x39,x40,x41,x42,x43,x44,x45,x46,x47,x48,x49,x50,x51,x52,x53
,x54,x55,x56,x57,x58,x59,x60,x61,x62,x63,x64,x65,objvar;
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,e39,e40,e41,e42,e43,e44,e45,e46;
e1.. - sqrt(sqr(x3 - x1) + sqr(x4 - x2)) + x21 =E= 0;
e2.. - sqrt(sqr(x5 - x1) + sqr(x6 - x2)) + x22 =E= 0;
e3.. - sqrt(sqr(x5 - x3) + sqr(x6 - x4)) + x23 =E= 0;
e4.. - sqrt(sqr(x7 - x1) + sqr(x8 - x2)) + x24 =E= 0;
e5.. - sqrt(sqr(x7 - x3) + sqr(x8 - x4)) + x25 =E= 0;
e6.. - sqrt(sqr(x7 - x5) + sqr(x8 - x6)) + x26 =E= 0;
e7.. - sqrt(sqr(x9 - x1) + sqr(x10 - x2)) + x27 =E= 0;
e8.. - sqrt(sqr(x9 - x3) + sqr(x10 - x4)) + x28 =E= 0;
e9.. - sqrt(sqr(x9 - x5) + sqr(x10 - x6)) + x29 =E= 0;
e10.. - sqrt(sqr(x9 - x7) + sqr(x10 - x8)) + x30 =E= 0;
e11.. - sqrt(sqr(x11 - x1) + sqr(x12 - x2)) + x31 =E= 0;
e12.. - sqrt(sqr(x11 - x3) + sqr(x12 - x4)) + x32 =E= 0;
e13.. - sqrt(sqr(x11 - x5) + sqr(x12 - x6)) + x33 =E= 0;
e14.. - sqrt(sqr(x11 - x7) + sqr(x12 - x8)) + x34 =E= 0;
e15.. - sqrt(sqr(x11 - x9) + sqr(x12 - x10)) + x35 =E= 0;
e16.. - sqrt(sqr(x13 - x1) + sqr(x14 - x2)) + x36 =E= 0;
e17.. - sqrt(sqr(x13 - x3) + sqr(x14 - x4)) + x37 =E= 0;
e18.. - sqrt(sqr(x13 - x5) + sqr(x14 - x6)) + x38 =E= 0;
e19.. - sqrt(sqr(x13 - x7) + sqr(x14 - x8)) + x39 =E= 0;
e20.. - sqrt(sqr(x13 - x9) + sqr(x14 - x10)) + x40 =E= 0;
e21.. - sqrt(sqr(x13 - x11) + sqr(x14 - x12)) + x41 =E= 0;
e22.. - sqrt(sqr(x15 - x1) + sqr(x16 - x2)) + x42 =E= 0;
e23.. - sqrt(sqr(x15 - x3) + sqr(x16 - x4)) + x43 =E= 0;
e24.. - sqrt(sqr(x15 - x5) + sqr(x16 - x6)) + x44 =E= 0;
e25.. - sqrt(sqr(x15 - x7) + sqr(x16 - x8)) + x45 =E= 0;
e26.. - sqrt(sqr(x15 - x9) + sqr(x16 - x10)) + x46 =E= 0;
e27.. - sqrt(sqr(x15 - x11) + sqr(x16 - x12)) + x47 =E= 0;
e28.. - sqrt(sqr(x15 - x13) + sqr(x16 - x14)) + x48 =E= 0;
e29.. - sqrt(sqr(x17 - x1) + sqr(x18 - x2)) + x49 =E= 0;
e30.. - sqrt(sqr(x17 - x3) + sqr(x18 - x4)) + x50 =E= 0;
e31.. - sqrt(sqr(x17 - x5) + sqr(x18 - x6)) + x51 =E= 0;
e32.. - sqrt(sqr(x17 - x7) + sqr(x18 - x8)) + x52 =E= 0;
e33.. - sqrt(sqr(x17 - x9) + sqr(x18 - x10)) + x53 =E= 0;
e34.. - sqrt(sqr(x17 - x11) + sqr(x18 - x12)) + x54 =E= 0;
e35.. - sqrt(sqr(x17 - x13) + sqr(x18 - x14)) + x55 =E= 0;
e36.. - sqrt(sqr(x17 - x15) + sqr(x18 - x16)) + x56 =E= 0;
e37.. - sqrt(sqr(x19 - x1) + sqr(x20 - x2)) + x57 =E= 0;
e38.. - sqrt(sqr(x19 - x3) + sqr(x20 - x4)) + x58 =E= 0;
e39.. - sqrt(sqr(x19 - x5) + sqr(x20 - x6)) + x59 =E= 0;
e40.. - sqrt(sqr(x19 - x7) + sqr(x20 - x8)) + x60 =E= 0;
e41.. - sqrt(sqr(x19 - x9) + sqr(x20 - x10)) + x61 =E= 0;
e42.. - sqrt(sqr(x19 - x11) + sqr(x20 - x12)) + x62 =E= 0;
e43.. - sqrt(sqr(x19 - x13) + sqr(x20 - x14)) + x63 =E= 0;
e44.. - sqrt(sqr(x19 - x15) + sqr(x20 - x16)) + x64 =E= 0;
e45.. - sqrt(sqr(x19 - x17) + sqr(x20 - x18)) + x65 =E= 0;
e46.. - (1/POWER(x21,12) - 2/POWER(x21,6) + 1/POWER(x22,12) - 2/POWER(x22,6)
+ 1/POWER(x23,12) - 2/POWER(x23,6) + 1/POWER(x24,12) - 2/POWER(x24,6) +
1/POWER(x25,12) - 2/POWER(x25,6) + 1/POWER(x26,12) - 2/POWER(x26,6) + 1/
POWER(x27,12) - 2/POWER(x27,6) + 1/POWER(x28,12) - 2/POWER(x28,6) + 1/
POWER(x29,12) - 2/POWER(x29,6) + 1/POWER(x30,12) - 2/POWER(x30,6) + 1/
POWER(x31,12) - 2/POWER(x31,6) + 1/POWER(x32,12) - 2/POWER(x32,6) + 1/
POWER(x33,12) - 2/POWER(x33,6) + 1/POWER(x34,12) - 2/POWER(x34,6) + 1/
POWER(x35,12) - 2/POWER(x35,6) + 1/POWER(x36,12) - 2/POWER(x36,6) + 1/
POWER(x37,12) - 2/POWER(x37,6) + 1/POWER(x38,12) - 2/POWER(x38,6) + 1/
POWER(x39,12) - 2/POWER(x39,6) + 1/POWER(x40,12) - 2/POWER(x40,6) + 1/
POWER(x41,12) - 2/POWER(x41,6) + 1/POWER(x42,12) - 2/POWER(x42,6) + 1/
POWER(x43,12) - 2/POWER(x43,6) + 1/POWER(x44,12) - 2/POWER(x44,6) + 1/
POWER(x45,12) - 2/POWER(x45,6) + 1/POWER(x46,12) - 2/POWER(x46,6) + 1/
POWER(x47,12) - 2/POWER(x47,6) + 1/POWER(x48,12) - 2/POWER(x48,6) + 1/
POWER(x49,12) - 2/POWER(x49,6) + 1/POWER(x50,12) - 2/POWER(x50,6) + 1/
POWER(x51,12) - 2/POWER(x51,6) + 1/POWER(x52,12) - 2/POWER(x52,6) + 1/
POWER(x53,12) - 2/POWER(x53,6) + 1/POWER(x54,12) - 2/POWER(x54,6) + 1/
POWER(x55,12) - 2/POWER(x55,6) + 1/POWER(x56,12) - 2/POWER(x56,6) + 1/
POWER(x57,12) - 2/POWER(x57,6) + 1/POWER(x58,12) - 2/POWER(x58,6) + 1/
POWER(x59,12) - 2/POWER(x59,6) + 1/POWER(x60,12) - 2/POWER(x60,6) + 1/
POWER(x61,12) - 2/POWER(x61,6) + 1/POWER(x62,12) - 2/POWER(x62,6) + 1/
POWER(x63,12) - 2/POWER(x63,6) + 1/POWER(x64,12) - 2/POWER(x64,6) + 1/
POWER(x65,12) - 2/POWER(x65,6)) + objvar =E= 0;
* set non default bounds
x1.lo = 1;
x2.lo = 1;
x3.lo = 1;
x4.lo = 1;
x5.lo = 1;
x6.lo = 1;
x7.lo = 1;
x8.lo = 1;
x9.lo = 1;
x10.lo = 1;
x11.lo = 1;
x12.lo = 1;
x13.lo = 1;
x14.lo = 1;
x15.lo = 1;
x16.lo = 1;
x17.lo = 1;
x18.lo = 1;
x19.lo = 1;
x20.lo = 1;
x21.lo = 0.8;
x22.lo = 0.8;
x23.lo = 0.8;
x24.lo = 0.8;
x25.lo = 0.8;
x26.lo = 0.8;
x27.lo = 0.8;
x28.lo = 0.8;
x29.lo = 0.8;
x30.lo = 0.8;
x31.lo = 0.8;
x32.lo = 0.8;
x33.lo = 0.8;
x34.lo = 0.8;
x35.lo = 0.8;
x36.lo = 0.8;
x37.lo = 0.8;
x38.lo = 0.8;
x39.lo = 0.8;
x40.lo = 0.8;
x41.lo = 0.8;
x42.lo = 0.8;
x43.lo = 0.8;
x44.lo = 0.8;
x45.lo = 0.8;
x46.lo = 0.8;
x47.lo = 0.8;
x48.lo = 0.8;
x49.lo = 0.8;
x50.lo = 0.8;
x51.lo = 0.8;
x52.lo = 0.8;
x53.lo = 0.8;
x54.lo = 0.8;
x55.lo = 0.8;
x56.lo = 0.8;
x57.lo = 0.8;
x58.lo = 0.8;
x59.lo = 0.8;
x60.lo = 0.8;
x61.lo = 0.8;
x62.lo = 0.8;
x63.lo = 0.8;
x64.lo = 0.8;
x65.lo = 0.8;
* set non default levels
x3.l = 2;
x4.l = 2;
x5.l = 3;
x6.l = 3;
x7.l = 4;
x8.l = 4;
x9.l = 5;
x10.l = 5;
x11.l = 6;
x12.l = 6;
x13.l = 7;
x14.l = 7;
x15.l = 8;
x16.l = 8;
x17.l = 9;
x18.l = 9;
x19.l = 10;
x20.l = 10;
* set non default marginals
Model m / all /;
m.limrow=0; m.limcol=0;
$if NOT '%gams.u1%' == '' $include '%gams.u1%'
Solve m using NLP minimizing objvar;
]]></plaintext></body>