# Optimizer returns status infeasible for a feasible problem

**URL:** <https://discourse.julialang.org/t/optimizer-returns-status-infeasible-for-a-feasible-problem/55064>\
**Category:** Optimization (Mathematical)\
**Tags:** jump\
**Created:** [February 11, 2021, 2:52pm UTC](https://discourse.julialang.org/t/optimizer-returns-status-infeasible-for-a-feasible-problem/55064 "2021-02-11T14:52:40Z")\
**Posts on this page:** 3\
**Page:** 1

<div class="post-metadata">

**Author:** ![htpusa](https://avatars.discourse-cdn.com/v4/letter/h/d2c977/32.png) [@htpusa](https://discourse.julialang.org/u/htpusa)\
**Post date:** [February 11, 2021, 2:52pm UTC](https://discourse.julialang.org/t/optimizer-returns-status-infeasible-for-a-feasible-problem/55064/1 "2021-02-11T14:52:40Z")

</div>

I’m running into the following issue while solving a set of LP-problems using JuMP:

I want to create a model, optimize it, and then in turn minimize and maximize its variables while requiring the original optimum to be achieved. For those in the know, yes, it’s FVA. The code would look something like this

```julia
using JuMP, CPLEX
model = JuMP.Model(CPLEX.Optimizer)
S = rand(5,5)
lower = zeros(5)
upper = ones(5)
c = [1; ones(4)]'
b = zeros(5)
@variable(model, x[i=1:5], lower_bound=lower[i], upper_bound=upper[i])
@objective(model, Max, c * x)
@constraint(model, S * x .== b)
JuMP.optimize!(model)
optimum = JuMP.objective_value(model)
@constraint(model, c * x ≥ optimum)

for i in 1:5
      @objective(model, MOI.MIN_SENSE, x[i])
      JuMP.optimize!(model)
      min = JuMP.objective_value(model)

      JuMP.set_objective_sense(model, MOI.MAX_SENSE)
      JuMP.optimize!(model)
      max = JuMP.objective_value(model)
end

```

Now by definition, if the original problem is feasible, all of the subproblems within the for-loop should be as well. However, it appears that sometimes the optimizer enters into a “degenerate” state, and return infeasible for one of the subproblems. If a new optimizer instance is given, it returns feasible as should. In other words

```julia
JuMP.optimize!(model)
# termination_status(model) will return
# INFEASIBLE::TerminationStatusCode = 2

set_optimizer(model, CPLEX.Optimizer);
JuMP.optimize!(model)
# termination_status(model) will return
# OPTIMAL::TerminationStatusCode = 1

```

I’m at a loss as to what is going on, I’m not even sure at what level the problem is. I originally thought the problem was related to CPLEX but it appears it happens with Gurobi as well (but not necessarily with the exact same model).

I’m on  
Julia Version 1.3.0  
Commit 46ce4d7933 (2019-11-26 06:09 UTC)  
Platform Info:  
OS: Linux (x86\_64-pc-linux-gnu)

with  
JuMP v0.21.3  
CPLEX v0.6.6

---

<div class="post-metadata">

**Author:** ![rocarvaj](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/rocarvaj/32/4769_2.png) [@rocarvaj](https://discourse.julialang.org/u/rocarvaj)\
**Post date:** [February 11, 2021, 7:17pm UTC](https://discourse.julialang.org/t/optimizer-returns-status-infeasible-for-a-feasible-problem/55064/2 "2021-02-11T19:17:28Z")

</div>

Is that a typo in `termination_status(optimization_model)`? It should say `model`, right? I just want to make sure.

---

<div class="post-metadata">

**Author:** ![odow](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/odow/32/28685_2.png) [@odow](https://discourse.julialang.org/u/odow)\
**Post date:** [February 11, 2021, 7:27pm UTC](https://discourse.julialang.org/t/optimizer-returns-status-infeasible-for-a-feasible-problem/55064/3 "2021-02-11T19:27:33Z")

</div>

Please read the first post of [Please read: make it easier to help you](https://discourse.julialang.org/t/psa-make-it-easier-to-help-you/14757). You should format your code, and provide a minimal working example. I can’t run your code because I don’t know what `lower`, `upper`, `c`, `b`, and `S` are.

I’ve seen infeasibilities like this occur because the presolve tolerance of the solver is different to the final tolerance. (So using an optimal solution in the next solve may be classified as infeasible.)

I suggest you try `@constraint(model, c * x ≥ optimum - 1e-4)` (or some other similar tolerance).