$title 'Test for known XPRESS bug' (MIP01,SEQ=215) $onText In this test we run a smallish MIP. XPRESS used to return a bogus nonbasic status for some of the continous free vars, resulting in a bad value for variable x43 in the model below Contributed by Steve Dirkse, February 2005 $offText * MIP written by GAMS Convert at 12/23/04 10:31:59 * * Equation counts * Total E G L N X C * 151 63 45 42 0 0 0 * * Variable counts * x b i s1s s2s sc si * Total cont binary integer sos1 sos2 scont sint * 153 132 21 0 0 0 0 0 * FX 5 5 0 0 0 0 0 0 * * Nonzero counts * Total const NL DLL * 490 490 0 0 * * Solve m using MIP minimizing x108; 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,b46,b47,b48,b49,b50,b51,b52,b53 ,b54,b55,b56,b57,b58,b59,b60,b61,b62,b63,b64,b65,b66,x67,x68,x69,x70 ,x71,x72,x73,x74,x75,x76,x77,x78,x79,x80,x81,x82,x83,x84,x85,x86,x87 ,x88,x89,x90,x91,x92,x93,x94,x95,x96,x97,x98,x99,x100,x101,x102,x103 ,x104,x105,x106,x107,x108,x109,x110,x111,x112,x113,x114,x115,x116 ,x117,x118,x119,x120,x121,x122,x123,x124,x125,x126,x127,x128,x129 ,x130,x131,x132,x133,x134,x135,x136,x137,x138,x139,x140,x141,x142 ,x143,x144,x145,x146,x147,x148,x149,x150,x151,x152,x153; Negative Variables x68,x69,x71,x72,x74,x75,x77,x78,x79; Positive Variables x10,x11,x22,x67,x70,x73,x76,x80,x81,x83,x84,x86,x87,x89 ,x90,x93,x96,x97,x98,x100,x101,x103,x104,x105,x106,x109,x110,x111 ,x112,x113,x114,x115,x116,x117,x118,x119,x120,x121,x122,x123,x124 ,x125,x126,x127,x128,x129,x130,x131,x132,x133,x134,x135,x136,x137 ,x138,x139,x140,x141,x142,x143,x144,x145,x146,x147,x148,x149,x150 ,x151,x152,x153; Binary Variables b46,b47,b48,b49,b50,b51,b52,b53,b54,b55,b56,b57,b58,b59,b60 ,b61,b62,b63,b64,b65,b66; 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,e47,e48,e49,e50,e51,e52,e53 ,e54,e55,e56,e57,e58,e59,e60,e61,e62,e63,e64,e65,e66,e67,e68,e69,e70 ,e71,e72,e73,e74,e75,e76,e77,e78,e79,e80,e81,e82,e83,e84,e85,e86,e87 ,e88,e89,e90,e91,e92,e93,e94,e95,e96,e97,e98,e99,e100,e101,e102,e103 ,e104,e105,e106,e107,e108,e109,e110,e111,e112,e113,e114,e115,e116 ,e117,e118,e119,e120,e121,e122,e123,e124,e125,e126,e127,e128,e129 ,e130,e131,e132,e133,e134,e135,e136,e137,e138,e139,e140,e141,e142 ,e143,e144,e145,e146,e147,e148,e149,e150,e151; e1.. - 2.28*x73 - 2.28*x76 - 1.68*x80 - 1.68*x82 - 2.28*x85 + x107 =E= 0; e2.. - x17 + x18 - x67 =E= 0; e3.. - x6 + x7 - x68 =E= 0; e4.. - x23 + x24 - x69 =E= 0; e5.. - x10 - x11 + x12 + x13 - x70 =E= 0; e6.. - x20 - x71 =E= 0; e7.. - x3 - x4 + x5 - x72 =E= 0; e8.. - x1 - x2 + x3 + x4 - x73 =E= 0; e9.. - x7 + x8 - x74 =E= 0; e10.. - x12 - x13 + x14 + x15 - x75 =E= 0; e11.. x6 - x76 =E= 0; e12.. - x19 + x20 - x77 =E= 0; e13.. - x16 + x17 - x78 =E= 0; e14.. - x24 - x79 =E= 0; e15.. - x9 - x18 + x19 - x80 =E= 0; e16.. - x22 + x23 - x81 =E= 0; e17.. x10 + x11 - x82 =E= 0; e18.. - x21 + x22 - x83 =E= 0; e19.. - x14 - x15 + x16 + x21 - x84 =E= 0; e20.. x1 + x2 - x85 =E= 0; e21.. - x5 - x8 + x9 - x86 =E= 0; e22.. 0.171441607948668*x1 + 29.6475712987489*x25 + 8.99076279639865*x93 - 8.99076279639865*x105 + x109 - x110 =E= 0.933492412620295; e23.. 0.171441601272821*x2 + 29.6475736076706*x26 + 8.99076279639865*x93 - 8.99076279639865*x105 + x111 - x112 =E= 0.933492412620286; e24.. 0.608059438332994*x3 + 83.285226048484*x27 + 5.99384186426577*x92 - 5.99384186426577*x93 - x113 + x114 =E= 5.54919541047039; e25.. 0.608059651973782*x4 + 83.2851675244689*x28 + 5.99384186426577*x92 - 5.99384186426577*x93 - x115 + x116 =E= 5.54919541047829; e26.. 3.31225570710095*x5 + 205.467704811367*x29 - 1.38319427636902*x92 + 1.38319427636902*x106 + x117 - x118 =E= 47.4783526503163; e27.. 0.912001255667401*x6 + 15.1621100827725*x30 + 0.0993769948466698*x88 - 0.0993769948466698*x96 - x119 + x120 =E= 3.55120103838917; e28.. 0.000919458994335796*x7 + 0.0196405039299483*x31 - 0.147352095807131*x88 + 0.147352095807131*x94 + x121 - x122 =E= -0.000128857196643815; e29.. 1.02401253338158*x8 + 29.118376447065*x32 - 0.224905830442463*x94 + 0.224905830442463*x106 + x123 - x124 =E= -5.52571969313141; e30.. 1.92072589956292*x9 + 79.8878462793372*x33 + 0.653873657919902*x100 - 0.653873657919902*x106 + x125 - x126 =E= 17.1674483352018; e31.. 11.212740307335*x12 + 382.080065171013*x34 - 1.79815255927973*x90 + 1.79815255927973*x95 + x127 - x128 =E= 219.173959455898; e32.. 1.31505154723441*x13 + 6.07687352689234*x35 - 0.0265962529480588*x90 + 0.0265962529480588*x95 + x129 - x130 =E= 3.24177505045168; e33.. 5.8448830543316*x14 + 191.736473189268*x36 - 1.43852204742378*x95 + 1.43852204742378*x104 + x131 - x132 =E= 80.9334760412807; e34.. 0.67247517674705*x15 + 3.16877725664789*x37 - 0.021277002358447*x95 + 0.021277002358447*x104 + x133 - x134 =E= 1.19707707533436; e35.. 3.39936806740346*x16 + 181.872194406949*x38 + 0.856263123466538*x98 - 0.856263123466538*x104 - x135 + x136 =E= 45.8438775562256; e36.. 2.0705203224792*x17 + 129.18595445796*x39 + 0.899076279639865*x87 - 0.899076279639865*x98 - x137 + x138 =E= 23.5335338448435; e37.. 1.01884973395169*x18 + 157.904354074667*x40 - 7.19261023711892*x87 + 7.19261023711892*x100 - x139 + x140 =E= 12.802865678784; e38.. 4.11626427812865*x19 + 504.631296788268*x41 + 3.59630511855946*x97 - 3.59630511855946*x100 - x141 + x142 =E= 92.4677608161024; e39.. 4.84719203271307*x20 + 243.859456273986*x42 + 1.43852204742378*x91 - 1.43852204742378*x97 - x143 + x144 =E= 75.6937508253697; e40.. 0.839534658280487*x21 + 4.58388107520205*x43 + 0.0509935163731392*x103 - 0.0509935163731392*x104 + x145 - x146 =E= 1.79744371812274; e41.. 3.5361193480341*x23 + 4.58388107520205*x44 + 0.0016882734118691*x89 - 0.0016882734118691*x101 - x147 + x148 =E= 7.57083158624364; e42.. 0.44721468857422*x24 + 3.68256103370239*x45 - 0.0275751323938619*x89 + 0.0275751323938619*x99 - x149 + x150 =E= 0.858204991301027; e43.. - 36.8878419417062*x10 + 7.19261023711892*x90 - 7.19261023711892*x102 + x151 =G= -340.178220779075; e44.. - 7.13615799159523*x11 + 0.106385011792235*x90 - 0.106385011792235*x102 + x152 =G= -12.7311877202521; e45.. - 4.28200003512473*x22 + 0.00636349209089121*x101 - 0.00636349209089121*x103 + x153 =G= -4.58388107520205; e46.. x25 - 2*b46 =E= -1; e47.. x26 - 2*b47 =E= -1; e48.. x27 - 2*b48 =E= -1; e49.. x28 - 2*b49 =E= -1; e50.. x29 - 2*b50 =E= -1; e51.. x30 - 2*b51 =E= -1; e52.. x31 - 2*b52 =E= -1; e53.. x32 - 2*b53 =E= -1; e54.. x33 - 2*b54 =E= -1; e55.. x34 - 2*b55 =E= -1; e56.. x35 - 2*b56 =E= -1; e57.. x36 - 2*b57 =E= -1; e58.. x37 - 2*b58 =E= -1; e59.. x38 - 2*b59 =E= -1; e60.. x39 - 2*b60 =E= -1; e61.. x40 - 2*b61 =E= -1; e62.. x41 - 2*b62 =E= -1; e63.. x42 - 2*b63 =E= -1; e64.. x43 - 2*b64 =E= -1; e65.. x44 - 2*b65 =E= -1; e66.. x45 - 2*b66 =E= -1; e67.. - 5929*b46 - x93 + x105 =G= -5929; e68.. - 5929*b47 - x93 + x105 =G= -5929; e69.. - 6400*b48 - x92 + x93 =G= -6400; e70.. - 6400*b49 - x92 + x93 =G= -6400; e71.. - 5500*b50 + x92 - x106 =G= -5500; e72.. - 6400*b51 - x88 + x96 =G= -6400; e73.. - 5500*b52 + x88 - x94 =G= -5500; e74.. - 5500*b53 + x94 - x106 =G= -5500; e75.. - 4382.44*b54 - x100 + x106 =G= -4382.44; e76.. - 4382.44*b55 + x90 - x95 =G= -4382.44; e77.. - 4382.44*b56 + x90 - x95 =G= -4382.44; e78.. - 3482.44*b57 + x95 - x104 =G= -3482.44; e79.. - 3482.44*b58 + x95 - x104 =G= -3482.44; e80.. - 4382.44*b59 - x98 + x104 =G= -4382.44; e81.. - 4382.44*b60 - x87 + x98 =G= -4382.44; e82.. - 4382.44*b61 + x87 - x100 =G= -4382.44; e83.. - 4382.44*b62 - x97 + x100 =G= -4382.44; e84.. - 4382.44*b63 - x91 + x97 =G= -4382.44; e85.. - 4382.44*b64 - x103 + x104 =G= -4382.44; e86.. - 4382.44*b65 - x89 + x101 =G= -4382.44; e87.. - 4382.44*b66 + x89 - x99 =G= -4382.44; e88.. - 5929*b46 - x93 + x105 =L= 0; e89.. - 5929*b47 - x93 + x105 =L= 0; e90.. - 5029*b48 - x92 + x93 =L= 0; e91.. - 5029*b49 - x92 + x93 =L= 0; e92.. - 6400*b50 + x92 - x106 =L= 0; e93.. - 5029*b51 - x88 + x96 =L= 0; e94.. - 5500*b52 + x88 - x94 =L= 0; e95.. - 6400*b53 + x94 - x106 =L= 0; e96.. - 6400*b54 - x100 + x106 =L= 0; e97.. - 3482.44*b55 + x90 - x95 =L= 0; e98.. - 3482.44*b56 + x90 - x95 =L= 0; e99.. - 4382.44*b57 + x95 - x104 =L= 0; e100.. - 4382.44*b58 + x95 - x104 =L= 0; e101.. - 4382.44*b59 - x98 + x104 =L= 0; e102.. - 4382.44*b60 - x87 + x98 =L= 0; e103.. - 4382.44*b61 + x87 - x100 =L= 0; e104.. - 4382.44*b62 - x97 + x100 =L= 0; e105.. - 1882.44*b63 - x91 + x97 =L= 0; e106.. - 4382.44*b64 - x103 + x104 =L= 0; e107.. - 3969*b65 - x89 + x101 =L= 0; e108.. - 3757.44*b66 + x89 - x99 =L= 0; e109.. x1 - 230.881425454383*b46 =G= -230.881425454383; e110.. x2 - 230.881425454383*b47 =G= -230.881425454383; e111.. x3 - 195.858591670881*b48 =G= -195.858591670881; e112.. x4 - 195.858591670881*b49 =G= -195.858591670881; e113.. x5 - 87.2213765084548*b50 =G= -87.2213765084548; e114.. x6 - 25.2192935471771*b51 =G= -25.2192935471771; e115.. x7 - 28.4681669051455*b52 =G= -28.4681669051455; e116.. x8 - 35.1707558553061*b53 =G= -35.1707558553061; e117.. x9 - 53.5309450076729*b54 =G= -53.5309450076729; e118.. x12 - 88.7710296317997*b55 =G= -88.7710296317997; e119.. x13 - 10.7961327691767*b56 =G= -10.7961327691767; e120.. x14 - 70.7782927092091*b57 =G= -70.7782927092091; e121.. x15 - 8.6078966125965*b58 =G= -8.6078966125965; e122.. x16 - 61.2578302162646*b59 =G= -61.2578302162646; e123.. x17 - 62.7705970255575*b60 =G= -62.7705970255575; e124.. x18 - 177.542059263599*b61 =G= -177.542059263599; e125.. x19 - 125.541194051115*b62 =G= -125.541194051115; e126.. x20 - 79.3992226757409*b63 =G= -79.3992226757409; e127.. x21 - 14.9491145521834*b64 =G= -14.9491145521834; e128.. x23 - 2.72006561154535*b65 =G= -2.72006561154535; e129.. x24 - 10.9930142912741*b66 =G= -10.9930142912741; e130.. x1 - 230.881425454383*b46 =L= 0; e131.. x2 - 230.881425454383*b47 =L= 0; e132.. x3 - 173.61748395652*b48 =L= 0; e133.. x4 - 173.61748395652*b49 =L= 0; e134.. x5 - 94.0874240733678*b50 =L= 0; e135.. x6 - 22.3554670513479*b51 =L= 0; e136.. x7 - 28.4681669051455*b52 =L= 0; e137.. x8 - 37.939390016601*b53 =L= 0; e138.. x9 - 64.6899637554959*b54 =L= 0; e139.. x12 - 79.1325369145847*b55 =L= 0; e140.. x13 - 9.62392098452798*b56 =L= 0; e141.. x14 - 79.3992226757409*b57 =L= 0; e142.. x15 - 9.65635470639685*b58 =L= 0; e143.. x16 - 61.2578302162646*b59 =L= 0; e144.. x17 - 62.7705970255575*b60 =L= 0; e145.. x18 - 177.542059263599*b61 =L= 0; e146.. x19 - 125.541194051115*b62 =L= 0; e147.. x20 - 52.0377886055166*b63 =L= 0; e148.. x21 - 14.9491145521834*b64 =L= 0; e149.. x23 - 2.58858207745253*b65 =L= 0; e150.. x24 - 10.1789933422708*b66 =L= 0; e151.. x107 - x108 + 2.56805524313174E-11*x109 + 2.55391279970258E-11*x111 + 2.01943285039219E-10*x113 + 1.97489763499426E-10*x115 + 4.72458685501564E-9*x117 + 4.02476419666052E-10*x119 + 2.8217280183436E-8*x121 + 3.60569518859414E-8*x123 + 2.75669248700102E-8*x125 + 3.93809302190442E-8*x127 + 1.96561974510839E-6*x129 + 5.46984377754984E-8*x131 + 2.18934836443633E-6*x133 + 2.61864978361483E-8*x135 + 3.32498218656729E-8*x137 + 5.23410981254842E-9*x139 + 8.43370081030883E-9*x141 + 3.25600083985981E-8*x143 + 2.61125935780483E-6*x145 + 1.35576920276201E-5*x147 + 4.07028955496529E-7*x149 + 7.7348160445682E-9*x151 + 6.60109445593242E-7*x152 + 1.92513685100548E-5*x153 =E= 0; * set non default bounds x67.up = 1.2; x68.up = -4.034; x69.up = -0.222; x70.fx = 0; x71.up = -15.616; x72.up = -3.918; x73.up = 8.4; x74.up = -5.256; x75.up = -6.385; x76.up = 4.8; x77.up = -6.848; x78.up = -2.12; x79.up = -1.919; x80.up = 0.96; x81.fx = 0; x82.lo = 20.344; x82.up = 22.012; x83.fx = 0; x84.fx = 0; x85.lo = 8.87; x85.up = 11.594; x86.fx = 0; x87.up = 4382.44; x88.lo = 900; x88.up = 6400; x89.up = 4382.44; x90.up = 4382.44; x91.lo = 2500; x91.up = 4382.44; x92.lo = 900; x92.up = 6400; x93.up = 5929; x94.lo = 900; x94.up = 6400; x95.lo = 900; x95.up = 4382.44; x96.up = 5929; x97.up = 4382.44; x98.up = 4382.44; x99.lo = 625; x99.up = 4382.44; x100.up = 4382.44; x101.up = 3969; x102.lo = 2500; x102.up = 4382.44; x103.up = 4382.44; x104.up = 4382.44; x105.up = 5929; x106.up = 6400; * set non default levels x1.l = 5.44495833765043; x2.l = 5.44495854967424; x3.l = 9.12607396685365; x4.l = 9.12607076043512; x5.l = 14.3341447185162; x6.l = 3.89385542653711; x7.l = -0.140144582235448; x8.l = -5.39614459100801; x9.l = 8.93800012750823; x10.l = 18.4439209708531; x11.l = 3.56807899579761; x12.l = 19.5468684236379; x13.l = 2.46513154352711; x14.l = 13.846893990685; x15.l = 1.78010596781424; x16.l = 13.4859999409369; x17.l = 11.3659999321643; x18.l = 12.5659999233912; x19.l = 22.4640000175451; x20.l = 15.6160000087726; x21.l = 2.14100001756236; x22.l = 2.14100001756236; x23.l = 2.14100001756236; x24.l = 1.91900000878124; x25.l = 0.0157431514914628; x26.l = 0.0157431502654026; x27.l = 0.0333144044493556; x28.l = 0.0333144278592461; x29.l = 0.115537263371644; x30.l = 0.11710774486541; x31.l = -0.00328039436012971; x32.l = -0.0948837189322137; x33.l = 0.107447184614126; x34.l = 0.286816795005778; x35.l = 0.266730501803737; x36.l = 0.211053939542816; x37.l = 0.188886276689687; x38.l = 0.126033222686167; x39.l = 0.0910839492713656; x40.l = 0.0405399387300305; x41.l = 0.0916191300506078; x42.l = 0.155199539894661; x43.l = 0.19606133848526; x44.l = 0.825810209955102; x45.l = 0.116522846932723; b46.l = 0.507871575745731; b47.l = 0.507871575132701; b48.l = 0.516657202224678; b49.l = 0.516657213929623; b50.l = 0.557768631685822; b51.l = 0.558553872432705; b52.l = 0.498359802819935; b53.l = 0.452558140533893; b54.l = 0.553723592307063; b55.l = 0.643408397502889; b56.l = 0.633365250901869; b57.l = 0.605526969771408; b58.l = 0.594443138344844; b59.l = 0.563016611343083; b60.l = 0.545541974635683; b61.l = 0.520269969365015; b62.l = 0.545809565025304; b63.l = 0.577599769947331; b64.l = 0.59803066924263; b65.l = 0.912905104977551; b66.l = 0.558261423466361; x67.l = 1.19999999122682; x68.l = -4.03400000877256; x69.l = -0.222000008781126; x71.l = -15.6160000087726; x72.l = -3.91800000877253; x73.l = 7.36222783968099; x74.l = -5.25600000877256; x75.l = -6.38500000877251; x76.l = 3.89385542647227; x77.l = -6.84800000877255; x78.l = -2.12000000877252; x79.l = -1.91900000878124; x80.l = 0.959999966645751; x82.l = 22.0119999666507; x85.l = 10.8899168872668; x87.l = 2909.59674752093; x88.l = 2909.54928693577; x89.l = 1595.41873634504; x90.l = 3038.5290822919; x91.l = 2869.54125150194; x92.l = 2938.99681832198; x93.l = 2939.4597263813; x94.l = 2909.54972417826; x95.l = 2977.58487634204; x96.l = 2927.41660644888; x97.l = 2895.85080684065; x98.l = 2922.68436554071; x99.l = 1579.8575238158; x100.l = 2908.70674616843; x101.l = 3837.6004609165; x102.l = 2642.44041110422; x103.l = 2931.82986493535; x104.l = 2949.45410273583; x105.l = 2939.51164035296; x106.l = 2921.83424273046; x107.l = 89.0858402375047; * set non default marginals Model m / all /; m.limrow=0; m.limcol=0; Solve m using MIP minimizing x108; * check that this equation is not violated * e64.. x43 - 2*b64 =E= -1; * this was a problem with XPRESS, * check PR XPRESS/60 or AB00710 abort$[round(x43.l - 2*b64.l + 1,5)] 'round(x43.l - 2*b64.l + 1,5) <> 0 ';