# NLP written by GAMS Convert at 02/17/22 17:20:16 # # Equation counts # Total E G L N X C B # 64 64 0 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 # # Nonzero counts # Total const NL # 187 97 90 # # Reformulation has removed 1 variable and 1 equation from pyomo.environ import * model = m = ConcreteModel() m.x1 = Var(within=Reals, bounds=(-3,3), initialize=0) m.x2 = Var(within=Reals, bounds=(-3,3), initialize=0) m.x3 = Var(within=Reals, bounds=(-1,1), initialize=0) m.x4 = Var(within=Reals, bounds=(-1,1), initialize=0) m.x5 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x6 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x7 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x8 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x9 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x10 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x11 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x12 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x13 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x14 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x15 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x16 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x17 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x18 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x19 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x20 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x21 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x22 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x23 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x24 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x25 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x26 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x27 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x28 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x29 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x30 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x31 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x32 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x33 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x34 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x35 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x36 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x37 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x38 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x39 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x40 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x41 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x42 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x43 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x44 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x45 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x46 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x47 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x48 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x49 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x50 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x51 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x52 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x53 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x54 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x55 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x56 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x57 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x58 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x59 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x60 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x61 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x62 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x63 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x64 = Var(within=Reals, bounds=(0,1e+07), initialize=0) m.x65 = Var(within=Reals, bounds=(-1e+07,1e+07), initialize=0) m.x66 = Var(within=Reals, bounds=(-1e+07,1e+07), initialize=0) m.obj = Objective(sense=minimize, expr= m.x66) m.e1 = Constraint(expr= 0.166666666666667 * m.x1 - m.x3 == -0.5) m.e2 = Constraint(expr= 0.166666666666667 * m.x2 - m.x4 == -0.5) m.e3 = Constraint(expr= 28.8031207707063 * (0.111229942702413 - m.x3)**2 + 32.7180515537385 * (0.59541072256986 - m.x4)**2 - m.x5 == 0) m.e4 = Constraint(expr= (0.267663478098588 * (2.23606797749979 * sqrt(m.x5) + 1.66666666666667 * m.x5) + 0.267663478098588) * exp(-2.23606797749979 * sqrt(m.x5)) - m.x6 == 0) m.e5 = Constraint(expr= 28.8031207707063 * (0.863858827199432 - m.x3)**2 + 32.7180515537385 * (0.683898783136516 - m.x4)**2 - m.x7 == 0) m.e6 = Constraint(expr= (0.267663478098588 * (2.23606797749979 * sqrt(m.x7) + 1.66666666666667 * m.x7) + 0.267663478098588) * exp(-2.23606797749979 * sqrt(m.x7)) - m.x8 == 0) m.e7 = Constraint(expr= 28.8031207707063 * (0.0951907195096362 - m.x3)**2 + 32.7180515537385 * (0.987881605457913 - m.x4)**2 - m.x9 == 0) m.e8 = Constraint(expr= (0.267663478098588 * (2.23606797749979 * sqrt(m.x9) + 1.66666666666667 * m.x9) + 0.267663478098588) * exp(-2.23606797749979 * sqrt(m.x9)) - m.x10 == 0) m.e9 = Constraint(expr= 28.8031207707063 * (0.730646280982898 - m.x3)**2 + 32.7180515537385 * (0.241742586444198 - m.x4)**2 - m.x11 == 0) m.e10 = Constraint(expr= (0.267663478098588 * (2.23606797749979 * sqrt(m.x11) + 1.66666666666667 * m.x11) + 0.267663478098588) * exp(-2.23606797749979 * sqrt(m.x11)) - m.x12 == 0) m.e11 = Constraint(expr= 28.8031207707063 * (0.320238787378679 - m.x3)**2 + 32.7180515537385 * (0.218296935351353 - m.x4)**2 - m.x13 == 0) m.e12 = Constraint(expr= (0.267663478098588 * (2.23606797749979 * sqrt(m.x13) + 1.66666666666667 * m.x13) + 0.267663478098588) * exp(-2.23606797749979 * sqrt(m.x13)) - m.x14 == 0) m.e13 = Constraint(expr= 28.8031207707063 * (0.87059458573363 - m.x3)**2 + 32.7180515537385 * (0.801085894330831 - m.x4)**2 - m.x15 == 0) m.e14 = Constraint(expr= (0.267663478098588 * (2.23606797749979 * sqrt(m.x15) + 1.66666666666667 * m.x15) + 0.267663478098588) * exp(-2.23606797749979 * sqrt(m.x15)) - m.x16 == 0) m.e15 = Constraint(expr= 28.8031207707063 * (0.0415641459792807 - m.x3)**2 + 32.7180515537385 * (0.659793606686887 - m.x4)**2 - m.x17 == 0) m.e16 = Constraint(expr= (0.267663478098588 * (2.23606797749979 * sqrt(m.x17) + 1.66666666666667 * m.x17) + 0.267663478098588) * exp(-2.23606797749979 * sqrt(m.x17)) - m.x18 == 0) m.e17 = Constraint(expr= 28.8031207707063 * (0.801433475166551 - m.x3)**2 + 32.7180515537385 * (0.402392282692781 - m.x4)**2 - m.x19 == 0) m.e18 = Constraint(expr= (0.267663478098588 * (2.23606797749979 * sqrt(m.x19) + 1.66666666666667 * m.x19) + 0.267663478098588) * exp(-2.23606797749979 * sqrt(m.x19)) - m.x20 == 0) m.e19 = Constraint(expr= 28.8031207707063 * (0.419134454621807 - m.x3)**2 + 32.7180515537385 * (0.520974338874469 - m.x4)**2 - m.x21 == 0) m.e20 = Constraint(expr= (0.267663478098588 * (2.23606797749979 * sqrt(m.x21) + 1.66666666666667 * m.x21) + 0.267663478098588) * exp(-2.23606797749979 * sqrt(m.x21)) - m.x22 == 0) m.e21 = Constraint(expr= 28.8031207707063 * (0.917201570287689 - m.x3)**2 + 32.7180515537385 * (0.84457790520867 - m.x4)**2 - m.x23 == 0) m.e22 = Constraint(expr= (0.267663478098588 * (2.23606797749979 * sqrt(m.x23) + 1.66666666666667 * m.x23) + 0.267663478098588) * exp(-2.23606797749979 * sqrt(m.x23)) - m.x24 == 0) m.e23 = Constraint(expr= 28.8031207707063 * (0.452550488432341 - m.x3)**2 + 32.7180515537385 * (0.465470910341813 - m.x4)**2 - m.x25 == 0) m.e24 = Constraint(expr= (0.267663478098588 * (2.23606797749979 * sqrt(m.x25) + 1.66666666666667 * m.x25) + 0.267663478098588) * exp(-2.23606797749979 * sqrt(m.x25)) - m.x26 == 0) m.e25 = Constraint(expr= 28.8031207707063 * (0.381086148033683 - m.x3)**2 + 32.7180515537385 * (0.0284410867012327 - m.x4)**2 - m.x27 == 0) m.e26 = Constraint(expr= (0.267663478098588 * (2.23606797749979 * sqrt(m.x27) + 1.66666666666667 * m.x27) + 0.267663478098588) * exp(-2.23606797749979 * sqrt(m.x27)) - m.x28 == 0) m.e27 = Constraint(expr= 28.8031207707063 * (0.664445710060874 - m.x3)**2 + 32.7180515537385 * (0.715687166073911 - m.x4)**2 - m.x29 == 0) m.e28 = Constraint(expr= (0.267663478098588 * (2.23606797749979 * sqrt(m.x29) + 1.66666666666667 * m.x29) + 0.267663478098588) * exp(-2.23606797749979 * sqrt(m.x29)) - m.x30 == 0) m.e29 = Constraint(expr= 28.8031207707063 * (0.35783709395474 - m.x3)**2 + 32.7180515537385 * (0.557803528332154 - m.x4)**2 - m.x31 == 0) m.e30 = Constraint(expr= (0.267663478098588 * (2.23606797749979 * sqrt(m.x31) + 1.66666666666667 * m.x31) + 0.267663478098588) * exp(-2.23606797749979 * sqrt(m.x31)) - m.x32 == 0) m.e31 = Constraint(expr= 28.8031207707063 * (0.290737438310084 - m.x3)**2 + 32.7180515537385 * (0.784795849934872 - m.x4)**2 - m.x33 == 0) m.e32 = Constraint(expr= (0.267663478098588 * (2.23606797749979 * sqrt(m.x33) + 1.66666666666667 * m.x33) + 0.267663478098588) * exp(-2.23606797749979 * sqrt(m.x33)) - m.x34 == 0) m.e33 = Constraint(expr= 28.8031207707063 * (0.184545159482649 - m.x3)**2 + 32.7180515537385 * (0.071764872448536 - m.x4)**2 - m.x35 == 0) m.e34 = Constraint(expr= (0.267663478098588 * (2.23606797749979 * sqrt(m.x35) + 1.66666666666667 * m.x35) + 0.267663478098588) * exp(-2.23606797749979 * sqrt(m.x35)) - m.x36 == 0) m.e35 = Constraint(expr= 28.8031207707063 * (0.52271755235026 - m.x3)**2 + 32.7180515537385 * (0.337358966433907 - m.x4)**2 - m.x37 == 0) m.e36 = Constraint(expr= (0.267663478098588 * (2.23606797749979 * sqrt(m.x37) + 1.66666666666667 * m.x37) + 0.267663478098588) * exp(-2.23606797749979 * sqrt(m.x37)) - m.x38 == 0) m.e37 = Constraint(expr= 28.8031207707063 * (0.245628297893699 - m.x3)**2 + 32.7180515537385 * (0.10677814322697 - m.x4)**2 - m.x39 == 0) m.e38 = Constraint(expr= (0.267663478098588 * (2.23606797749979 * sqrt(m.x39) + 1.66666666666667 * m.x39) + 0.267663478098588) * exp(-2.23606797749979 * sqrt(m.x39)) - m.x40 == 0) m.e39 = Constraint(expr= 28.8031207707063 * (0.948138629452501 - m.x3)**2 + 32.7180515537385 * (0.382163884479295 - m.x4)**2 - m.x41 == 0) m.e40 = Constraint(expr= (0.267663478098588 * (2.23606797749979 * sqrt(m.x41) + 1.66666666666667 * m.x41) + 0.267663478098588) * exp(-2.23606797749979 * sqrt(m.x41)) - m.x42 == 0) m.e41 = Constraint(expr= 28.8031207707063 * (0.556510835307559 - m.x3)**2 + 32.7180515537385 * (0.272182954878091 - m.x4)**2 - m.x43 == 0) m.e42 = Constraint(expr= (0.267663478098588 * (2.23606797749979 * sqrt(m.x43) + 1.66666666666667 * m.x43) + 0.267663478098588) * exp(-2.23606797749979 * sqrt(m.x43)) - m.x44 == 0) m.e43 = Constraint(expr= 28.8031207707063 * (0.022380181406173 - m.x3)**2 + 32.7180515537385 * (0.610730424339593 - m.x4)**2 - m.x45 == 0) m.e44 = Constraint(expr= (0.267663478098588 * (2.23606797749979 * sqrt(m.x45) + 1.66666666666667 * m.x45) + 0.267663478098588) * exp(-2.23606797749979 * sqrt(m.x45)) - m.x46 == 0) m.e45 = Constraint(expr= 28.8031207707063 * (0.497905175502968 - m.x3)**2 + 32.7180515537385 * (0.0602044467956143 - m.x4)**2 - m.x47 == 0) m.e46 = Constraint(expr= (0.267663478098588 * (2.23606797749979 * sqrt(m.x47) + 1.66666666666667 * m.x47) + 0.267663478098588) * exp(-2.23606797749979 * sqrt(m.x47)) - m.x48 == 0) m.e47 = Constraint(expr= 28.8031207707063 * (0.792485426296966 - m.x3)**2 + 32.7180515537385 * (0.151938325910849 - m.x4)**2 - m.x49 == 0) m.e48 = Constraint(expr= (0.267663478098588 * (2.23606797749979 * sqrt(m.x49) + 1.66666666666667 * m.x49) + 0.267663478098588) * exp(-2.23606797749979 * sqrt(m.x49)) - m.x50 == 0) m.e49 = Constraint(expr= 28.8031207707063 * (0.670655556838736 - m.x3)**2 + 32.7180515537385 * (0.876936461235205 - m.x4)**2 - m.x51 == 0) m.e50 = Constraint(expr= (0.267663478098588 * (2.23606797749979 * sqrt(m.x51) + 1.66666666666667 * m.x51) + 0.267663478098588) * exp(-2.23606797749979 * sqrt(m.x51)) - m.x52 == 0) m.e51 = Constraint(expr= 28.8031207707063 * (0.207973182236053 - m.x3)**2 + 32.7180515537385 * (0.331173292175481 - m.x4)**2 - m.x53 == 0) m.e52 = Constraint(expr= (0.267663478098588 * (2.23606797749979 * sqrt(m.x53) + 1.66666666666667 * m.x53) + 0.267663478098588) * exp(-2.23606797749979 * sqrt(m.x53)) - m.x54 == 0) m.e53 = Constraint(expr= 28.8031207707063 * (0.966799636494498 - m.x3)**2 + 32.7180515537385 * (0.903623280335378 - m.x4)**2 - m.x55 == 0) m.e54 = Constraint(expr= (0.267663478098588 * (2.23606797749979 * sqrt(m.x55) + 1.66666666666667 * m.x55) + 0.267663478098588) * exp(-2.23606797749979 * sqrt(m.x55)) - m.x56 == 0) m.e55 = Constraint(expr= 28.8031207707063 * (0.73547573937596 - m.x3)**2 + 32.7180515537385 * (0.750945430575571 - m.x4)**2 - m.x57 == 0) m.e56 = Constraint(expr= (0.267663478098588 * (2.23606797749979 * sqrt(m.x57) + 1.66666666666667 * m.x57) + 0.267663478098588) * exp(-2.23606797749979 * sqrt(m.x57)) - m.x58 == 0) m.e57 = Constraint(expr= 28.8031207707063 * (0.610566185432721 - m.x3)**2 + 32.7180515537385 * (0.198845982666349 - m.x4)**2 - m.x59 == 0) m.e58 = Constraint(expr= (0.267663478098588 * (2.23606797749979 * sqrt(m.x59) + 1.66666666666667 * m.x59) + 0.267663478098588) * exp(-2.23606797749979 * sqrt(m.x59)) - m.x60 == 0) m.e59 = Constraint(expr= 28.8031207707063 * (0.141944660623143 - m.x3)**2 + 32.7180515537385 * (0.489288865821674 - m.x4)**2 - m.x61 == 0) m.e60 = Constraint(expr= (0.267663478098588 * (2.23606797749979 * sqrt(m.x61) + 1.66666666666667 * m.x61) + 0.267663478098588) * exp(-2.23606797749979 * sqrt(m.x61)) - m.x62 == 0) m.e61 = Constraint(expr= 28.8031207707063 * (0.578558541443415 - m.x3)**2 + 32.7180515537385 * (0.962551216897219 - m.x4)**2 - m.x63 == 0) m.e62 = Constraint(expr= (0.267663478098588 * (2.23606797749979 * sqrt(m.x63) + 1.66666666666667 * m.x63) + 0.267663478098588) * exp(-2.23606797749979 * sqrt(m.x63)) - m.x64 == 0) m.e63 = Constraint(expr= -0.141170866470016 * m.x6 - 0.0410227184433634 * m.x8 + 0.00923158102538404 * m.x10 + 0.0322471218273555 * m.x12 + 0.159305500449696 * m.x14 - 0.123117661099077 * m.x16 + 0.0380530713899642 * m.x18 + 0.889594746101726 * m.x20 + 0.514139990210469 * m.x22 - 0.0738007765330698 * m.x24 + 1.44710982730193 * m.x26 + 0.132568049860004 * m.x28 + 1.32552100239566 * m.x30 - 0.880281256081741 * m.x32 + 0.67317507860953 * m.x34 + 0.0699447658622974 * m.x36 - 0.180419508821865 * m.x38 + 0.114486082175227 * m.x40 + 0.0157657511439586 * m.x42 - 2.42272624857116 * m.x44 + 0.0619293699366356 * m.x46 - 0.102616266190338 * m.x48 + 0.236512593081257 * m.x50 + 0.39189121299016 * m.x52 + 0.0949736077002411 * m.x54 - 0.013659007622883 * m.x56 + 0.408191007714036 * m.x58 - 1.79436698216086 * m.x60 - 0.393018856049253 * m.x62 + 0.140369132198771 * m.x64 - m.x65 == 0) m.e64 = Constraint(expr= 1.70327473357547 * m.x65 - m.x66 == 0.104337166479966)