variables - 找出Pyomo模型不可行的原因

Question

我得到了一个包含大量变量和约束的 pyomo 具体模型。

不知何故，我的模型中的一个变量违反了一个约束，这使得我的模型不可行:

WARNING: Loading a SolverResults object with a warning status into model=xxxx;
    message from solver=Model was proven to be infeasible.

有没有办法询问求解器不可行的原因？

例如，假设我有一个名为 x 的变量，如果我定义了以下 2 个约束，模型将是 ofc。不可行。

const1:
    x >= 10

const2:
    x <= 5

我想要实现的是指出导致这种不可行性的约束和变量，以便我可以修复它。否则，对于我的大模型，很难知道是什么导致了这种不可行性。

IN: write_some_comment
OUT: variable "x" cannot fulfill "const1" and "const2" at the same time.

Answer 1

许多求解器(包括IPOPT)会在求解器终止时将变量的值交还给您，即使发现问题不可行。那时，您确实有一些选择。

在 pyomo.util.infeasible 中有贡献的代码这可能会帮助你。 https://github.com/Pyomo/pyomo/blob/master/pyomo/util/infeasible.py

用法:

from pyomo.util.infeasible import log_infeasible_constraints
...
SolverFactory('your_solver').solve(model)
...
log_infeasible_constraints(model)