# Unable to solve the problem using Lagrangian duality handler

**URL:** https://discourse.julialang.org/t/unable-to-solve-the-problem-using-lagrangian-duality-handler/129079
**Category:** Optimization (Mathematical)
**Tags:** sddp
**Created:** [May 16, 2025, 8:15pm UTC](https://discourse.julialang.org/t/unable-to-solve-the-problem-using-lagrangian-duality-handler/129079 "2025-05-16T20:15:51Z")
**Posts on this page:** 3
**Page:** 1

<div class="post-metadata">

### Author: ![Engr\_Moiz\_Ahmad](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/engr_moiz_ahmad/32/214621_2.png) [@Engr\_Moiz\_Ahmad](https://discourse.julialang.org/u/Engr_Moiz_Ahmad)
#### Post date: [May 16, 2025, 8:15pm UTC](https://discourse.julialang.org/t/unable-to-solve-the-problem-using-lagrangian-duality-handler/129079/1 "2025-05-16T20:15:52Z")

</div>

I am dealing with a multi-stage mixed-integer stochastic global supply chain optimization problem and want to solve it using the SDDiP algorithm in the [SDDP.jl](https://juliaregistries.github.io/General/packages/redirect_to_repo/SDDP) package in Julia.

In this regard, I tried solving the mixed-integer linear program using the Lagrangian duality handler. I am getting the following error as a result. Is it the issue with the nature of the formulation and the scale of the problem instance? The SDDiP algorithm (with CPLEX solver) was executing with other duality handlers.

 ![image](https://global.discourse-cdn.com/julialang/original/3X/1/2/12aca3eecb5ccada112d46455a6075381f5ab1ec.png)

---

<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: [May 16, 2025, 8:42pm UTC](https://discourse.julialang.org/t/unable-to-solve-the-problem-using-lagrangian-duality-handler/129079/2 "2025-05-16T20:42:16Z")

</div>

3019 means that there were unrecoverable numerical issues. It’s hard to say why without a reproducible example. Does it work with Gurobi instead?

I’ve previously found that CPLEX can be faster in some instances but it is more likely to fail with errors like this because it is less numerically stable.

---

<div class="post-metadata">

### Author: ![Engr\_Moiz\_Ahmad](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/engr_moiz_ahmad/32/214621_2.png) [@Engr\_Moiz\_Ahmad](https://discourse.julialang.org/u/Engr_Moiz_Ahmad)
#### Post date: [May 17, 2025, 9:40pm UTC](https://discourse.julialang.org/t/unable-to-solve-the-problem-using-lagrangian-duality-handler/129079/3 "2025-05-17T21:40:17Z")

</div>

Got it, Gurobi works, but it is taking a significantly longer time to execute under Lagrangian duality handler.
