MINLPLib

A Library of Mixed-Integer and Continuous Nonlinear Programming Instances

Home // Instances // Documentation // Download // Statistics


Instance cesam2cent

Illustrates a cross entropy technique for estimating the cells of a consistent SAM assuming that the initial data are inconsistent and measured with error.
Formats ams gms mod nl osil py
Primal Bounds (infeas ≤ 1e-08)
0.50796042 p1 ( gdx sol )
(infeas: 4e-15)
Other points (infeas > 1e-08)  
Dual Bounds  
References Robinson, S, Cattaneo, A, and El-Said, M, Updating and Estimating a Social Accounting Matrix Using Cross Entropy Methods, Economic Systems Research, 13:1, 2001, 47-64.
Golan, A, Judge, G, and Miller, D, Maximum Entropy Econometrics, John Wiley and Sons, 1996.
Judge, G and Mittelhammer, R C, An Information Theoretic Approach to Econometrics, Cambridge University Press, New York, NY, 2012.
Source GAMS Model Library model cesam2
Application Social Accounting Matrix Balancing
Added to library 18 Aug 2014
Problem type NLP
#Variables 316
#Binary Variables 0
#Integer Variables 0
#Nonlinear Variables 207
#Nonlinear Binary Variables 0
#Nonlinear Integer Variables 0
Objective Sense min
Objective type nonlinear
Objective curvature nonconcave
#Nonzeros in Objective 157
#Nonlinear Nonzeros in Objective 157
#Constraints 165
#Linear Constraints 124
#Quadratic Constraints 28
#Polynomial Constraints 0
#Signomial Constraints 0
#General Nonlinear Constraints 13
Operands in Gen. Nonlin. Functions centropy exp
Constraints curvature indefinite
#Nonzeros in Jacobian 663
#Nonlinear Nonzeros in Jacobian 69
#Nonzeros in (Upper-Left) Hessian of Lagrangian 226
#Nonzeros in Diagonal of Hessian of Lagrangian 170
#Blocks in Hessian of Lagrangian 179
Minimal blocksize in Hessian of Lagrangian 1
Maximal blocksize in Hessian of Lagrangian 6
Average blocksize in Hessian of Lagrangian 1.156425
#Semicontinuities 0
#Nonlinear Semicontinuities 0
#SOS type 1 0
#SOS type 2 0
Minimal coefficient 2.4525e-04
Maximal coefficient 1.1121e+01
Infeasibility of initial point 0.3965
Sparsity Jacobian Sparsity of Objective Gradient and Jacobian
Sparsity Hessian of Lagrangian Sparsity of Hessian of Lagrangian

$offlisting
*  
*  Equation counts
*      Total        E        G        L        N        X        C        B
*        166      166        0        0        0        0        0        0
*  
*  Variable counts
*                   x        b        i      s1s      s2s       sc       si
*      Total     cont   binary  integer     sos1     sos2    scont     sint
*        317      317        0        0        0        0        0        0
*  FX     53
*  
*  Nonzero counts
*      Total    const       NL      DLL
*        821      595      226        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,x66,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,x154,x155
          ,x156,x157,x158,x159,x160,x161,x162,x163,x164,x165,x166,x167,x168
          ,x169,x170,x171,x172,x173,x174,x175,x176,x177,x178,x179,x180,x181
          ,x182,x183,x184,x185,x186,x187,x188,x189,x190,x191,x192,x193,x194
          ,x195,x196,x197,x198,x199,x200,x201,x202,x203,x204,x205,x206,x207
          ,x208,x209,x210,x211,x212,x213,x214,x215,x216,x217,x218,x219,x220
          ,x221,x222,x223,x224,x225,x226,x227,x228,x229,x230,x231,x232,x233
          ,x234,x235,x236,x237,x238,x239,x240,x241,x242,x243,x244,x245,x246
          ,x247,x248,x249,x250,x251,x252,x253,x254,x255,x256,x257,x258,x259
          ,x260,x261,x262,x263,x264,x265,x266,x267,x268,x269,x270,x271,x272
          ,x273,x274,x275,x276,x277,x278,x279,x280,x281,x282,x283,x284,x285
          ,x286,x287,x288,x289,x290,x291,x292,x293,x294,x295,x296,x297,x298
          ,x299,x300,x301,x302,x303,x304,x305,x306,x307,x308,x309,x310,x311
          ,x312,x313,x314,x315,x316,objvar;

