# MINLP written by GAMS Convert at 02/17/22 17:20:42 # # Equation counts # Total E G L N X C B # 9 0 9 0 0 0 0 0 # # Variable counts # x b i s1s s2s sc si # Total cont binary integer sos1 sos2 scont sint # 9 0 0 9 0 0 0 0 # FX 0 # # Nonzero counts # Total const NL # 81 0 81 # # Reformulation has removed 1 variable and 1 equation from pyomo.environ import * model = m = ConcreteModel() m.i1 = Var(within=Integers, bounds=(0,200), initialize=100) m.i2 = Var(within=Integers, bounds=(0,200), initialize=100) m.i3 = Var(within=Integers, bounds=(0,200), initialize=100) m.i4 = Var(within=Integers, bounds=(0,200), initialize=100) m.i5 = Var(within=Integers, bounds=(0,200), initialize=100) m.i6 = Var(within=Integers, bounds=(0,200), initialize=100) m.i7 = Var(within=Integers, bounds=(0,200), initialize=100) m.i8 = Var(within=Integers, bounds=(0,200), initialize=100) m.i9 = Var(within=Integers, bounds=(0,200), initialize=100) m.obj = Objective(sense=minimize, expr= 7 * m.i1**2 + 6 * m.i2**2 + 24.4 * m.i1 - 0.2 * m.i2 + 8 * m.i3**2 - 6 * m.i3 * m.i1 + 4 * m.i3 * m.i2 + m.i3 + 6 * m.i4**2 + 2 * m.i4 * m.i1 + 2 * m.i4 * m.i3 - 39.2 * m.i4 + 7 * m.i5**2 - 4 * m.i5 * m.i1 - 2 * m.i5 * m.i2 - 6 * m.i5 * m.i3 - 118.4 * m.i5 + 4 * m.i6**2 + 2 * m.i6 * m.i1 - 4 * m.i6 * m.i2 - 4 * m.i6 * m.i3 - 2 * m.i6 * m.i4 + 6 * m.i6 * m.i5 - 73 * m.i6 + 6 * m.i7**2 - 2 * m.i7 * m.i1 - 6 * m.i7 * m.i2 - 2 * m.i7 * m.i3 + 4 * m.i7 * m.i5 + 4 * m.i7 * m.i6 - 110.8 * m.i7 + 7 * m.i8**2 - 4 * m.i8 * m.i1 - 2 * m.i8 * m.i2 + 6 * m.i8 * m.i3 + 4 * m.i8 * m.i4 - 4 * m.i8 * m.i5 - 2 * m.i8 * m.i6 + 4 * m.i8 * m.i7 - 17.8 * m.i8 + 8 * m.i9**2 - 2 * m.i9 * m.i1 - 4 * m.i9 * m.i2 + 4 * m.i9 * m.i3 + 4 * m.i9 * m.i4 - 4 * m.i9 * m.i5 - 4 * m.i9 * m.i6 + 8 * m.i9 * m.i7 + 4 * m.i9 * m.i8 - 29.4 * m.i9) m.e1 = Constraint(expr= -9 * m.i1**2 - 10 * m.i1 * m.i2 - 8 * m.i2**2 - 5 * m.i3**2 - 6 * m.i3 * m.i1 - 10 * m.i3 * m.i2 - 7 * m.i4**2 - 10 * m.i4 * m.i1 - 6 * m.i4 * m.i2 - 2 * m.i4 * m.i3 - 2 * m.i5 * m.i2 - 7 * m.i5**2 - 6 * m.i6 * m.i1 - 2 * m.i6 * m.i2 - 2 * m.i6 * m.i4 - 5 * m.i6**2 + 6 * m.i7 * m.i1 + 2 * m.i7 * m.i2 + 4 * m.i7 * m.i3 + 2 * m.i7 * m.i4 - 4 * m.i7 * m.i5 + 4 * m.i7 * m.i6 - 8 * m.i7**2 - 2 * m.i8 * m.i1 - 8 * m.i8 * m.i2 - 2 * m.i8 * m.i3 + 6 * m.i8 * m.i5 - 2 * m.i8 * m.i7 - 6 * m.i8**2 + 2 * m.i9 * m.i3 - 4 * m.i9 * m.i4 + 8 * m.i9 * m.i5 + 4 * m.i9 * m.i6 - 6 * m.i9 * m.i8 - 6 * m.i9**2 >= -1850) m.e2 = Constraint(expr= -6 * m.i1**2 - 8 * m.i1 * m.i2 - 6 * m.i2**2 - 4 * m.i3 **2 - 2 * m.i3 * m.i1 - 2 * m.i3 * m.i2 - 8 * m.i4**2 + 2 * m.i4 * m.i1 + 10 * m.i4 * m.i2 - 2 * m.i5 * m.i1 - 6 * m.i5 * m.i2 + 6 * m.i5 * m.i4 + 7 * m.i5**2 - 2 * m.i6 * m.i2 + 8 * m.i6 * m.i3 + 2 * m.i6 * m.i4 - 4 * m.i6 * m.i5 - 8 * m.i6**2 - 6 * m.i7 * m.i1 - 10 * m.i7 * m.i2 - 2 * m.i7 * m.i3 + 10 * m.i7 * m.i4 - 10 * m.i7 * m.i5 - 8 * m.i7**2 - 2 * m.i8 * m.i1 - 4 * m.i8 * m.i2 - 2 * m.i8 * m.i3 - 8 * m.i8 * m.i5 - 8 * m.i8 * m.i7 - 5 * m.i8**2 - 2 * m.i9 * m.i1 - 2 * m.i9 * m.i2 + 4 * m.i9 * m.i6 + 2 * m.i9 * m.i7 - 6 * m.i9**2 >= -3170) m.e3 = Constraint(expr= -9 * m.i1**2 - 6 * m.i2**2 - 8 * m.i3**2 + 2 * m.i2 * m.i1 + 2 * m.i3 * m.i2 - 6 * m.i4**2 + 4 * m.i4 * m.i1 + 4 * m.i4 * m.i2 - 2 * m.i4 * m.i3 - 6 * m.i5 * m.i1 - 2 * m.i5 * m.i2 + 4 * m.i5 * m.i4 + 6 * m.i5**2 + 2 * m.i6 * m.i1 + 4 * m.i6 * m.i2 - 6 * m.i6 * m.i4 - 2 * m.i6 * m.i5 - 5 * m.i6**2 + 2 * m.i7 * m.i2 - 4 * m.i7 * m.i3 - 6 * m.i7 * m.i5 - 4 * m.i7 * m.i6 - 7 * m.i7**2 - 2 * m.i8 * m.i1 + 4 * m.i8 * m.i3 + 2 * m.i8 * m.i4 - 4 * m.i8**2 + 10 * m.i9 * m.i1 + 6 * m.i9 * m.i2 - 4 * m.i9 * m.i3 - 10 * m.i9 * m.i4 + 8 * m.i9 * m.i5 - 6 * m.i9 * m.i6 - 2 * m.i9 * m.i7 - 8 * m.i9**2 >= -1770) m.e4 = Constraint(expr= -8 * m.i1**2 - 4 * m.i2**2 - 9 * m.i3**2 - 7 * m.i4**2 - 2 * m.i2 * m.i1 - 2 * m.i3 * m.i1 - 4 * m.i3 * m.i2 + 6 * m.i4 * m.i1 + 2 * m.i4 * m.i2 - 2 * m.i4 * m.i3 - 6 * m.i5 * m.i1 - 4 * m.i5 * m.i2 - 2 * m.i5 * m.i3 + 6 * m.i5 * m.i4 + 6 * m.i5**2 - 10 * m.i6 * m.i1 - 10 * m.i6 * m.i3 + 4 * m.i6 * m.i4 - 2 * m.i6 * m.i5 - 7 * m.i6**2 + 6 * m.i7 * m.i1 - 2 * m.i7 * m.i2 - 2 * m.i7 * m.i3 + 6 * m.i7 * m.i5 + 2 * m.i7 * m.i6 - 6 * m.i7**2 + 4 * m.i8 * m.i1 - 4 * m.i8 * m.i2 + 2 * m.i8 * m.i3 - 4 * m.i8 * m.i4 - 4 * m.i8 * m.i5 + 8 * m.i8 * m.i6 + 6 * m.i8 * m.i6 - 8 * m.i8**2 - 4 * m.i9 * m.i1 + 4 * m.i9 * m.i2 + 6 * m.i9 * m.i3 - 2 * m.i9 * m.i4 + 2 * m.i9 * m.i6 + 8 * m.i9 * m.i7 - 4 * m.i9 * m.i8 - 10 * m.i9**2 >= -1460) m.e5 = Constraint(expr= 2 * m.i2 * m.i1 - 4 * m.i1**2 - 5 * m.i2**2 - 6 * m.i3 * m.i1 - 8 * m.i3**2 - 2 * m.i4 * m.i1 + 6 * m.i4 * m.i2 - 2 * m.i4 * m.i3 - 6 * m.i4**2 - 4 * m.i5 * m.i1 + 2 * m.i5 * m.i2 - 6 * m.i5 * m.i3 - 8 * m.i5 * m.i4 - 7 * m.i5**2 + 4 * m.i6 * m.i1 - 4 * m.i6 * m.i2 + 6 * m.i6 * m.i3 + 4 * m.i6 * m.i5 - 7 * m.i6**2 + 4 * m.i7 * m.i1 - 4 * m.i7 * m.i2 - 4 * m.i7 * m.i3 + 4 * m.i7 * m.i4 + 4 * m.i7 * m.i5 + 4 * m.i7 * m.i6 - 8 * m.i7**2 - 2 * m.i8 * m.i1 + 4 * m.i8 * m.i4 + 2 * m.i8 * m.i6 + 2 * m.i8 * m.i7 - 4 * m.i8**2 - 2 * m.i9 * m.i2 + 4 * m.i9 * m.i3 + 4 * m.i9 * m.i4 - 2 * m.i9 * m.i5 + 2 * m.i9 * m.i6 + 6 * m.i9 * m.i7 - 6 * m.i9 * m.i8 - 7 * m.i9**2 >= -1140) m.e6 = Constraint(expr= 2 * m.i2 * m.i1 - 7 * m.i1**2 - 7 * m.i2**2 - 6 * m.i3 * m.i1 - 2 * m.i3 * m.i2 - 6 * m.i3**2 - 2 * m.i4 * m.i1 + 2 * m.i4 * m.i2 - 2 * m.i4 * m.i3 - 5 * m.i4**2 - 2 * m.i5 * m.i1 - 4 * m.i5 * m.i3 + 2 * m.i5 * m.i4 - 5 * m.i5**2 + 2 * m.i6 * m.i1 - 4 * m.i6 * m.i2 + 4 * m.i6 * m.i3 + 2 * m.i6 * m.i4 + 6 * m.i6 * m.i5 - 9 * m.i6**2 + 4 * m.i7 * m.i2 - 4 * m.i7 * m.i3 + 4 * m.i7 * m.i4 - 4 * m.i7 * m.i5 + 8 * m.i7 * m.i6 - 5 * m.i7**2 + 4 * m.i8 * m.i1 + 8 * m.i8 * m.i2 + 2 * m.i8 * m.i3 - 4 * m.i8 * m.i4 - 2 * m.i8 * m.i5 + 4 * m.i8 * m.i6 - 9 * m.i8**2 - 4 * m.i9 * m.i1 + 2 * m.i9 * m.i4 + 6 * m.i9 * m.i5 - 4 * m.i9 * m.i6 - 2 * m.i9 * m.i7 + 2 * m.i9 * m.i8 - 6 * m.i9**2 >= -940) m.e7 = Constraint(expr= -9 * m.i1**2 - 4 * m.i2 * m.i1 - 8 * m.i2**2 + 4 * m.i3 * m.i1 + 2 * m.i3 * m.i2 - 7 * m.i3**2 + 4 * m.i4 * m.i1 + 4 * m.i4 * m.i3 - 7 * m.i4**2 - 2 * m.i5 * m.i1 - 12 * m.i5 * m.i2 - 4 * m.i5 * m.i3 - 8 * m.i5**2 - 8 * m.i6 * m.i1 + 2 * m.i6 * m.i2 - 2 * m.i6 * m.i5 - 6 * m.i6 **2 - 4 * m.i7 * m.i1 - 6 * m.i7 * m.i2 - 2 * m.i7 * m.i3 + 10 * m.i7 * m.i4 - 2 * m.i7 * m.i5 + 2 * m.i7 * m.i6 - 7 * m.i7**2 - 2 * m.i8 * m.i1 + 2 * m.i8 * m.i2 + 2 * m.i8 * m.i3 + 2 * m.i8 * m.i4 - 6 * m.i8 * m.i6 - 2 * m.i8 * m.i7 - 6 * m.i8**2 + 4 * m.i9 * m.i1 + 2 * m.i9 * m.i2 + 4 * m.i9 * m.i3 + 4 * m.i9 * m.i4 + 2 * m.i9 * m.i5 - 2 * m.i9 * m.i6 - 8 * m.i9**2 >= -2720) m.e8 = Constraint(expr= 4 * m.i2 * m.i1 - 7 * m.i1**2 - 8 * m.i2**2 + 4 * m.i3 * m.i1 - 8 * m.i3**2 + 4 * m.i4 * m.i1 + 8 * m.i4 * m.i2 - 6 * m.i4 * m.i3 - 7 * m.i4**2 - 2 * m.i5 * m.i2 + 2 * m.i5 * m.i4 - 5 * m.i5**2 - 2 * m.i6 * m.i1 - 2 * m.i6 * m.i2 + 4 * m.i6 * m.i4 - 4 * m.i6 * m.i5 - 7 * m.i6**2 - 2 * m.i7 * m.i1 + 8 * m.i7 * m.i2 - 2 * m.i7 * m.i3 - 2 * m.i7 * m.i4 + 6 * m.i7 * m.i5 + 2 * m.i7 * m.i6 - 7 * m.i7**2 + 2 * m.i8 * m.i1 - 6 * m.i8 * m.i2 + 6 * m.i8 * m.i3 + 4 * m.i8 * m.i4 + 2 * m.i8 * m.i5 - 4 * m.i8 * m.i6 - 6 * m.i8**2 + 4 * m.i9 * m.i1 - 6 * m.i9 * m.i2 + 2 * m.i9 * m.i3 - 2 * m.i9 * m.i4 + 2 * m.i9 * m.i5 + 6 * m.i9 * m.i6 + 2 * m.i9 * m.i7 - 4 * m.i9 * m.i8 - 6 * m.i9**2 >= -870) m.e9 = Constraint(expr= 2 * m.i2 * m.i1 - 4 * m.i1**2 - 7 * m.i2**2 + 8 * m.i3 * m.i1 - 4 * m.i3 * m.i2 - 9 * m.i3**2 - 2 * m.i4 * m.i1 - 4 * m.i4 * m.i2 - 2 * m.i4 * m.i3 - 6 * m.i4**2 + 4 * m.i5 * m.i1 + 2 * m.i5 * m.i2 + 4 * m.i5 * m.i3 + 6 * m.i5 * m.i4 - 6 * m.i5**2 + 4 * m.i6 * m.i3 - 6 * m.i6 * m.i4 - 7 * m.i6**2 - 2 * m.i7 * m.i2 - 4 * m.i7 * m.i3 + 4 * m.i7 * m.i5 + 8 * m.i7 * m.i6 - 7 * m.i7**2 + 2 * m.i8 * m.i2 - 4 * m.i8 * m.i3 + 2 * m.i8 * m.i4 + 2 * m.i8 * m.i5 + 6 * m.i8 * m.i7 - 7 * m.i8**2 + 4 * m.i9 * m.i1 + 2 * m.i9 * m.i2 - 10 * m.i9 * m.i3 + 2 * m.i9 * m.i5 + 2 * m.i9 * m.i6 - 8 * m.i9 * m.i8 - 6 * m.i9**2 >= -670)