add infeasibility tests

Author: joaquimg

src/infeasibility.jl

@@ -261,7 +261,9 @@ function ModelAnalyzer.analyze(
     end
     iis = iis_elastic_filter(model, optimizer)
     # for now, only one iis is computed
-    push!(out.iis, IrreducibleInfeasibleSubset(iis))
+    if iis !== nothing
+        push!(out.iis, IrreducibleInfeasibleSubset(iis))
+    end
 
     return out
 end
@@ -278,7 +280,7 @@ function iis_elastic_filter(original_model::JuMP.GenericModel, optimizer)
         println(
             "iis resolver cannot continue because model is found to be $(status) by the solver",
         )
-        return
+        return nothing
     end
 
     model, reference_map = JuMP.copy_model(original_model)

test/infeasibility.jl

@@ -34,6 +34,31 @@ function test_bounds()
     @test length(ret) == 1
     @test ret[] ==
           ModelAnalyzer.Infeasibility.InfeasibleBounds{Float64}(y, 2.0, 1.0)
+    #
+    buf = IOBuffer()
+    ModelAnalyzer.summarize(
+        buf,
+        ModelAnalyzer.Infeasibility.InfeasibleBounds{Float64},
+    )
+    str = String(take!(buf))
+    @test startswith(str, "# `InfeasibleBounds`")
+    ModelAnalyzer.summarize(
+        buf,
+        ModelAnalyzer.Infeasibility.InfeasibleBounds{Float64},
+        verbose = false,
+    )
+    str = String(take!(buf))
+    @test str == "# InfeasibleBounds"
+    #
+    ModelAnalyzer.summarize(buf, ret[1], verbose = true)
+    str = String(take!(buf))
+    @test startswith(str, "Variable: ")
+    @test contains(str, " with lower bound ")
+    @test contains(str, " and upper bound ")
+    ModelAnalyzer.summarize(buf, ret[1], verbose = false)
+    str = String(take!(buf))
+    @test contains(str, " : ")
+    @test contains(str, " !<= ")
     return
 end
 
@@ -54,6 +79,33 @@ function test_integrality()
         2.9,
         MOI.Integer(),
     )
+    #
+    buf = IOBuffer()
+    ModelAnalyzer.summarize(
+        buf,
+        ModelAnalyzer.Infeasibility.InfeasibleIntegrality{Float64},
+    )
+    str = String(take!(buf))
+    @test startswith(str, "# `InfeasibleIntegrality`")
+    ModelAnalyzer.summarize(
+        buf,
+        ModelAnalyzer.Infeasibility.InfeasibleIntegrality{Float64},
+        verbose = false,
+    )
+    str = String(take!(buf))
+    @test str == "# InfeasibleIntegrality"
+    #
+    ModelAnalyzer.summarize(buf, ret[1], verbose = true)
+    str = String(take!(buf))
+    @test startswith(str, "Variable: ")
+    @test contains(str, " with lower bound ")
+    @test contains(str, " and upper bound ")
+    @test contains(str, " and integrality constraint: ")
+    ModelAnalyzer.summarize(buf, ret[1], verbose = false)
+    str = String(take!(buf))
+    @test contains(str, " : [")
+    @test contains(str, "; ")
+    @test contains(str, "], ")
     return
 end
 
@@ -95,9 +147,95 @@ function test_range()
         22.0,
         MOI.LessThan{Float64}(1.0),
     )
+    #
+    buf = IOBuffer()
+    ModelAnalyzer.summarize(
+        buf,
+        ModelAnalyzer.Infeasibility.InfeasibleConstraintRange{Float64},
+    )
+    str = String(take!(buf))
+    @test startswith(str, "# `InfeasibleConstraintRange`")
+    ModelAnalyzer.summarize(
+        buf,
+        ModelAnalyzer.Infeasibility.InfeasibleConstraintRange{Float64},
+        verbose = false,
+    )
+    str = String(take!(buf))
+    @test str == "# InfeasibleConstraintRange"
+    #
+    ModelAnalyzer.summarize(buf, ret[1], verbose = true)
+    str = String(take!(buf))
+    @test startswith(str, "Constraint: ")
+    @test contains(str, " with computed lower bound ")
+    @test contains(str, " and computed upper bound ")
+    @test contains(str, " and set: ")
+    ModelAnalyzer.summarize(buf, ret[1], verbose = false)
+    str = String(take!(buf))
+    @test contains(str, " : [")
+    @test contains(str, "; ")
+    @test contains(str, "], !in ")
+    return
+end
+
+function test_range_greaterthan()
+    model = Model()
+    @variable(model, 10 <= x <= 11)
+    @variable(model, 1 <= y <= 11)
+    @constraint(model, c, x + y >= 100)
+    @objective(model, Max, x + y)
+    data = ModelAnalyzer.analyze(ModelAnalyzer.Infeasibility.Analyzer(), model)
+    list = ModelAnalyzer.list_of_issue_types(data)
+    @test length(list) == 1
+    ret = ModelAnalyzer.list_of_issues(data, list[1])
+    @test length(ret) == 1
+    @test ret[] ==
+          ModelAnalyzer.Infeasibility.InfeasibleConstraintRange{Float64}(
+        c,
+        11.0,
+        22.0,
+        MOI.GreaterThan{Float64}(100.0),
+    )
     return
 end
 
+function test_range_equalto()
+    model = Model()
+    @variable(model, 10 <= x <= 11)
+    @variable(model, 1 <= y <= 11)
+    @constraint(model, c, x + y == 100)
+    @objective(model, Max, x + y)
+    data = ModelAnalyzer.analyze(ModelAnalyzer.Infeasibility.Analyzer(), model)
+    list = ModelAnalyzer.list_of_issue_types(data)
+    @test length(list) == 1
+    ret = ModelAnalyzer.list_of_issues(data, list[1])
+    @test length(ret) == 1
+    @test ret[] ==
+          ModelAnalyzer.Infeasibility.InfeasibleConstraintRange{Float64}(
+        c,
+        11.0,
+        22.0,
+        MOI.EqualTo{Float64}(100.0),
+    )
+    return
+end
+
+function test_iis_feasible()
+    model = Model(HiGHS.Optimizer)
+    set_silent(model)
+    @variable(model, 0 <= x <= 10)
+    @variable(model, 0 <= y <= 20)
+    @constraint(model, c1, x + y <= 1)
+    @objective(model, Max, x + y)
+    optimize!(model)
+    data = ModelAnalyzer.analyze(
+        ModelAnalyzer.Infeasibility.Analyzer(),
+        model,
+        optimizer = HiGHS.Optimizer,
+    )
+    list = ModelAnalyzer.list_of_issue_types(data)
+    @test length(list) == 0
+end
+
 function test_iis()
     model = Model(HiGHS.Optimizer)
     set_silent(model)
@@ -107,6 +245,9 @@ function test_iis()
     @constraint(model, c2, x + y >= 2)
     @objective(model, Max, x + y)
     optimize!(model)
+    data = ModelAnalyzer.analyze(ModelAnalyzer.Infeasibility.Analyzer(), model)
+    list = ModelAnalyzer.list_of_issue_types(data)
+    @test length(list) == 0
     data = ModelAnalyzer.analyze(
         ModelAnalyzer.Infeasibility.Analyzer(),
         model,