AC Optimal Power Flow in Various Nonlinear Optimization Frameworks

tmigot · February 16, 2024, 9:51pm

Congrat’s @sshin23 that’s impressive new option!

I agree with the idea of renaming NPLModels as ADNLPModels, especially since ExaModels is also implementing the NLPModel API.

amontoison · February 17, 2024, 6:11am

I suggest to also specify which backend of BLAS / LAPACK is used for the benchmarks.

Is it OpenBLAS32_jll.jl?

Ipopt_jll.jl and HSL_jll.jl (from libHSL) are compiled with LBT and it could highly impact the performances is you use MKL.jl or AppleAccelerate.jl for example.

amontoison · February 18, 2024, 4:55am

I checked the Manifest.toml of rosetta-opf and an old version of Ipopt (artifact Ipopt_jll.jl) is used:

github.com

lanl-ansi/rosetta-opf/blob/main/Manifest.toml#L667-L671


      
          [[deps.Ipopt_jll]]
          deps = ["ASL_jll", "Artifacts", "CompilerSupportLibraries_jll", "JLLWrappers", "Libdl", "MUMPS_seq_jll", "OpenBLAS32_jll", "Pkg"]
          git-tree-sha1 = "e3e202237d93f18856b6ff1016166b0f172a49a8"
          uuid = "9cc047cb-c261-5740-88fc-0cf96f7bdcc7"
          version = "300.1400.400+0"

The lastest version is 300.1400.1400.
The version 300.1400.400 is the last version supported by Julia 1.6 and it was compiled 2021.

ccoffrin · February 18, 2024, 4:50pm

@tmigot, I have updated NPLModels to ADNLPModels. This will be reflected in the next update to the table.

@amontoison, I forget exactly which dependency combination is holding the Ipopt_jll version down. @odow had looked into it once and may remember more details. For the linear solver, Ipopt_jll is pulling an HSL DLL that is installed separately from Julia, as per these instructions. I am interested to update to the latest best practices for using HSL in Julia, but probably need some help in getting it setup.

odow · February 18, 2024, 8:13pm

Ah. It’s because of OptimizationMOI. Let me take a look. They shouldn’t need that compat

github.com

SciML/Optimization.jl/blob/4246bb0c58b3b5147613247ef6895251f3395029/lib/OptimizationMOI/Project.toml#L17


      
          Ipopt_jll = "9cc047cb-c261-5740-88fc-0cf96f7bdcc7"
          MathOptInterface = "b8f27783-ece8-5eb3-8dc8-9495eed66fee"
          ModelingToolkit = "961ee093-0014-501f-94e3-6117800e7a78"
          Optimization = "7f7a1694-90dd-40f0-9382-eb1efda571ba"
          Reexport = "189a3867-3050-52da-a836-e630ba90ab69"
          SparseArrays = "2f01184e-e22b-5df5-ae63-d93ebab69eaf"
          SymbolicIndexingInterface = "2efcf032-c050-4f8e-a9bb-153293bab1f5"
          Symbolics = "0c5d862f-8b57-4792-8d23-62f2024744c7"
          
          [compat]
          Ipopt_jll = "=300.1400.400"
          Juniper = "0.9"
          MathOptInterface = "1"
          ModelingToolkit = "8.74"
          Optimization = "3.21"
          Reexport = "1.2"
          SymbolicIndexingInterface = "0.3"
          Symbolics = "5"
          julia = "1"
          
          [extras]

See [OptimizationMOI] fix compat bounds in Project.toml by odow · Pull Request #705 · SciML/Optimization.jl · GitHub

odow · February 18, 2024, 8:21pm

@ccoffrin, you can change the BLAS backend in addition to HSL.

See GitHub - jump-dev/Ipopt.jl: Julia interface to the Ipopt nonlinear solver

This requires the updated Ipopt_jll, but it should “just” be a matter of adding using MKL before using Ipopt.

amontoison · February 20, 2024, 5:03am

@ccoffrin
For the HSL linear solvers, you just need to download and install HSL_jll.jl now (require an academic licence).
Please take the version precompiled with LBT.
With a using HSL_jll, you can direclty use the HSL linear solvers in Ipopt:

odow · February 20, 2024, 5:21am

Actually, we probably can’t update anyway, because of:

ccoffrin · February 20, 2024, 5:29am

That only impacts the variants Project.toml, not the main one right?

odow · February 20, 2024, 11:19pm

Now fixed in AmplNLWriter v1.2.1 and OptimizationMOI v0.4.0. I’ll make a PR