Positive Variables  x160,x161,x162,x163,x164,x165,x166,x167,x168,x169,x170
          ,x171,x172,x173,x174,x175,x176,x177,x178,x179,x180,x181,x182,x183
          ,x184,x185,x186,x187,x188,x189,x190,x191,x192,x193,x194,x195,x196
          ,x197,x198,x199,x200,x201,x202,x203,x204,x205,x206,x207,x208,x209
          ,x210,x211,x212,x213,x214,x215,x216,x217,x218,x219,x220,x221,x222
          ,x223,x224,x225,x226,x227,x228,x229,x230,x231,x232,x233,x234,x235
          ,x236,x237,x238,x239,x240,x241,x242,x243,x244,x245,x246,x247,x248
          ,x249,x250,x251,x252,x253,x254,x255,x256,x257,x258,x259,x260,x261
          ,x262,x263,x264,x265,x266,x267,x268,x269,x270,x271,x272,x273,x274
          ,x275,x276,x277,x278,x279,x280,x281,x282,x283,x284,x285,x286,x287
          ,x288,x289,x290,x291,x292,x293,x294,x295,x296,x297,x298,x299,x300
          ,x301,x302,x303,x304,x305,x306,x307,x308,x309,x310,x311,x312,x313
          ,x314,x315,x316;

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,e152,e153,e154,e155
          ,e156,e157,e158,e159,e160,e161,e162,e163,e164,e165,e166;


e1..    x112 - x121 =E= 18.4364105;

e2..    x113 - x122 =E= 21.1551365;

e3..    x114 - x123 =E= 9.78976;

e4..    x115 - x124 =E= 3.673953;

e5..    x116 - x125 =E= 9.6863185;

e6..    x117 - x126 =E= 1.3701;

e7..    x118 - x127 =E= 1.9123;

e8..    x119 - x128 =E= 2.398969;

e9..    x120 - x129 =E= 5.5690645;

e10..    x29 + x30 + x31 + x32 + x33 + x34 + x35 + x36 + x37 - x112 =E= 0;

e11..    x38 + x39 + x40 + x41 + x42 + x43 + x44 + x45 + x46 - x113 =E= 0;

e12..    x47 + x48 + x49 + x50 + x51 + x52 + x53 + x54 + x55 - x114 =E= 0;

e13..    x56 + x57 + x58 + x59 + x60 + x61 + x62 + x63 + x64 - x115 =E= 0;

e14..    x65 + x66 + x67 + x68 + x69 + x70 + x71 + x72 + x73 - x116 =E= 0;

e15..    x74 + x75 + x76 + x77 + x78 + x79 + x80 + x81 + x82 - x117 =E= 0;

e16..    x83 + x84 + x85 + x86 + x87 + x88 + x89 + x90 + x91 - x118 =E= 0;

e17..    x92 + x93 + x94 + x95 + x96 + x97 + x98 + x99 + x100 - x119 =E= 0;

e18..    x101 + x102 + x103 + x104 + x105 + x106 + x107 + x108 + x109 - x120
       =E= 0;

e19..    x29 + x38 + x47 + x56 + x65 + x74 + x83 + x92 + x101 - x112 =E= 0;

e20..    x30 + x39 + x48 + x57 + x66 + x75 + x84 + x93 + x102 - x113 =E= 0;

e21..    x31 + x40 + x49 + x58 + x67 + x76 + x85 + x94 + x103 - x114 =E= 0;

e22..    x32 + x41 + x50 + x59 + x68 + x77 + x86 + x95 + x104 - x115 =E= 0;

e23..    x33 + x42 + x51 + x60 + x69 + x78 + x87 + x96 + x105 - x116 =E= 0;

e24..    x34 + x43 + x52 + x61 + x70 + x79 + x88 + x97 + x106 - x117 =E= 0;

e25..    x35 + x44 + x53 + x62 + x71 + x80 + x89 + x98 + x107 - x118 =E= 0;

e26..    x36 + x45 + x54 + x63 + x72 + x81 + x90 + x99 + x108 - x119 =E= 0;

e27..    x37 + x46 + x55 + x64 + x73 + x82 + x91 + x100 + x109 - x120 =E= 0;

e28.. -x1*x113 + x30 =E= 0;

e29.. -x2*x116 + x33 =E= 0;

e30.. -x3*x117 + x34 =E= 0;

e31.. -x4*x120 + x37 =E= 0;

e32.. -x5*x112 + x38 =E= 0;

e33.. -x6*x116 + x42 =E= 0;

e34.. -x7*x117 + x43 =E= 0;

e35.. -x8*x118 + x44 =E= 0;

e36.. -x9*x119 + x45 =E= 0;

e37.. -x10*x112 + x47 =E= 0;

e38.. -x11*x114 + x58 =E= 0;

e39.. -x12*x117 + x61 =E= 0;

e40.. -x13*x114 + x67 =E= 0;

e41.. -x14*x115 + x68 =E= 0;

e42.. -x15*x117 + x70 =E= 0;

e43.. -x16*x120 + x73 =E= 0;

e44.. -x17*x112 + x74 =E= 0;

e45.. -x18*x113 + x75 =E= 0;

e46.. -x19*x114 + x76 =E= 0;

e47.. -x20*x115 + x77 =E= 0;

e48.. -x21*x116 + x78 =E= 0;

e49.. -x22*x120 + x91 =E= 0;

e50.. -x23*x115 + x95 =E= 0;

e51.. -x24*x116 + x96 =E= 0;

e52.. -x25*x117 + x97 =E= 0;

e53.. -x26*x118 + x98 =E= 0;

