more tests

Author: joaquimg

src/infeasibility.jl

@@ -317,19 +317,20 @@ function iis_elastic_filter(original_model::JuMP.GenericModel, optimizer)
                 var = collect(keys(func.terms))
                 coef1 = func.terms[var[1]]
                 coef2 = func.terms[var[2]]
-                if JuMP.value(var1) > tolerance && JuMP.value(var2) > tolerance
+                if JuMP.value(var[1]) > tolerance &&
+                   JuMP.value(var[2]) > tolerance
                     error("IIS failed due numerical instability")
                 elseif JuMP.value(var[1]) > tolerance
                     has_lower = JuMP.has_lower_bound(var[1])
                     JuMP.fix(var[1], 0.0; force = true)
-                    # or delete(model, var1)
+                    # or delete(model, var[1])
                     delete!(constraint_to_affine, con)
                     constraint_to_affine[con] = coef2 * var[2]
                     push!(de_elastisized, (con, var[1], has_lower))
                 elseif JuMP.value(var[2]) > tolerance
                     has_lower = JuMP.has_lower_bound(var[2])
                     JuMP.fix(var[2], 0.0; force = true)
-                    # or delete(model, var2)
+                    # or delete(model, var[2])
                     delete!(constraint_to_affine, con)
                     constraint_to_affine[con] = coef1 * var[1]
                     push!(de_elastisized, (con, var[2], has_lower))
@@ -572,7 +573,7 @@ function ModelAnalyzer._summarize(
 end
 
 function ModelAnalyzer._summarize(io::IO, issue::IrreducibleInfeasibleSubset)
-    return print(io, "IIS: ", join(map(issue.constraint, _name), ", "))
+    return print(io, "IIS: ", join(map(_name, issue.constraint), ", "))
 end
 
 function ModelAnalyzer._verbose_summarize(
@@ -631,7 +632,7 @@ function ModelAnalyzer._verbose_summarize(
     return print(
         io,
         "Irreducible Infeasible Subset: ",
-        join(map(issue.constraint, _name), ", "),
+        join(map(_name, issue.constraint), ", "),
     )
 end
 

test/infeasibility.jl

@@ -260,6 +260,114 @@ function test_iis()
     @test length(ret[].constraint) == 2
     @test ret[].constraint[1] == c2
     @test ret[].constraint[2] == c1
+    #
+    buf = IOBuffer()
+    ModelAnalyzer.summarize(
+        buf,
+        ModelAnalyzer.Infeasibility.IrreducibleInfeasibleSubset,
+    )
+    str = String(take!(buf))
+    @test startswith(str, "# `IrreducibleInfeasibleSubset`")
+    ModelAnalyzer.summarize(
+        buf,
+        ModelAnalyzer.Infeasibility.IrreducibleInfeasibleSubset,
+        verbose = false,
+    )
+    str = String(take!(buf))
+    @test str == "# IrreducibleInfeasibleSubset"
+    #
+    ModelAnalyzer.summarize(buf, ret[1], verbose = true)
+    str = String(take!(buf))
+    @test startswith(str, "Irreducible Infeasible Subset: ")
+    @test contains(str, ", ")
+    ModelAnalyzer.summarize(buf, ret[1], verbose = false)
+    str = String(take!(buf))
+    @test startswith(str, "IIS: ")
+    @test contains(str, ", ")
+
+    buf = IOBuffer()
+    Base.show(buf, data)
+    str = String(take!(buf))
+    return
+end
+
+function test_iis_multiple()
+    model = Model(HiGHS.Optimizer)
+    set_silent(model)
+    @variable(model, 0 <= x <= 10)
+    @variable(model, 0 <= y <= 20)
+    @constraint(model, c1, x + y <= 1)
+    @constraint(model, c3, x + y <= 1.5)
+    @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,
+        optimizer = HiGHS.Optimizer,
+    )
+    list = ModelAnalyzer.list_of_issue_types(data)
+    @test length(list) == 1
+    ret = ModelAnalyzer.list_of_issues(data, list[1])
+    @test length(ret) == 1
+    @test length(ret[].constraint) == 2
+    @test c2 in Set([ret[].constraint[1], ret[].constraint[2]])
+    @test Set([ret[].constraint[1], ret[].constraint[2]]) ⊆ Set([c3, c2, c1])
+    return
+end
+
+function test_iis_interval_right()
+    model = Model(HiGHS.Optimizer)
+    set_silent(model)
+    @variable(model, 0 <= x <= 10)
+    @variable(model, 0 <= y <= 20)
+    @constraint(model, c1, 0 <= x + y <= 1)
+    @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,
+        optimizer = HiGHS.Optimizer,
+    )
+    list = ModelAnalyzer.list_of_issue_types(data)
+    @test length(list) == 1
+    ret = ModelAnalyzer.list_of_issues(data, list[1])
+    @test length(ret) == 1
+    @test length(ret[].constraint) == 2
+    @test Set([ret[].constraint[1], ret[].constraint[2]]) == Set([c2, c1])
+    return
+end
+
+function test_iis_interval_left()
+    model = Model(HiGHS.Optimizer)
+    set_silent(model)
+    @variable(model, 0 <= x <= 10)
+    @variable(model, 0 <= y <= 20)
+    @constraint(model, c1, x + y <= 1)
+    @constraint(model, c2, 2 <= x + y <= 5)
+    @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,
+        optimizer = HiGHS.Optimizer,
+    )
+    list = ModelAnalyzer.list_of_issue_types(data)
+    @test length(list) == 1
+    ret = ModelAnalyzer.list_of_issues(data, list[1])
+    @test length(ret) == 1
+    @test length(ret[].constraint) == 2
+    @test Set([ret[].constraint[1], ret[].constraint[2]]) == Set([c2, c1])
     return
 end