Fix computing primal status with a VectorNonlinearOracle (#511)

Author: odow

ext/IpoptMathOptInterfaceExt/MOI_wrapper.jl

@@ -1466,14 +1466,19 @@ end
 
 function _manually_evaluated_primal_status(model::Optimizer)
     x, g = model.inner.x, model.inner.g
-    m, n = length(g), length(x)
     x_L, x_U = model.variables.lower, model.variables.upper
     g_L, g_U = copy(model.qp_data.g_L), copy(model.qp_data.g_U)
-    # Assuming constraints are guaranteed to be in the order [qp_cons, nlp_cons]
+    # Assuming constraints are guaranteed to be in the order:
+    # [qp_cons, nlp_cons, oracle]
     for bound in model.nlp_data.constraint_bounds
         push!(g_L, bound.lower)
         push!(g_U, bound.upper)
     end
+    for (_, cache) in model.vector_nonlinear_oracle_constraints
+        append!(g_L, cache.set.l)
+        append!(g_U, cache.set.u)
+    end
+    m, n = length(g_L), length(x)
     # 1e-8 is the default tolerance
     tol = get(model.options, "tol", 1e-8)
     if all(x_L[i] - tol <= x[i] <= x_U[i] + tol for i in 1:n) &&

test/MOI_wrapper.jl

@@ -936,6 +936,9 @@ function test_vector_nonlinear_oracle()
         @test 0.9 < cache.eval_jacobian_timer < 4
         @test 0.9 < cache.eval_hessian_lagrangian_timer < 2
     end
+    MOI.set(model, MOI.RawOptimizerAttribute("max_iter"), 0)
+    MOI.optimize!(model)
+    @test MOI.get(model, MOI.PrimalStatus()) == MOI.INFEASIBLE_POINT
     return
 end