e54.. -x27*x120 + x100 =E= 0;

e55.. -x28*x113 + x102 =E= 0;

e56..    x30 - x132 =E= 14.827424;

e57..    x34 - x134 =E= -0.000327;

e58..    x37 - x135 =E= 1.488157;

e59..    x43 - x138 =E= 1.5645;

e60..    x44 - x139 =E= 2.5185;

e61..    x45 - x140 =E= 2.597798;

e62..    x61 - x143 =E= 0.033;

e63..    x70 - x146 =E= 0.0296;

e64..    x73 - x147 =E= 0.2;

e65..    x75 - x149 =E= 0.3574;

e66..    x91 - x153 =E= 1.7123;

e67..    x97 - x156 =E= -0.356673;

e68..    x98 - x157 =E= -0.4062;

e69..    x100 - x158 =E= 2.163857;

e70..    x102 - x159 =E= 5.573815;

e71.. -0.213455359357076*exp(x133) + x2 =E= 0;

e72.. -0.428981457932639*exp(x136) + x5 =E= 0;

e73.. -0.706421402256235*exp(x137) + x6 =E= 0;

e74.. -0.531271066405917*exp(x141) + x10 =E= 0;

e75.. -0.37852116602787*exp(x142) + x11 =E= 0;

e76.. -0.613866884603052*exp(x144) + x13 =E= 0;

e77.. -0.912812569152467*exp(x145) + x14 =E= 0;

e78.. -0.0397474756614438*exp(x148) + x17 =E= 0;

e79.. -0.00761194936907785*exp(x150) + x19 =E= 0;

e80.. -0.0456959504315114*exp(x151) + x20 =E= 0;

e81.. -0.0141724551070975*exp(x152) + x21 =E= 0;

e82.. -0.0414914804160212*exp(x154) + x23 =E= 0;

e83.. -0.0659507832795914*exp(x155) + x24 =E= 0;

e84..  - x47 + x110 =E= 0;

e85..    x34 - x47 - x74 - x75 + x111 =E= 0;

e86..    x110 - x130 =E= 9.805414;

e87..    x111 - x131 =E= 10.896741;

e88..    x121 + 2.765461575*x160 + 1.84364105*x161 + 0.921820525*x162
       - 0.921820525*x164 - 1.84364105*x165 - 2.765461575*x166 =E= 0;

e89..    x122 + 3.173270475*x167 + 2.11551365*x168 + 1.057756825*x169
       - 1.057756825*x171 - 2.11551365*x172 - 3.173270475*x173 =E= 0;

e90..    x123 + 1.468464*x174 + 0.978976*x175 + 0.489488*x176 - 0.489488*x178
       - 0.978976*x179 - 1.468464*x180 =E= 0;

e91..    x124 + 0.55109295*x181 + 0.3673953*x182 + 0.18369765*x183
       - 0.18369765*x185 - 0.3673953*x186 - 0.55109295*x187 =E= 0;

e92..    x125 + 1.452947775*x188 + 0.96863185*x189 + 0.484315925*x190
       - 0.484315925*x192 - 0.96863185*x193 - 1.452947775*x194 =E= 0;

e93..    x126 + 0.205515*x195 + 0.13701*x196 + 0.068505*x197 - 0.068505*x199
       - 0.13701*x200 - 0.205515*x201 =E= 0;

e94..    x127 + 0.286845*x202 + 0.19123*x203 + 0.095615*x204 - 0.095615*x206
       - 0.19123*x207 - 0.286845*x208 =E= 0;

e95..    x128 + 0.35984535*x209 + 0.2398969*x210 + 0.11994845*x211
       - 0.11994845*x213 - 0.2398969*x214 - 0.35984535*x215 =E= 0;

e96..    x129 + 0.835359675*x216 + 0.55690645*x217 + 0.278453225*x218
       - 0.278453225*x220 - 0.55690645*x221 - 0.835359675*x222 =E= 0;

e97..    x130 + 1.4708121*x223 + 0.73540605*x224 - 0.73540605*x226
       - 1.4708121*x227 =E= 0;

e98..    x131 + 1.63451115*x228 + 0.817255575*x229 - 0.817255575*x231
       - 1.63451115*x232 =E= 0;

e99..    x132 + 11.120568*x233 - 11.120568*x235 =E= 0;

e100..    x133 + 0.75*x236 - 0.75*x238 =E= 0;

e101..    x134 + 0.00024525*x239 - 0.00024525*x241 =E= 0;

e102..    x135 + 1.11611775*x242 - 1.11611775*x244 =E= 0;

e103..    x136 + 0.75*x245 - 0.75*x247 =E= 0;

e104..    x137 + 0.75*x248 - 0.75*x250 =E= 0;

e105..    x138 + 1.173375*x251 - 1.173375*x253 =E= 0;

e106..    x139 + 1.888875*x254 - 1.888875*x256 =E= 0;

e107..    x140 + 1.9483485*x257 - 1.9483485*x259 =E= 0;

