使用 JumP 获得不可约不一致子系统 (IIS) - Julia

Question

我正在尝试获取 IIS 中包含的约束，例如以下不可行的问题

using JuMP, Gurobi
model = direct_model(Gurobi.Optimizer())

# defining variables
@variable(model, z)
@variable(model, x1[a=1:10, b=1:5] => 0)
@variable(model, x2[a=1:10, b=1:5, c=1:7] => 0)
# and so on...

# objective and constraints

@objective(model, Max, z)

@constraint(model, const1[a=1:10, b=1:5, c=1:7], equation1)
@constraint(model, const2[b=1:5, c=1:7], equation2)
# and so on...

optimize!(model)

if termination_status(model) == MOI.INFEASIBLE_OR_UNBOUNDED
    @assert termination_status(model) == MOI.INFEASIBLE_OR_UNBOUNDED
    compute_conflict!(model)
    MOI.get(model, MOI.ConstraintConflictStatus(), const1)
end

但是，它总是返回错误

MethodError: no method matching get(::Model, ::MathOptInterface.ConstraintConflictStatus, ::Array{ConstraintRef{Model,MathOptInterface.ConstraintIndex{MathOptInterface.ScalarAffineFunction{Float64},MathOptInterface.EqualTo{Float64}},ScalarShape},3})

我这样做正确吗？

我正在使用 JuliaPro_v1.53-1、JuMP v0.21.6、Gurobi v0.9.11

如有任何帮助或建议，我们将不胜感激。

提前致谢

Answer 1

在社区论坛上回答：https://discourse.julialang.org/t/obtaining-irreducible-inconcient-subsystem-iis-with-jump-gurobi/58491/2

const1

是一个数组。使用广播：

MOI.get.(model, MOI.ConstraintConflictStatus(), const1)

Answer 2

为了快速检查，您还可以将 IIS 打印到文件中，如下所示：

GRBwrite(backend(model), "conflict.ilp")