如何获取每次迭代后的即时值？Pyomo迭代访问模型值及定期检查方案

1. 获取每次迭代后的即时值 #

不同求解器的实现方式略有差异，以常用的IPOPT为例：

• 利用IPOPT的迭代回调函数，在每次迭代完成后直接读取模型变量的当前值。
• 示例代码：

from pyomo.environ import *
from pyomo.opt import SolverFactory

def ipopt_iteration_callback(alg_mod, iter_count, obj_value, inf_pr, inf_du, mu,
                             d_norm, regularization_size, alpha_du, alpha_pr,
                             ls_trials):
    # 输出当前迭代的变量值
    print(f"第{iter_count}次迭代，变量x的值：{model.x.value}")
    return True  # 返回True继续迭代，返回False终止求解

model = ConcreteModel()
model.x = Var(initialize=1.0)
model.obj = Objective(expr=(model.x - 2)**2)

solver = SolverFactory('ipopt')
solver.options['max_iter'] = 1000
# 注册迭代回调函数
solver.set_callback(ipopt_iteration_callback)
solver.solve(model)

• 其他求解器（如Gurobi、CPLEX）也有类似回调机制，可通过求解器对应接口在回调中直接访问变量的value属性（如model.x.value）获取即时值，具体参数需参考对应求解器的文档。

2. 传入自定义函数访问模型值，每50次迭代检查结果 #

Pyomo完全支持通过求解器回调接口传入自定义函数，实现每50次迭代的检查逻辑：

• 核心思路：在自定义回调函数中维护迭代计数器，当计数器达到50的倍数时执行结果检查。
• 示例代码：

from pyomo.environ import *
from pyomo.opt import SolverFactory

# 初始化迭代计数器
iter_counter = 0

def check_every_50_iterations(alg_mod, iter_count, obj_value, inf_pr, inf_du, mu,
                             d_norm, regularization_size, alpha_du, alpha_pr,
                             ls_trials):
    global iter_counter
    iter_counter += 1
    if iter_counter % 50 == 0:
        print(f"===== 第{iter_counter}次迭代检查 =====")
        print(f"当前目标函数值：{obj_value}")
        print(f"变量x当前值：{model.x.value}")
        # 可添加自定义检查逻辑，比如判断是否满足提前终止条件
        # if abs(obj_value) < 1e-6:
        #     return False  # 返回False触发提前终止
    return True

model = ConcreteModel()
model.x = Var(initialize=1.0)
model.obj = Objective(expr=(model.x - 2)**2)

solver = SolverFactory('ipopt')
solver.options['max_iter'] = 1000
solver.set_callback(check_every_50_iterations)
solver.solve(model)

• 注意事项：
  • 不同求解器的回调函数参数可能不同，需对应调整函数定义。
  • 若不想用全局变量，可通过类封装计数器和检查逻辑，更符合Python编码规范。
  • 部分求解器需开启特定选项启用回调，IPOPT默认支持无需额外配置。

内容的提问来源于stack exchange，提问作者zz z