def define_model(**kwds):
model = ConcreteModel()
model.x1 = Var(bounds=(-5,4)) # domain variable
model.x2 = Var(bounds=(-5,4)) # domain variable
model.x3 = Var(bounds=(-5,4)) # domain variable
model.x4 = Var(bounds=(-5,4)) # domain variable
model.x5 = Var(bounds=(-5,4)) # domain variable
model.x6 = Var(bounds=(-5,4)) # domain variable
model.x7 = Var(bounds=(-5,4)) # domain variable
model.Fx1 = Var() # range variable
model.Fx2 = Var() # range variable
model.Fx3 = Var() # range variable
model.Fx4 = Var() # range variable
model.Fx5 = Var() # range variable
model.Fx6 = Var() # range variable
model.Fx7 = Var() # range variable
model.obj = Objective(expr=model.Fx1+model.Fx2+model.Fx3+model.Fx4+model.Fx5+model.Fx6+model.Fx7, sense=kwds.pop('sense',maximize))
model.piecewise1 = Piecewise(model.Fx1,model.x1,
pw_pts=DOMAIN_PTS,
f_rule=F, **kwds)
model.piecewise2 = Piecewise(model.Fx2,model.x2,
pw_pts=DOMAIN_PTS,
f_rule=F, **kwds)
model.piecewise3 = Piecewise(model.Fx3,model.x3,
pw_pts=DOMAIN_PTS,
f_rule=F, **kwds)
model.piecewise4 = Piecewise(model.Fx4,model.x4,
pw_pts=DOMAIN_PTS,
f_rule=F, **kwds)
model.piecewise5 = Piecewise(model.Fx5,model.x5,
pw_pts=DOMAIN_PTS,
f_rule=F, **kwds)
model.piecewise6 = Piecewise(model.Fx6,model.x6,
pw_pts=DOMAIN_PTS,
f_rule=F, **kwds)
model.piecewise7 = Piecewise(model.Fx7,model.x7,
pw_pts=DOMAIN_PTS,
f_rule=F, **kwds)
#Fix the answer for testing purposes
model.set_answer_constraint1 = Constraint(expr= model.x1 == -5.0)
model.set_answer_constraint2 = Constraint(expr= model.x2 == -3.0)
model.set_answer_constraint3 = Constraint(expr= model.x3 == -2.5)
model.set_answer_constraint4 = Constraint(expr= model.x4 == -1.5)
model.set_answer_constraint5 = Constraint(expr= model.x5 == 2.0)
model.set_answer_constraint6 = Constraint(expr= model.x6 == 3.5)
model.set_answer_constraint7 = Constraint(expr= model.x7 == 4.0)
return model
def define_model(**kwds):
model = ConcreteModel()
model.x1 = Var(bounds=(0, 3)) # domain variable
model.x2 = Var(bounds=(0, 3)) # domain variable
model.x3 = Var(bounds=(0, 3)) # domain variable
model.x4 = Var(bounds=(0, 3)) # domain variable
model.Fx1 = Var() # range variable
model.Fx2 = Var() # range variable
model.Fx3 = Var() # range variable
model.Fx4 = Var() # range variable
model.p = Param(initialize=1.0)
model.obj = Objective(expr=model.Fx1 + model.Fx2 + model.Fx3 + model.Fx4 +
model.x1 + model.x2 + model.x3 + model.x4,
sense=kwds.pop('sense', maximize))
model.piecewise1 = Piecewise(model.Fx1,
model.x1,
pw_pts=DOMAIN_PTS,
f_rule=F,
**kwds)
model.piecewise2 = Piecewise(model.Fx2,
model.x2,
pw_pts=DOMAIN_PTS,
f_rule=F,
**kwds)
model.piecewise3 = Piecewise(model.Fx3,
model.x3,
pw_pts=DOMAIN_PTS,
f_rule=F,
**kwds)
model.piecewise4 = Piecewise(model.Fx4,
model.x4,
pw_pts=DOMAIN_PTS,
f_rule=F,
**kwds)
#Fix the answer for testing purposes
model.set_answer_constraint1 = Constraint(
expr=model.x1 == 0.5) # Fx1 should solve to 0
model.set_answer_constraint2 = Constraint(expr=model.x2 == 1.0) #
model.set_answer_constraint3 = Constraint(expr=model.Fx2 == 0.5) #
model.set_answer_constraint4 = Constraint(
expr=model.x3 == 1.5) # Fx3 should solve to 1
model.set_answer_constraint5 = Constraint(
expr=model.x4 == 2.5) # Fx4 should solve to 1.5
return model
