Why do I get "GurobiError: Variable not in model" after using Model.copy()?

Question

我需要优化具有不同约束集的模型，但这些约束的子集对于每个模型都是相同的。我的想法是构建一个包含所有变量和每次都需要的所有约束的基础模型。但是，Gurobi Model.copy() 方法看起来并不像我想象的那样复制变量。这是我希望做的事情的一个非常简化的版本，它显示了我遇到的错误。

>>> from gurobipy import Model
>>> m0 = Model("test")
>>> v = m0.addVar(lb=-1, ub=1)
>>> m0.update()
>>> print m0
<gurobi.Model Continuous instance test: 0 constrs, 1 vars, Parameter changes: LogFile=gurobi.log>
>>> print v in m0.getVars()
True

>>> m = m0.copy()
>>> print m
<gurobi.Model Continuous instance test_copy: 0 constrs, 1 vars, Parameter changes: LogFile=gurobi.log>
>>> print v in m.getVars()
True

>>> m0.addConstr(v <= 0)
<gurobi.Constr *Awaiting Model Update*>
>>> m0.update()
>>> print m0
<gurobi.Model Continuous instance test: 1 constrs, 1 vars, Parameter changes: LogFile=gurobi.log>

>>> m.addConstr(v >= 0)
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
  File "model.pxi", line 2196, in gurobipy.Model.addConstr (../../src/python/gurobipy.c:66304)
  File "model.pxi", line 2089, in gurobipy.Model.__addConstr (../../src/python/gurobipy.c:64663)
gurobipy.GurobiError: Variable not in model

>>> print m
<gurobi.Model Continuous instance test_copy: 0 constrs, 1 vars, Parameter changes: LogFile=gurobi.log>

1. 为什么 m.addConstr(v >= 0) 在 m.getVars() 时找不到变量 v 清楚地显示 v 在模型中？
2. 是否有任何其他方法可以在 为了避免从头开始重建每个模型？

软件：Python 2.7.11 和 Gurobi 6.5 OS X El Capitan（也在 Ubuntu 15.04 上测试）

Answer 1

关于问题 1 为什么 v in m.getVars() returns 在原始示例中为真，我找不到明确的答案。我想 "in" 只是比较对象内的值而不是对象本身，但我没有这方面的来源。

关于问题2，我找到了解决办法：

1. 创建带名称的变量，保存名称。
2. 复制模型。
3. 使用名称映射副本中的变量。

简单示例：

from gurobipy import Model
m0 = Model("test")
v_name = 'v'
v0 = m0.addVar(lb=-1, ub=1, name=v_name)
m0.update()

m = m0.copy()
# returns True
print v0 in m.getVars()

m0.addConstr(v0 <= 0)
m0.update()
# won't work
#m.addConstr(v0 >= 0)
# will work
v = m.getVarByName(v_name)
m.addConstr(v >= 0)
m.update()

更复杂的例子：

from gurobipy import Model
v_name = {}
m0 = Model("test")
v0 = {}
for i in range(10):
    v_name[i] = "v{0}".format(i)
    v0[i] = m0.addVar(lb=-1, ub=1, name=v_name[i])
m0.update()

m = m0.copy()
v = {}
for i in range(10):
    m0.addConstr(v0[i] <= 0)
    v[i] = m.getVarByName(v_name[i])
    m.addConstr(v[i] >= 0)
m0.update()
m.update()

...等等。