Fix issue 287 (#288)

Author: odow

src/MOI_wrapper.jl

@@ -2077,22 +2077,6 @@ function _store_solution(model::Optimizer, ret::HighsInt)
     resize!(x.rowdual, numRows)
     x.model_status = Highs_getModelStatus(model)
     statusP = Ref{HighsInt}()
-    certificates = MOI.get(model, ComputeInfeasibilityCertificate())
-    if certificates && x.model_status == kHighsModelStatusInfeasible
-        ret = Highs_getDualRay(model, statusP, x.rowdual)
-        # Don't `_check_ret(ret)` here, just bail is there isn't a dual ray.
-        x.has_dual_ray = (ret == kHighsStatusOk) && (statusP[] == 1)
-        if x.has_dual_ray
-            _compute_farkas_variable_dual(model, x.coldual)
-        end
-    elseif certificates && x.model_status == kHighsModelStatusUnbounded
-        ret = Highs_getPrimalRay(model, statusP, x.colvalue)
-        # Don't `_check_ret(ret)` here, just bail is there isn't a dual ray.
-        x.has_primal_ray = (ret == kHighsStatusOk) && (statusP[] == 1)
-    end
-    if x.has_dual_ray || x.has_primal_ray
-        return  # If a ray is present, we don't query the solution
-    end
     Highs_getIntInfoValue(model, "primal_solution_status", statusP)
     x.primal_solution_status = statusP[]
     Highs_getIntInfoValue(model, "dual_solution_status", statusP)
@@ -2106,6 +2090,25 @@ function _store_solution(model::Optimizer, ret::HighsInt)
             Highs_getBasis(model, x.colstatus, x.rowstatus)
         end
     end
+    if !MOI.get(model, ComputeInfeasibilityCertificate())
+        return
+    end
+    if x.model_status == kHighsModelStatusInfeasible
+        ret = Highs_getDualRay(model, statusP, x.rowdual)
+        # Don't `_check_ret(ret)` here, just bail is there isn't a dual ray.
+        x.has_dual_ray = (ret == kHighsStatusOk) && (statusP[] == 1)
+        if x.has_dual_ray
+            _compute_farkas_variable_dual(model, x.coldual)
+            x.dual_solution_status = kHighsSolutionStatusNone
+        end
+    elseif x.model_status == kHighsModelStatusUnbounded
+        ret = Highs_getPrimalRay(model, statusP, x.colvalue)
+        # Don't `_check_ret(ret)` here, just bail is there isn't a dual ray.
+        x.has_primal_ray = (ret == kHighsStatusOk) && (statusP[] == 1)
+        if x.has_primal_ray
+            x.primal_solution_status = kHighsSolutionStatusNone
+        end
+    end
     return
 end
 

test/MOI_wrapper.jl

@@ -870,11 +870,8 @@ function test_infeasible_point()
     MOI.optimize!(model)
     @test MOI.get(model, MOI.TerminationStatus()) == MOI.INFEASIBLE
     @test MOI.get(model, MOI.PrimalStatus()) == MOI.INFEASIBLE_POINT
-    dual_stat = MOI.get(model, MOI.DualStatus())
-    @test dual_stat in (MOI.INFEASIBLE_POINT, MOI.NO_SOLUTION)
     @test MOI.get(model, MOI.ResultCount()) == 1
     @test MOI.get(model, MOI.VariablePrimal(), x) isa Vector{Float64}
-    @test MOI.get(model, MOI.ConstraintDual(), ci) isa Float64
     return
 end
 
@@ -1178,6 +1175,16 @@ function test_ConstraintName()
     return
 end
 
+function test_issue_287()
+    model = HiGHS.Optimizer()
+    x, _ = MOI.add_constrained_variable(model, MOI.ZeroOne())
+    MOI.add_constraint(model, 1.0 * x, MOI.EqualTo(0.5))
+    MOI.optimize!(model)
+    @test MOI.get(model, MOI.TerminationStatus()) == MOI.INFEASIBLE
+    @test MOI.get(model, MOI.PrimalStatus()) == MOI.NO_SOLUTION
+    return
+end
+
 end  # module
 
 TestMOIHighs.runtests()