Pyomo Abstract exemplo 1

Exemplo de modelo abstrato no Pyomo.

Link para download.

A grande vantagem é que separa claramente os dados do modelo.

Uma desvantagem é que o formato de dados tem que carregar os índices das matrizes, de forma até redundante.

Portanto, creio que vale a pena para modelos grandes, onde a chance de infeasible é alta e uma organização melhor do problema se paga.

Dados do arquivo dat:

param : N : c :=
1 1
2 2;

param : M : b :=
1 1
2 2;

param a :=
1 1 3
1 2 4
2 1 2
2 2 5;

Modelo abstrato:

– coding: utf-8 –

“””
Created on Mon Jul 13 05:48:48 2020

@author: asgun
“””

from pyomo.environ import *

model = AbstractModel()

model.N = Set()
model.M = Set()
model.c = Param(model.N)
model.a = Param(model.M, model.N)
model.b = Param(model.M)

model.x = Var(model.N, within= NonNegativeReals)

def obj_rule(model):
return sum(model.c[i]*model.x[i] for i in model.N)
model.obj = Objective(rule = obj_rule)

def con_rule(model, m):
return sum(model.a[m,i]*model.x[i] for i in model.N) >= model.b[m]

model.con = Constraint (model.M, rule = con_rule)

opt = SolverFactory(‘glpk’)

instance = model.create_instance(“Abstract_data.dat”)
results = opt.solve(instance)

instance.display()

Pyomo concrete 02

Este é um exemplo do Pyomo, para um modelo concreto.

Porém, uma vantagem bastante grande é que os dados podem vir em listas e dicionários do Python, e as regras da função objetivo e restrições podem ser escritas como funções do Python.

Como os projetos normalmente serão de tamanho médio, eu acho este o método mais interessante, por manipular via Python o arquivo inteiro.

Projetos muito mais pesados teriam que ser resolvidos com solvers mais poderosos, como o CPLEX, e aí a coisa muda de nível – é melhor usar AIMMS ou algo mais robusto, até para fins de debug.

Link para download.

Código:

from pyomo.environ import *

#dados

N = [1,2]
M = [1,2]

c = {1:1, 2:2}
a = {(1,1):3, (1,2):4, (2,1):2, (2,2):5}
b = {1:1, 2:2}

#Model

model = ConcreteModel()

model.x = Var(N, within = NonNegativeReals)

def obj_rule(model):
return sum(c[i] + model.x[i] for i in N)

model.obj = Objective(rule = obj_rule)

def cons_1(model,m):
return sum(a[m,i]*model.x[i] for i in N) >= b[m]

model.con = Constraint(M, rule = cons_1)

opt = SolverFactory(‘cbc’) #Ou glpk
results = opt.solve(model)

print(results)
model.display()

Pyomo Exemplo 1

Link para download.

Exemplo mais simples possível do Pyomo.

É um modelo concreto, que resolve

Minimizar x1 + 2*x2

s.a.

2*x1+4*x2>=1

2*x1+5*x2>=2

Código:

from pyomo.environ import *

model = ConcreteModel()
model.x1 = Var(within=NonNegativeReals)
model.x2 = Var(within=NonNegativeReals)
model.obj= Objective(expr = model.x1 + 2*model.x2, sense=minimize)

model.con1= Constraint(expr = 2model.x1 +4model.x2 >=1)
model.con2 = Constraint(expr = 2model.x1 + 5model.x2>=2)

opt = SolverFactory(‘cbc’) #Ou glpk
results = opt.solve(model)

print(results)
model.display()