The Mixed-Integer Nonlinear Decomposition Toolbox in Pyomo (MindtPy) solver allows users to solve Mixed-Integer Nonlinear Programs (MINLP) using decomposition algorithms. These decomposition algorithms usually rely on the solution of Mixed-Intger Linear Programs (MILP) and Nonlinear Programs (NLP).
MindtPy currently implements the Outer Approximation (OA) algorithm originally described in Duran & Grossmann. Usage and implementation details for MindtPy can be found in the PSE 2018 paper Bernal et al., (ref, preprint).
Usage of MindtPy to solve a Pyomo concrete model involves:
An example which includes the modeling approach may be found below.
Required imports >>> from pyomo.environ import * Create a simple model >>> model = ConcreteModel() >>> model.x = Var(bounds=(1.0,10.0),initialize=5.0) >>> model.y = Var(within=Binary) >>> model.c1 = Constraint(expr=(model.x-3.0)**2 <= 50.0*(1-model.y)) >>> model.c2 = Constraint(expr=model.x*log(model.x)+5.0 <= 50.0*(model.y)) >>> model.objective = Objective(expr=model.x, sense=minimize) Solve the model using MindtPy >>> SolverFactory('mindtpy').solve(model, mip_solver='glpk', nlp_solver='ipopt')
The solution may then be displayed by using the commands
>>> model.objective.display() >>> model.display() >>> model.pprint()
When troubleshooting, it can often be helpful to turn on verbose
output using the
>>> SolverFactory('mindtpy').solve(model, mip_solver='glpk', nlp_solver='ipopt', tee=True)
MindtPy implementation and optional arguments¶
MindtPy optional arguments should be considered beta code and are subject to change.
A decomposition-based MINLP solver.
Check if solver is available. TODO: For now, it is always available. However, sub-solvers may not always be available, and so this should reflect that possibility.
Solve the model. Warning: this solver is still in beta. Keyword arguments subject to change. Undocumented keyword arguments definitely subject to change. Warning: at this point in time, if you try to use PSC or GBD with anything other than IPOPT as the NLP solver, bad things will happen. This is because the suffixes are not in place to extract dual values from the variable bounds for any other solver. TODO: fix needed with the GBD implementation. :param model: a Pyomo model or block to be solved :type model: Block
Return a 3-tuple describing the solver version.