将索引变量的值设置为索引表达式 - pyomo

Question

我正在研究一个需要我使用python lib编码的项目。称为pyomo。

我需要将索引变量的值设置为索引表达式。

因此，总共有2个模型和1个创建其实例的脚本。

第一种模式：

m.i = [1,2,3] # i=1,2,3 
m.j = [1,2,3] # j=1,2,3 

m.index = Set(within=m.i*m.j) # m.index = [i]*[j] (pyomo set) 

m.a = Var(m.index) # m.a[i,j] 

m.Exp = Expression(m.index) # pyomo expression with [i,j] index 

m.Const = Constraint(m.index, rule=const_rule) 

def const_rule(m, index): 
    return m.a[index] == m.Exp[index] 

第二种模式：

m.index = ...   # same index with first model's m.index 
m.y = Var(m.index) # lets assume it has values or I can put with init command 

的run.py脚本：

instance1 = create_first_model() 
instance2 = create_second model() 

solve(instance2) 

for k in instance1.index:          
    instance1.Exp[k].expr = instance2.y[k].value 

solve(instance1) 

不久：

我需要为这个更好的方法for循环。或更有效的方法（循环很慢...）

问候 奥坎

PS：所有的增值经销商，集，表达和约束的是pyomo对象。

Answer 1

是m.Exp一个包含Pyomo变量的表达式，还是用作一个常量（即一个可变参数）？

假设它是一个表达式，你要根据你的for循环，下面应该工作表达添加约束条件：

def _linking_con(m, k): 
    return m.Exp[k] == value(instance2.y[k]) 
instance1.linking_con = Constraint(instance1.index, rule=_linking_con) 

如果m.Exp是一个恒定值，那我就有一个可变的帕拉姆取代它

m.Exp = Param(m.index, mutable=True) 

，然后改变你的for循环到：

for k in instance1.index: 
    instance1.Exp[k] = value(instance2.y[k]) 

我不认为有一种方法可以避免for循环