e108..    x141 + 0.75*x260 - 0.75*x262 =E= 0;

e109..    x142 + 0.75*x263 - 0.75*x265 =E= 0;

e110..    x143 + 0.02475*x266 - 0.02475*x268 =E= 0;

e111..    x144 + 0.75*x269 - 0.75*x271 =E= 0;

e112..    x145 + 0.75*x272 - 0.75*x274 =E= 0;

e113..    x146 + 0.0222*x275 - 0.0222*x277 =E= 0;

e114..    x147 + 0.15*x278 - 0.15*x280 =E= 0;

e115..    x148 + 0.75*x281 - 0.75*x283 =E= 0;

e116..    x149 + 0.26805*x284 - 0.26805*x286 =E= 0;

e117..    x150 + 0.75*x287 - 0.75*x289 =E= 0;

e118..    x151 + 0.75*x290 - 0.75*x292 =E= 0;

e119..    x152 + 0.75*x293 - 0.75*x295 =E= 0;

e120..    x153 + 1.284225*x296 - 1.284225*x298 =E= 0;

e121..    x154 + 0.75*x299 - 0.75*x301 =E= 0;

e122..    x155 + 0.75*x302 - 0.75*x304 =E= 0;

e123..    x156 + 0.26750475*x305 - 0.26750475*x307 =E= 0;

e124..    x157 + 0.30465*x308 - 0.30465*x310 =E= 0;

e125..    x158 + 1.62289275*x311 - 1.62289275*x313 =E= 0;

e126..    x159 + 4.18036125*x314 - 4.18036125*x316 =E= 0;

e127..    x160 + x161 + x162 + x163 + x164 + x165 + x166 =E= 1;

e128..    x167 + x168 + x169 + x170 + x171 + x172 + x173 =E= 1;

e129..    x174 + x175 + x176 + x177 + x178 + x179 + x180 =E= 1;

e130..    x181 + x182 + x183 + x184 + x185 + x186 + x187 =E= 1;

e131..    x188 + x189 + x190 + x191 + x192 + x193 + x194 =E= 1;

e132..    x195 + x196 + x197 + x198 + x199 + x200 + x201 =E= 1;

e133..    x202 + x203 + x204 + x205 + x206 + x207 + x208 =E= 1;

e134..    x209 + x210 + x211 + x212 + x213 + x214 + x215 =E= 1;

e135..    x216 + x217 + x218 + x219 + x220 + x221 + x222 =E= 1;

e136..    x223 + x224 + x225 + x226 + x227 =E= 1;

e137..    x228 + x229 + x230 + x231 + x232 =E= 1;

e138..    x233 + x234 + x235 =E= 1;

e139..    x236 + x237 + x238 =E= 1;

e140..    x239 + x240 + x241 =E= 1;

e141..    x242 + x243 + x244 =E= 1;

e142..    x245 + x246 + x247 =E= 1;

e143..    x248 + x249 + x250 =E= 1;

e144..    x251 + x252 + x253 =E= 1;

e145..    x254 + x255 + x256 =E= 1;

e146..    x257 + x258 + x259 =E= 1;

e147..    x260 + x261 + x262 =E= 1;

e148..    x263 + x264 + x265 =E= 1;

e149..    x266 + x267 + x268 =E= 1;

e150..    x269 + x270 + x271 =E= 1;

e151..    x272 + x273 + x274 =E= 1;

e152..    x275 + x276 + x277 =E= 1;

e153..    x278 + x279 + x280 =E= 1;

e154..    x281 + x282 + x283 =E= 1;

e155..    x284 + x285 + x286 =E= 1;

e156..    x287 + x288 + x289 =E= 1;

e157..    x290 + x291 + x292 =E= 1;

e158..    x293 + x294 + x295 =E= 1;

e159..    x296 + x297 + x298 =E= 1;

e160..    x299 + x300 + x301 =E= 1;

e161..    x302 + x303 + x304 =E= 1;

e162..    x305 + x306 + x307 =E= 1;

e163..    x308 + x309 + x310 =E= 1;

e164..    x311 + x312 + x313 =E= 1;

e165..    x314 + x315 + x316 =E= 1;

