Skip to content

Instantly share code, notes, and snippets.

@kersulis
Created July 19, 2018 15:43
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save kersulis/263bd12f30dcdfc158e589238e756854 to your computer and use it in GitHub Desktop.
Save kersulis/263bd12f30dcdfc158e589238e756854 to your computer and use it in GitHub Desktop.
Timing results for SDPDecompPowerModel
using PowerModels, CSV, Mosek, Ipopt
mosek_solver = MosekSolver()
ipopt_solver = IpoptSolver(tol=1e-6)
nesta_dir = "../nesta/opf"
cases = [
"nesta_case3_lmbd",
"nesta_case5_pjm",
"nesta_case14_ieee",
"nesta_case24_ieee_rts",
"nesta_case30_as",
"nesta_case30_fsr",
"nesta_case30_ieee",
"nesta_case39_epri",
"nesta_case57_ieee",
"nesta_case73_ieee_rts",
"nesta_case89_pegase",
"nesta_case118_ieee",
"nesta_case162_ieee_dtc",
"nesta_case300_ieee",
"nesta_case1354_pegase"
]
gaps = Float64[]
times_ac = Float64[]
times_sdp = Float64[]
for case in cases
fname = joinpath(nesta_dir, case) * ".m"
result_ac = run_opf(fname, ACPPowerModel, ipopt_solver)
result_sdp = run_opf(fname, SDPDecompPowerModel, mosek_solver)
push!(times_ac, result_ac["solve_time"])
push!(times_sdp, result_sdp["solve_time"])
push!(gaps, 100*(result_ac["objective"] - result_sdp["objective"])/result_ac["objective"])
end
case gap time_ac time_sdp
nesta_case3_lmbd 0.380212 0.0115185 0.00657511
nesta_case5_pjm 5.21946 0.0187869 0.026509
nesta_case14_ieee 0.000147327 0.0249393 0.0432889
nesta_case24_ieee_rts 4.66891e-5 0.0591923 0.166873
nesta_case30_as 0.000902287 0.0485134 0.110496
nesta_case30_fsr 2.73208e-5 0.0556282 0.234215
nesta_case30_ieee 2.83224e-5 0.0627696 0.127613
nesta_case39_epri 0.0151207 0.0661744 0.257961
nesta_case57_ieee 0.000367954 0.105879 0.340923
nesta_case73_ieee_rts 0.00105411 0.222324 0.584769
nesta_case89_pegase NaN 0.355663 0.00903893
nesta_case118_ieee 0.0644316 0.403797 1.06378
nesta_case162_ieee_dtc 1.07609 0.665934 5.35741
nesta_case300_ieee 0.111634 1.18858 4.20815
nesta_case1354_pegase NaN 5.7544 0.499721
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment