Pyomo结合Xpress求解器热启动实现及tee日志解读问题

Xpress的运行日志结构清晰，你可以通过以下特征快速判断热启动是否生效，以及求解进度：

热启动生效标识 #

日志开头如果出现*Loading user MIP start solution*、*Initial feasible solution objective: [你的初始解目标值]*这类输出，说明你传入的初始值已经被求解器成功接收。如果没有上述输出，说明初始值根本没进入求解器。

常规MIP求解日志字段说明 #

默认日志每行按顺序对应以下字段：

• Node：当前分支定界搜索树已探索的节点数，为0时代表仍在处理根节点
• BestSol：当前搜索到的最优可行解目标值，最小化问题无可行解时显示inf，最大化问题无可行解时显示-inf。如果热启动生效，该行第一行就会显示你传入的初始解目标值，而非无穷值
• BestBnd：当前全局最优理论边界（最小化问题为下界，最大化问题为上界）
• Gap：当前可行解与理论边界的相对间隙，计算公式为|BestSol - BestBnd| / |BestSol| * 100%，热启动生效时初始间隙不会是100%或无穷大
• Itt：累计单纯形迭代次数
• NodeInf：当前节点的不可行约束数量
• Time：已运行求解时间

你当前的代码只开启了warmstart=True参数，但未完成前置的变量赋值和求解器参数配置，导致初始值没有被传递到Xpress内核，正确操作流程如下：

步骤1：上一次求解完成后保存变量值 #

确保你保存的是所有激活变量的取值，变量名是后续匹配的唯一标识：

import pickle
from pyomo.environ import Var

# 求解完成后执行，保存所有激活变量的取值到字典
var_value_cache = {var.name: var.value for var in model.component_data_objects(Var, active=True)}
# 持久化到本地文件，下次运行直接读取
with open("var_cache.pkl", "wb") as f:
    pickle.dump(var_value_cache, f)

步骤2：新模型构建完成后加载初始值 #

要保证新模型的变量名称、维度和旧模型完全一致，否则对应变量的初始值会无法匹配：

# 读取缓存的变量值
with open("var_cache.pkl", "rb") as f:
    var_value_cache = pickle.load(f)

# 给新模型的变量赋值
for var in model.component_data_objects(Var, active=True):
    if var.name in var_value_cache:
        var.set_value(var_value_cache[var.name], skip_validation=True)

步骤3：正确配置求解器参数调用 #

仅传warmstart=True在部分Pyomo版本中不会把参数传递到Xpress内核，需要手动配置Xpress原生参数强制启用初始解读取：

from pyomo.environ import SolverFactory

opt = SolverFactory('xpress_direct')
# 强制启用用户初始解
opt.options['mipstart'] = 1
# 初始解处理模式：1=直接使用（不可行则丢弃），2=尝试修复不可行的初始解
opt.options['mipstartmode'] = 1

results = opt.solve(model, warmstart=True, tee=True)

常见失效原因排查 #

• 两次运行的模型结构不一致：变量增删、约束修改都会导致初始值无法匹配，Xpress会自动丢弃无效的初始解
• 初始解不可行：日志会输出*MIP start solution is infeasible, discarded*，可以把mipstartmode设为2让求解器尝试修复
• Pyomo版本过低：5.7版本以前的Pyomo对xpress_direct的热启动支持存在bug，升级到最新稳定版即可解决

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