Add odow's bench

Author: blegat

perf/opf/Project.toml

@@ -3,4 +3,5 @@ BenchmarkTools = "6e4b80f9-dd63-53aa-95a3-0cdb28fa8baf"
 Ipopt = "b6b21f68-93f8-5de0-b562-5493be1d77c9"
 JuMP = "4076af6c-e467-56ae-b986-b466b2749572"
 MathOptInterface = "b8f27783-ece8-5eb3-8dc8-9495eed66fee"
+PGLib = "07a8691f-3d11-4330-951b-3c50f98338be"
 PowerModels = "c36e90e8-916a-50a6-bd94-075b64ef4655"

perf/opf/opf2.jl

@@ -0,0 +1,32 @@
+using Revise
+import Ipopt
+import JuMP
+import MathOptInterface as MOI
+import PGLib
+import PowerModels
+
+model = JuMP.direct_model(Ipopt.Optimizer())
+pm = PowerModels.instantiate_model(
+    PGLib.pglib("pglib_opf_case10000_goc"),
+    PowerModels.ACPPowerModel,
+    PowerModels.build_opf;
+    jump_model = model,
+);
+
+ipopt = JuMP.backend(model)
+x = MOI.get(ipopt, MOI.ListOfVariableIndices())
+m, n = length(ipopt.nlp_model.constraints), length(x)
+
+evaluator = MOI.Nonlinear.Evaluator(
+    ipopt.nlp_model,
+    MOI.Nonlinear.SparseReverseMode(),
+    x,
+)
+MOI.initialize(evaluator, [:Grad, :Jac, :Hess])
+
+H_struct = MOI.hessian_lagrangian_structure(evaluator)
+H = zeros(length(H_struct))
+mu = rand(m)
+sigma = 0.0
+x_v = rand(n)
+@time MOI.eval_hessian_lagrangian(evaluator, H, x_v, sigma, mu)