e166.. -(Centropy(x233,0.0555555555555556) + Centropy(x234,0.888888888888889)
        + Centropy(x235,0.0555555555555556) + Centropy(x236,0.0555555555555556)
        + Centropy(x237,0.888888888888889) + Centropy(x238,0.0555555555555556)
        + Centropy(x239,0.0555555555555556) + Centropy(x240,0.888888888888889)
        + Centropy(x241,0.0555555555555556) + Centropy(x242,0.0555555555555556)
        + Centropy(x243,0.888888888888889) + Centropy(x244,0.0555555555555556)
        + Centropy(x245,0.0555555555555556) + Centropy(x246,0.888888888888889)
        + Centropy(x247,0.0555555555555556) + Centropy(x248,0.0555555555555556)
        + Centropy(x249,0.888888888888889) + Centropy(x250,0.0555555555555556)
        + Centropy(x251,0.0555555555555556) + Centropy(x252,0.888888888888889)
        + Centropy(x253,0.0555555555555556) + Centropy(x254,0.0555555555555556)
        + Centropy(x255,0.888888888888889) + Centropy(x256,0.0555555555555556)
        + Centropy(x257,0.0555555555555556) + Centropy(x258,0.888888888888889)
        + Centropy(x259,0.0555555555555556) + Centropy(x260,0.0555555555555556)
        + Centropy(x261,0.888888888888889) + Centropy(x262,0.0555555555555556)
        + Centropy(x263,0.0555555555555556) + Centropy(x264,0.888888888888889)
        + Centropy(x265,0.0555555555555556) + Centropy(x266,0.0555555555555556)
        + Centropy(x267,0.888888888888889) + Centropy(x268,0.0555555555555556)
        + Centropy(x269,0.0555555555555556) + Centropy(x270,0.888888888888889)
        + Centropy(x271,0.0555555555555556) + Centropy(x272,0.0555555555555556)
        + Centropy(x273,0.888888888888889) + Centropy(x274,0.0555555555555556)
        + Centropy(x275,0.0555555555555556) + Centropy(x276,0.888888888888889)
        + Centropy(x277,0.0555555555555556) + Centropy(x278,0.0555555555555556)
        + Centropy(x279,0.888888888888889) + Centropy(x280,0.0555555555555556)
        + Centropy(x281,0.0555555555555556) + Centropy(x282,0.888888888888889)
        + Centropy(x283,0.0555555555555556) + Centropy(x284,0.0555555555555556)
        + Centropy(x285,0.888888888888889) + Centropy(x286,0.0555555555555556)
        + Centropy(x287,0.0555555555555556) + Centropy(x288,0.888888888888889)
        + Centropy(x289,0.0555555555555556) + Centropy(x290,0.0555555555555556)
        + Centropy(x291,0.888888888888889) + Centropy(x292,0.0555555555555556)
        + Centropy(x293,0.0555555555555556) + Centropy(x294,0.888888888888889)
        + Centropy(x295,0.0555555555555556) + Centropy(x296,0.0555555555555556)
        + Centropy(x297,0.888888888888889) + Centropy(x298,0.0555555555555556)
        + Centropy(x299,0.0555555555555556) + Centropy(x300,0.888888888888889)
        + Centropy(x301,0.0555555555555556) + Centropy(x302,0.0555555555555556)
        + Centropy(x303,0.888888888888889) + Centropy(x304,0.0555555555555556)
        + Centropy(x305,0.0555555555555556) + Centropy(x306,0.888888888888889)
        + Centropy(x307,0.0555555555555556) + Centropy(x308,0.0555555555555556)
        + Centropy(x309,0.888888888888889) + Centropy(x310,0.0555555555555556)
        + Centropy(x311,0.0555555555555556) + Centropy(x312,0.888888888888889)
        + Centropy(x313,0.0555555555555556) + Centropy(x314,0.0555555555555556)
        + Centropy(x315,0.888888888888889) + Centropy(x316,0.0555555555555556)
        + Centropy(x160,0.142857142857143) + Centropy(x161,0.142857142857143)
        + Centropy(x162,0.142857142857143) + Centropy(x163,0.142857142857143)
        + Centropy(x164,0.142857142857143) + Centropy(x165,0.142857142857143)
        + Centropy(x166,0.142857142857143) + Centropy(x167,0.142857142857143)
        + Centropy(x168,0.142857142857143) + Centropy(x169,0.142857142857143)
        + Centropy(x170,0.142857142857143) + Centropy(x171,0.142857142857143)
        + Centropy(x172,0.142857142857143) + Centropy(x173,0.142857142857143)
        + Centropy(x174,0.142857142857143) + Centropy(x175,0.142857142857143)
        + Centropy(x176,0.142857142857143) + Centropy(x177,0.142857142857143)
        + Centropy(x178,0.142857142857143) + Centropy(x179,0.142857142857143)
        + Centropy(x180,0.142857142857143) + Centropy(x181,0.142857142857143)
        + Centropy(x182,0.142857142857143) + Centropy(x183,0.142857142857143)
        + Centropy(x184,0.142857142857143) + Centropy(x185,0.142857142857143)
        + Centropy(x186,0.142857142857143) + Centropy(x187,0.142857142857143)
        + Centropy(x188,0.142857142857143) + Centropy(x189,0.142857142857143)
        + Centropy(x190,0.142857142857143) + Centropy(x191,0.142857142857143)
        + Centropy(x192,0.142857142857143) + Centropy(x193,0.142857142857143)
        + Centropy(x194,0.142857142857143) + Centropy(x195,0.142857142857143)
        + Centropy(x196,0.142857142857143) + Centropy(x197,0.142857142857143)
        + Centropy(x198,0.142857142857143) + Centropy(x199,0.142857142857143)
        + Centropy(x200,0.142857142857143) + Centropy(x201,0.142857142857143)
        + Centropy(x202,0.142857142857143) + Centropy(x203,0.142857142857143)
        + Centropy(x204,0.142857142857143) + Centropy(x205,0.142857142857143)
        + Centropy(x206,0.142857142857143) + Centropy(x207,0.142857142857143)
        + Centropy(x208,0.142857142857143) + Centropy(x209,0.142857142857143)
        + Centropy(x210,0.142857142857143) + Centropy(x211,0.142857142857143)
        + Centropy(x212,0.142857142857143) + Centropy(x213,0.142857142857143)
        + Centropy(x214,0.142857142857143) + Centropy(x215,0.142857142857143)
        + Centropy(x216,0.142857142857143) + Centropy(x217,0.142857142857143)
        + Centropy(x218,0.142857142857143) + Centropy(x219,0.142857142857143)
        + Centropy(x220,0.142857142857143) + Centropy(x221,0.142857142857143)
        + Centropy(x222,0.142857142857143) + Centropy(x223,0.00617283950617284)
        + Centropy(x224,0.197530864197531) + Centropy(x225,0.592592592592593)
        + Centropy(x226,0.197530864197531) + Centropy(x227,0.00617283950617284)
        + Centropy(x228,0.00617283950617284) + Centropy(x229,0.197530864197531)
        + Centropy(x230,0.592592592592593) + Centropy(x231,0.197530864197531)
        + Centropy(x232,0.00617283950617284)) + objvar =E= 0;

