I’m having some issues with Ipopt not finding a solution with a script I’m updating to Julia 1.1 / JuMP 0.19 from a functioning one on Julia 0.64 / JuMP 0.18.
Unsure exactly where the problem is, so I’ve put the script in a gist so you can see the changes I’ve made.
Diff between the two versions
The changes here, as far as I can tell, are the only ones that need to be made. Arrays need to conform to the v1.1 syntax, and then the new JuMP syntax is also altered, although nothing else is done.
Is there something obvious that I’ve missed? Are there some pointers you could give me to help solve this problem?
I’m expecting the optimization to converge quite quickly (as you can see, the old version takes 5 seconds), but instead it just hits this -1.2028687e+12 objective value and goes on forever, with no changes to the output. I’ve waited for over an hour with no change. So somehow my solution has become infeasible. This to me suggests either an oversight in my implementation, or has identified some edge case in the new JuMP version.
No, the problem changes. The initial run sets the \psi_2 parameter to 0, then is altered based on the result. Strictly speaking this will ultimately be two runs rather than three. The initial call is more to do with another portion of the script (see below).
Yep, this was just to try and drop down the complexity of the example. This isn’t an issue in the real codebase.
Also, try simplifying the problem as much as possible so that it still demonstrates the problem. It could be for any number of reasons, although it is probably a problem handling pψ₂.
I’ve been cutting down the job to not run multiple times in an attempt to simplify the example for you. If I hard code the result of CPRICE from the first run and only optimize! one time, I get the correct solution in half the time as the v0.64/v0.18 version.
Thanks for your help so far. I’ll work on it a bit more—see if I can identify why the parameter isn’t doing what I expect, then will come back to you if I get stumped.