Pyomo中model.clone()行为疑问及模型安全性验证咨询

我需要克隆已构建的Pyomo模型，但测试后发现model.clone()的行为不符合预期——部分变量被克隆，部分未被克隆。以下是测试代码：

from pyomo.environ import *
from pyomo import *
cplex = SolverFactory('cplex')
model2 = ConcreteModel()
model3 = ConcreteModel()
model2.x = Var(within=Reals, bounds=(-10, 0))
model3.y = Var(within=Reals, bounds=(2, 3))
model3.obj = Objective(expr=(model2.x - model3.y) ** 2, sense=minimize)
model4 = model3.clone()
cplex.solve(model4)

assert model3.y.value is None 
assert abs(model4.y.value - 2) < 1e-8
assert abs(model2.x.value - 0) < 1e-8

疑问与解答 #

• 当被克隆模型引用其他模型的变量时，clone()函数是否会创建新变量？
  clone()不会为外部模型的变量创建新副本。你的代码中model3的目标函数引用了model2.x，克隆model3得到model4后，model4的目标仍然指向原始的model2.x变量实例。这就是求解model4后model2.x的值也被更新的原因——两者共享同一个变量。

• 若模型的约束和目标涉及的所有变量均属于该模型本身，是否可以认为clone()是安全的？
  完全安全。这种情况下，克隆得到的新模型会包含原模型所有组件（变量、约束、目标等）的独立副本，新模型与原模型完全隔离，不会出现变量值互相干扰的情况。

• Pyomo中是否存在可针对特定模型验证该判断的函数？
  Pyomo没有内置的专门函数，但可以自行实现一个检查函数，遍历模型的目标和约束，验证其中的变量是否都属于当前模型：

def model_uses_only_own_vars(model):
    # 检查所有激活的目标函数
    for obj in model.component_objects(Objective, active=True):
        for var in obj.expr.variables():
            if var.model() is not model:
                return False
    # 检查所有激活的约束
    for con in model.component_objects(Constraint, active=True):
        for var in con.expr.variables():
            if var.model() is not model:
                return False
    return True

调用该函数时，若返回True，说明模型所有组件仅引用自身变量，此时使用clone()不会有混淆风险。

内容的提问来源于stack exchange，提问作者若秋雨