* set non-default bounds
x29.fx = 0;
x31.fx = 0;
x32.fx = 0;
x35.fx = 0;
x36.fx = 0;
x39.fx = 0;
x40.fx = 0;
x41.fx = 0;
x46.fx = 0;
x48.fx = 0;
x49.fx = 0;
x50.fx = 0;
x51.fx = 0;
x52.fx = 0;
x53.fx = 0;
x54.fx = 0;
x55.fx = 0;
x56.fx = 0;
x57.fx = 0;
x59.fx = 0;
x60.fx = 0;
x62.fx = 0;
x63.fx = 0;
x64.fx = 0;
x65.fx = 0;
x66.fx = 0;
x69.fx = 0;
x71.fx = 0;
x72.fx = 0;
x79.fx = 0;
x80.fx = 0;
x81.fx = 0;
x82.fx = 0;
x83.fx = 0;
x84.fx = 0;
x85.fx = 0;
x86.fx = 0;
x87.fx = 0;
x88.fx = 0;
x89.fx = 0;
x90.fx = 0;
x92.fx = 0;
x93.fx = 0;
x94.fx = 0;
x99.fx = 0;
x101.fx = 0;
x103.fx = 0;
x104.fx = 0;
x105.fx = 0;
x106.fx = 0;
x107.fx = 0;
x108.fx = 0;
x109.fx = 0;
x160.up = 1;
x161.up = 1;
x162.up = 1;
x163.up = 1;
x164.up = 1;
x165.up = 1;
x166.up = 1;
x167.up = 1;
x168.up = 1;
x169.up = 1;
x170.up = 1;
x171.up = 1;
x172.up = 1;
x173.up = 1;
x174.up = 1;
x175.up = 1;
x176.up = 1;
x177.up = 1;
x178.up = 1;
x179.up = 1;
x180.up = 1;
x181.up = 1;
x182.up = 1;
x183.up = 1;
x184.up = 1;
x185.up = 1;
x186.up = 1;
x187.up = 1;
x188.up = 1;
x189.up = 1;
x190.up = 1;
x191.up = 1;
x192.up = 1;
x193.up = 1;
x194.up = 1;
x195.up = 1;
x196.up = 1;
x197.up = 1;
x198.up = 1;
x199.up = 1;
x200.up = 1;
x201.up = 1;
x202.up = 1;
x203.up = 1;
x204.up = 1;
x205.up = 1;
x206.up = 1;
x207.up = 1;
x208.up = 1;
x209.up = 1;
x210.up = 1;
x211.up = 1;
x212.up = 1;
x213.up = 1;
x214.up = 1;
x215.up = 1;
x216.up = 1;
x217.up = 1;
x218.up = 1;
x219.up = 1;
x220.up = 1;
x221.up = 1;
x222.up = 1;
x223.up = 1;
x224.up = 1;
x225.up = 1;
x226.up = 1;
x227.up = 1;
x228.up = 1;
x229.up = 1;
x230.up = 1;
x231.up = 1;
x232.up = 1;
x233.up = 1;
x234.up = 1;
x235.up = 1;
x236.up = 1;
x237.up = 1;
x238.up = 1;
x239.up = 1;
x240.up = 1;
x241.up = 1;
x242.up = 1;
x243.up = 1;
x244.up = 1;
x245.up = 1;
x246.up = 1;
x247.up = 1;
x248.up = 1;
x249.up = 1;
x250.up = 1;
x251.up = 1;
x252.up = 1;
x253.up = 1;
x254.up = 1;
x255.up = 1;
x256.up = 1;
x257.up = 1;
x258.up = 1;
x259.up = 1;
x260.up = 1;
x261.up = 1;
x262.up = 1;
x263.up = 1;
x264.up = 1;
x265.up = 1;
x266.up = 1;
x267.up = 1;
x268.up = 1;
x269.up = 1;
x270.up = 1;
x271.up = 1;
x272.up = 1;
x273.up = 1;
x274.up = 1;
x275.up = 1;
x276.up = 1;
x277.up = 1;
x278.up = 1;
x279.up = 1;
x280.up = 1;
x281.up = 1;
x282.up = 1;
x283.up = 1;
x284.up = 1;
x285.up = 1;
x286.up = 1;
x287.up = 1;
x288.up = 1;
x289.up = 1;
x290.up = 1;
x291.up = 1;
x292.up = 1;
x293.up = 1;
x294.up = 1;
x295.up = 1;
x296.up = 1;
x297.up = 1;
x298.up = 1;
x299.up = 1;
x300.up = 1;
x301.up = 1;
x302.up = 1;
x303.up = 1;
x304.up = 1;
x305.up = 1;
x306.up = 1;
x307.up = 1;
x308.up = 1;
x309.up = 1;
x310.up = 1;
x311.up = 1;
x312.up = 1;
x313.up = 1;
x314.up = 1;
x315.up = 1;
x316.up = 1;

* set non-default levels
x1.l = 0.714277270296959;
x2.l = 0.213455359357076;
x3.l = -0.000257460042516337;
x4.l = 0.267446625046681;
x5.l = 0.428981457932639;
x6.l = 0.706421402256235;
x7.l = 1.23179277222266;
x8.l = 1.1923022297969;
x9.l = 1;
x10.l = 0.531271066405917;
x11.l = 0.37852116602787;
x12.l = 0.0259822061255019;
x13.l = 0.613866884603052;
x14.l = 0.912812569152467;
x15.l = 0.0233052515549957;
x16.l = 0.0359433346141142;
x17.l = 0.0397474756614438;
x18.l = 0.0172169283352343;
x19.l = 0.00761194936907785;
x20.l = 0.0456959504315114;
x21.l = 0.0141724551070975;
x22.l = 0.307728859298738;
x23.l = 0.0414914804160212;
x24.l = 0.0659507832795914;
x25.l = -0.280822769860641;
x26.l = -0.192302229796904;
x27.l = 0.388881181040466;
x28.l = 0.268505801367806;
x30.l = 14.827424;
x33.l = 2.101049;
x34.l = -0.000327;
x37.l = 1.488157;
x38.l = 7.917504;
x42.l = 6.953332;
x43.l = 1.5645;
x44.l = 2.5185;
x45.l = 2.597798;
x47.l = 9.805414;
x58.l = 3.699706;
x61.l = 0.033;
x67.l = 6;
x68.l = 3.3;
x70.l = 0.0296;
x73.l = 0.2;
x74.l = 0.7336;
x75.l = 0.3574;
x76.l = 0.0744;
x77.l = 0.1652;
x78.l = 0.1395;
x91.l = 1.7123;
x95.l = 0.15;
x96.l = 0.649156;
x97.l = -0.356673;
x98.l = -0.4062;
x100.l = 2.163857;
x102.l = 5.573815;
x110.l = 9.805414;
x111.l = 10.896741;
x112.l = 18.4364105;
x113.l = 21.1551365;
x114.l = 9.78976;
x115.l = 3.673953;
x116.l = 9.6863185;
x117.l = 1.3701;
x118.l = 1.9123;
x119.l = 2.398969;
x120.l = 5.5690645;
x160.l = 0.142857142857143;
x161.l = 0.142857142857143;
x162.l = 0.142857142857143;
x163.l = 0.142857142857143;
x164.l = 0.142857142857143;
x165.l = 0.142857142857143;
x166.l = 0.142857142857143;
x167.l = 0.142857142857143;
x168.l = 0.142857142857143;
x169.l = 0.142857142857143;
x170.l = 0.142857142857143;
x171.l = 0.142857142857143;
x172.l = 0.142857142857143;
x173.l = 0.142857142857143;
x174.l = 0.142857142857143;
x175.l = 0.142857142857143;
x176.l = 0.142857142857143;
x177.l = 0.142857142857143;
x178.l = 0.142857142857143;
x179.l = 0.142857142857143;
x180.l = 0.142857142857143;
x181.l = 0.142857142857143;
x182.l = 0.142857142857143;
x183.l = 0.142857142857143;
x184.l = 0.142857142857143;
x185.l = 0.142857142857143;
x186.l = 0.142857142857143;
x187.l = 0.142857142857143;
x188.l = 0.142857142857143;
x189.l = 0.142857142857143;
x190.l = 0.142857142857143;
x191.l = 0.142857142857143;
x192.l = 0.142857142857143;
x193.l = 0.142857142857143;
x194.l = 0.142857142857143;
x195.l = 0.142857142857143;
x196.l = 0.142857142857143;
x197.l = 0.142857142857143;
x198.l = 0.142857142857143;
x199.l = 0.142857142857143;
x200.l = 0.142857142857143;
x201.l = 0.142857142857143;
x202.l = 0.142857142857143;
x203.l = 0.142857142857143;
x204.l = 0.142857142857143;
x205.l = 0.142857142857143;
x206.l = 0.142857142857143;
x207.l = 0.142857142857143;
x208.l = 0.142857142857143;
x209.l = 0.142857142857143;
x210.l = 0.142857142857143;
x211.l = 0.142857142857143;
x212.l = 0.142857142857143;
x213.l = 0.142857142857143;
x214.l = 0.142857142857143;
x215.l = 0.142857142857143;
x216.l = 0.142857142857143;
x217.l = 0.142857142857143;
x218.l = 0.142857142857143;
x219.l = 0.142857142857143;
x220.l = 0.142857142857143;
x221.l = 0.142857142857143;
x222.l = 0.142857142857143;
x223.l = 0.00617283950617284;
x224.l = 0.197530864197531;
x225.l = 0.592592592592593;
x226.l = 0.197530864197531;
x227.l = 0.00617283950617284;
x228.l = 0.00617283950617284;
x229.l = 0.197530864197531;
x230.l = 0.592592592592593;
x231.l = 0.197530864197531;
x232.l = 0.00617283950617284;
x233.l = 0.0555555555555556;
x234.l = 0.888888888888889;
x235.l = 0.0555555555555556;
x236.l = 0.0555555555555556;
x237.l = 0.888888888888889;
x238.l = 0.0555555555555556;
x239.l = 0.0555555555555556;
x240.l = 0.888888888888889;
x241.l = 0.0555555555555556;
x242.l = 0.0555555555555556;
x243.l = 0.888888888888889;
x244.l = 0.0555555555555556;
x245.l = 0.0555555555555556;
x246.l = 0.888888888888889;
x247.l = 0.0555555555555556;
x248.l = 0.0555555555555556;
x249.l = 0.888888888888889;
x250.l = 0.0555555555555556;
x251.l = 0.0555555555555556;
x252.l = 0.888888888888889;
x253.l = 0.0555555555555556;
x254.l = 0.0555555555555556;
x255.l = 0.888888888888889;
x256.l = 0.0555555555555556;
x257.l = 0.0555555555555556;
x258.l = 0.888888888888889;
x259.l = 0.0555555555555556;
x260.l = 0.0555555555555556;
x261.l = 0.888888888888889;
x262.l = 0.0555555555555556;
x263.l = 0.0555555555555556;
x264.l = 0.888888888888889;
x265.l = 0.0555555555555556;
x266.l = 0.0555555555555556;
x267.l = 0.888888888888889;
x268.l = 0.0555555555555556;
x269.l = 0.0555555555555556;
x270.l = 0.888888888888889;
x271.l = 0.0555555555555556;
x272.l = 0.0555555555555556;
x273.l = 0.888888888888889;
x274.l = 0.0555555555555556;
x275.l = 0.0555555555555556;
x276.l = 0.888888888888889;
x277.l = 0.0555555555555556;
x278.l = 0.0555555555555556;
x279.l = 0.888888888888889;
x280.l = 0.0555555555555556;
x281.l = 0.0555555555555556;
x282.l = 0.888888888888889;
x283.l = 0.0555555555555556;
x284.l = 0.0555555555555556;
x285.l = 0.888888888888889;
x286.l = 0.0555555555555556;
x287.l = 0.0555555555555556;
x288.l = 0.888888888888889;
x289.l = 0.0555555555555556;
x290.l = 0.0555555555555556;
x291.l = 0.888888888888889;
x292.l = 0.0555555555555556;
x293.l = 0.0555555555555556;
x294.l = 0.888888888888889;
x295.l = 0.0555555555555556;
x296.l = 0.0555555555555556;
x297.l = 0.888888888888889;
x298.l = 0.0555555555555556;
x299.l = 0.0555555555555556;
x300.l = 0.888888888888889;
x301.l = 0.0555555555555556;
x302.l = 0.0555555555555556;
x303.l = 0.888888888888889;
x304.l = 0.0555555555555556;
x305.l = 0.0555555555555556;
x306.l = 0.888888888888889;
x307.l = 0.0555555555555556;
x308.l = 0.0555555555555556;
x309.l = 0.888888888888889;
x310.l = 0.0555555555555556;
x311.l = 0.0555555555555556;
x312.l = 0.888888888888889;
x313.l = 0.0555555555555556;
x314.l = 0.0555555555555556;
x315.l = 0.888888888888889;
x316.l = 0.0555555555555556;

Model m / all /;

m.limrow=0; m.limcol=0;
m.tolproj=0.0;

$if NOT '%gams.u1%' == '' $include '%gams.u1%'

$if not set NLP $set NLP NLP
Solve m using %NLP% minimizing objvar;


Last updated: 2024-03-25 Git hash: 1dae024f
Imprint / Privacy Policy / License: CC-BY 4.0