From e4b1235d3c463addc6ee9bb9829d44bd6c41b3c1 Mon Sep 17 00:00:00 2001 From: Carleton Coffrin Date: Sat, 5 Mar 2022 17:32:59 -0700 Subject: [PATCH] rename test solver variables (#814) --- test/data-basic.jl | 12 +- test/data.jl | 50 ++++---- test/matpower.jl | 16 +-- test/model.jl | 6 +- test/modify.jl | 10 +- test/multinetwork.jl | 42 +++---- test/opb.jl | 38 +++--- test/opf-obj.jl | 24 ++-- test/opf-ptdf.jl | 24 ++-- test/opf-var.jl | 152 ++++++++++++------------ test/opf.jl | 272 +++++++++++++++++++++---------------------- test/ots.jl | 62 +++++----- test/output.jl | 38 +++--- test/pf-native.jl | 34 +++--- test/pf.jl | 74 ++++++------ test/psse.jl | 66 +++++------ test/pti.jl | 4 +- test/runtests.jl | 10 +- test/tnep.jl | 28 ++--- test/util.jl | 52 ++++----- test/warmstart.jl | 24 ++-- 21 files changed, 519 insertions(+), 519 deletions(-) diff --git a/test/data-basic.jl b/test/data-basic.jl index a804b2edb..48cc86180 100644 --- a/test/data-basic.jl +++ b/test/data-basic.jl @@ -15,7 +15,7 @@ @test length(data["switch"]) == 0 @test length(data["storage"]) == 0 - result = run_opf(data, ACPPowerModel, ipopt_solver) + result = run_opf(data, ACPPowerModel, nlp_solver) @test isapprox(result["objective"], 1036.52; atol=1e0) end @@ -31,7 +31,7 @@ @test length(data["switch"]) == 0 @test length(data["storage"]) == 0 - result = run_opf(data, ACPPowerModel, ipopt_solver) + result = run_opf(data, ACPPowerModel, nlp_solver) @test isapprox(result["objective"], 16551.7; atol=1e0) end @@ -43,7 +43,7 @@ end @testset "basic bus injection" begin data = make_basic_network(PowerModels.parse_file("../test/data/matpower/case14.m")) - result = run_opf(data, DCPPowerModel, ipopt_solver) + result = run_opf(data, DCPPowerModel, nlp_solver) update_data!(data, result["solution"]) bi = calc_basic_bus_injection(data) @@ -51,7 +51,7 @@ end @test isapprox(real(sum(bi)), 0.0; atol=1e-6) - result = run_opf(data, ACPPowerModel, ipopt_solver) + result = run_opf(data, ACPPowerModel, nlp_solver) update_data!(data, result["solution"]) bi = calc_basic_bus_injection(data) @@ -95,7 +95,7 @@ end @test isapprox(SM_1, SM_2; atol=1e-6) - result = run_opf(data, DCPPowerModel, ipopt_solver) + result = run_opf(data, DCPPowerModel, nlp_solver) update_data!(data, result["solution"]) va = angle.(calc_basic_bus_voltage(data)) @@ -155,7 +155,7 @@ end @test isapprox(sum(P), 0.9894736; atol=1e-6) - result = run_opf(data, DCPPowerModel, ipopt_solver) + result = run_opf(data, DCPPowerModel, nlp_solver) update_data!(data, result["solution"]) bi = real(calc_basic_bus_injection(data)) diff --git a/test/data.jl b/test/data.jl index eb64cc10b..e32ef8ef1 100644 --- a/test/data.jl +++ b/test/data.jl @@ -30,7 +30,7 @@ TESTLOG = Memento.getlogger(PowerModels) end @testset "5-bus solution summary from dict" begin - result = run_ac_opf("../test/data/matpower/case5.m", ipopt_solver) + result = run_ac_opf("../test/data/matpower/case5.m", nlp_solver) output = sprint(PowerModels.summary, result["solution"]) line_count = count(c -> c == '\n', output) @@ -128,7 +128,7 @@ end @testset "3-bus case solution" begin - result = run_ac_opf("../test/data/matpower/case3.m", ipopt_solver) + result = run_ac_opf("../test/data/matpower/case3.m", nlp_solver) result_base = deepcopy(result) PowerModels.make_mixed_units!(result["solution"]) @@ -137,7 +137,7 @@ end @test InfrastructureModels.compare_dict(result, result_base) end @testset "5-bus case solution" begin - result = run_ac_opf("../test/data/matpower/case5_asym.m", ipopt_solver) + result = run_ac_opf("../test/data/matpower/case5_asym.m", nlp_solver) result_base = deepcopy(result) PowerModels.make_mixed_units!(result["solution"]) @@ -146,7 +146,7 @@ end @test InfrastructureModels.compare_dict(result, result_base) end @testset "24-bus case solution" begin - result = run_ac_opf("../test/data/matpower/case24.m", ipopt_solver) + result = run_ac_opf("../test/data/matpower/case24.m", nlp_solver) result_base = deepcopy(result) PowerModels.make_mixed_units!(result["solution"]) @@ -157,7 +157,7 @@ end @testset "5-bus case solution with duals" begin - result = run_dc_opf("../test/data/matpower/case5.m", ipopt_solver, setting = Dict("output" => Dict("branch_flows" => true, "duals" => true))) + result = run_dc_opf("../test/data/matpower/case5.m", nlp_solver, setting = Dict("output" => Dict("branch_flows" => true, "duals" => true))) result_base = deepcopy(result) PowerModels.make_mixed_units!(result["solution"]) @@ -437,7 +437,7 @@ end @testset "output values" begin data = PowerModels.parse_file("../test/data/matpower/case7_tplgy.m") PowerModels.simplify_network!(data) - result = run_opf(data, ACPPowerModel, ipopt_solver) + result = run_opf(data, ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 1778; atol = 1e0) @@ -594,7 +594,7 @@ end data = PowerModels.parse_file("../test/data/matpower/case5.m") data["branch"]["4"]["br_status"] = 0 data["buspairs"] = PowerModels.calc_buspair_parameters(data["bus"], data["branch"], 1:1, false) - result = run_opf(data, ACPPowerModel, ipopt_solver) + result = run_opf(data, ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 16642; atol = 1e0) @@ -608,7 +608,7 @@ end @testset "5-bus ac polar flow" begin data = PowerModels.parse_file("../test/data/matpower/case5.m") data["branch"]["4"]["br_status"] = 0 - result = run_opf(data, ACPPowerModel, ipopt_solver) + result = run_opf(data, ACPPowerModel, nlp_solver) PowerModels.update_data!(data, result["solution"]) ac_flows = PowerModels.calc_branch_flow_ac(data) @@ -626,7 +626,7 @@ end @testset "5-bus ac rect flow" begin data = PowerModels.parse_file("../test/data/matpower/case5.m") data["branch"]["4"]["br_status"] = 0 - result = run_opf(data, ACRPowerModel, ipopt_solver, solution_processors=[sol_data_model!]) + result = run_opf(data, ACRPowerModel, nlp_solver, solution_processors=[sol_data_model!]) PowerModels.update_data!(data, result["solution"]) ac_flows = PowerModels.calc_branch_flow_ac(data) @@ -644,7 +644,7 @@ end @testset "5-bus dc flow" begin data = PowerModels.parse_file("../test/data/matpower/case5.m") data["branch"]["4"]["br_status"] = 0 - result = run_opf(data, DCPPowerModel, ipopt_solver) + result = run_opf(data, DCPPowerModel, nlp_solver) PowerModels.update_data!(data, result["solution"]) dc_flows = PowerModels.calc_branch_flow_dc(data) @@ -666,7 +666,7 @@ end @testset "5-bus polynomial gen cost" begin data = PowerModels.parse_file("../test/data/matpower/case5.m") - result = run_opf(data, ACPPowerModel, ipopt_solver) + result = run_opf(data, ACPPowerModel, nlp_solver) PowerModels.update_data!(data, result["solution"]) gen_cost = PowerModels.calc_gen_cost(data) @@ -681,7 +681,7 @@ end @test isa(gen["ncost"], Int) end - result = run_opf(data, ACPPowerModel, ipopt_solver) + result = run_opf(data, ACPPowerModel, nlp_solver) PowerModels.update_data!(data, result["solution"]) gen_cost = PowerModels.calc_gen_cost(data) @@ -700,7 +700,7 @@ end @testset "5-bus polynomial gen and dcline cost" begin data = PowerModels.parse_file("../test/data/matpower/case5_dc.m") - result = run_opf(data, ACPPowerModel, ipopt_solver) + result = run_opf(data, ACPPowerModel, nlp_solver) PowerModels.update_data!(data, result["solution"]) gen_cost = PowerModels.calc_gen_cost(data) @@ -713,7 +713,7 @@ end data["gen"]["1"]["gen_status"] = 0 data["dcline"]["1"]["br_status"] = 0 - result = run_opf(data, ACPPowerModel, ipopt_solver) + result = run_opf(data, ACPPowerModel, nlp_solver) PowerModels.update_data!(data, result["solution"]) gen_cost = PowerModels.calc_gen_cost(data) @@ -729,7 +729,7 @@ end @testset "5-bus ac polar balance" begin data = PowerModels.parse_file("../test/data/matpower/case5_dc.m") data["branch"]["4"]["br_status"] = 0 - result = run_opf(data, ACPPowerModel, ipopt_solver) + result = run_opf(data, ACPPowerModel, nlp_solver) PowerModels.update_data!(data, result["solution"]) balance = PowerModels.calc_power_balance(data) @@ -743,7 +743,7 @@ end @testset "5-bus dc balance" begin data = PowerModels.parse_file("../test/data/matpower/case5_dc.m") data["branch"]["4"]["br_status"] = 0 - result = run_opf(data, DCPPowerModel, ipopt_solver) + result = run_opf(data, DCPPowerModel, nlp_solver) PowerModels.update_data!(data, result["solution"]) balance = PowerModels.calc_power_balance(data) @@ -758,7 +758,7 @@ end @testset "5-bus ac polar balance with storage" begin data = PowerModels.parse_file("../test/data/matpower/case5_strg.m") data["branch"]["4"]["br_status"] = 0 - result = PowerModels._solve_opf_strg(data, ACPPowerModel, ipopt_solver) + result = PowerModels._solve_opf_strg(data, ACPPowerModel, nlp_solver) PowerModels.update_data!(data, result["solution"]) balance = PowerModels.calc_power_balance(data) @@ -772,7 +772,7 @@ end @testset "5-bus dc balance with storage" begin data = PowerModels.parse_file("../test/data/matpower/case5_strg.m") data["branch"]["4"]["br_status"] = 0 - result = PowerModels._solve_opf_strg(data, DCPPowerModel, ipopt_solver) + result = PowerModels._solve_opf_strg(data, DCPPowerModel, nlp_solver) PowerModels.update_data!(data, result["solution"]) balance = PowerModels.calc_power_balance(data) @@ -787,7 +787,7 @@ end @testset "5-bus balance from flow ac" begin data = PowerModels.parse_file("../test/data/matpower/case5_dc.m") data["branch"]["4"]["br_status"] = 0 - result = run_opf(data, ACPPowerModel, ipopt_solver) + result = run_opf(data, ACPPowerModel, nlp_solver) PowerModels.update_data!(data, result["solution"]) flows = PowerModels.calc_branch_flow_ac(data) @@ -804,7 +804,7 @@ end @testset "5-bus balance from flow dc" begin data = PowerModels.parse_file("../test/data/matpower/case5_dc.m") data["branch"]["4"]["br_status"] = 0 - result = run_opf(data, DCPPowerModel, ipopt_solver) + result = run_opf(data, DCPPowerModel, nlp_solver) PowerModels.update_data!(data, result["solution"]) flows = PowerModels.calc_branch_flow_dc(data) @@ -825,22 +825,22 @@ end @testset "test renumber bus ids" begin @testset "5-bus with dcline" begin data = PowerModels.parse_file("../test/data/matpower/case5_dc.m") - result_1 = run_opf(data, ACPPowerModel, ipopt_solver) + result_1 = run_opf(data, ACPPowerModel, nlp_solver) update_bus_ids!(data, Dict(1 => 10, 2=>20, 3=>30, 4=>40, 5=>50)) - result_2 = run_opf(data, ACPPowerModel, ipopt_solver) + result_2 = run_opf(data, ACPPowerModel, nlp_solver) @test isapprox(result_1["objective"], result_2["objective"]; atol=1e-6) end @testset "5-bus with switches" begin data = PowerModels.parse_file("../test/data/matpower/case5_sw.m") - result_1 = PowerModels._solve_opf_sw(data, DCPPowerModel, ipopt_solver) + result_1 = PowerModels._solve_opf_sw(data, DCPPowerModel, nlp_solver) update_bus_ids!(data, Dict(1 => 10, 2=>20, 3=>30, 4=>40, 10=>100)) - result_2 = PowerModels._solve_opf_sw(data, DCPPowerModel, ipopt_solver) + result_2 = PowerModels._solve_opf_sw(data, DCPPowerModel, nlp_solver) @test isapprox(result_1["objective"], result_2["objective"]; atol=1e-6) end @@ -854,7 +854,7 @@ end @test length(data["switch"]) == 0 @test length(data["bus"]) == 4 - result = run_opf(data, ACPPowerModel, ipopt_solver) + result = run_opf(data, ACPPowerModel, nlp_solver) @test isapprox(result["objective"], 16641.20; atol=1e0) end diff --git a/test/matpower.jl b/test/matpower.jl index f052e527e..7ebaffc71 100644 --- a/test/matpower.jl +++ b/test/matpower.jl @@ -1,6 +1,6 @@ @testset "test matpower parser" begin @testset "30-bus case file" begin - result = run_opf("../test/data/matpower/case30.m", ACPPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case30.m", ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 204.96; atol = 1e-1) @@ -10,7 +10,7 @@ data = PowerModels.parse_file("../test/data/matpower/case30.m") @test isa(JSON.json(data), String) - result = run_opf(data, ACPPowerModel, ipopt_solver) + result = run_opf(data, ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 204.96; atol = 1e-1) @@ -20,7 +20,7 @@ data = PowerModels.parse_matpower("../test/data/matpower/case30.m") @test isa(JSON.json(data), String) - result = run_opf(data, ACPPowerModel, ipopt_solver) + result = run_opf(data, ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 204.96; atol = 1e-1) @@ -31,7 +31,7 @@ data = PowerModels.parse_matpower(f) @test isa(JSON.json(data), String) - result = run_opf(data, ACPPowerModel, ipopt_solver) + result = run_opf(data, ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 204.96; atol = 1e-1) @@ -57,7 +57,7 @@ end @testset "2-bus case file with spaces" begin - result = run_pf("../test/data/matpower/case2.m", ACPPowerModel, ipopt_solver) + result = run_pf("../test/data/matpower/case2.m", ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 0.0; atol = 1e-1) @@ -76,20 +76,20 @@ end @testset "test matpower data coercion" begin @testset "ACP Model" begin - result = run_opf("../test/data/matpower/case14.m", ACPPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case14.m", ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 8081.5; atol = 1e0) #@test result["status"] = bus_name end @testset "DC Model" begin - result = run_opf("../test/data/matpower/case14.m", DCPPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case14.m", DCPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 7642.6; atol = 1e0) end @testset "QC Model" begin - result = run_opf("../test/data/matpower/case14.m", QCRMPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case14.m", QCRMPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 8075.1; atol = 1e0) diff --git a/test/model.jl b/test/model.jl index bddd3e5ee..52b5eb330 100644 --- a/test/model.jl +++ b/test/model.jl @@ -3,7 +3,7 @@ @testset "run with user provided JuMP model" begin m = JuMP.Model() x = JuMP.@variable(m, my_var >= 0, start=0.0) - result = run_ac_opf("../test/data/matpower/case5.m", ipopt_solver, jump_model=m) + result = run_ac_opf("../test/data/matpower/case5.m", nlp_solver, jump_model=m) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 18269; atol = 1e0) @@ -49,7 +49,7 @@ end @testset "relax integrality" begin @testset "relax OTS model" begin - result = run_ots("../test/data/matpower/case5.m", DCPPowerModel, ipopt_solver, relax_integrality=true) + result = run_ots("../test/data/matpower/case5.m", DCPPowerModel, nlp_solver, relax_integrality=true) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 14810.0; atol = 1e0) @@ -59,7 +59,7 @@ end end @testset "relax TNEP model" begin - result = run_tnep("../test/data/matpower/case5_tnep.m", SOCWRPowerModel, ipopt_solver, relax_integrality=true) + result = run_tnep("../test/data/matpower/case5_tnep.m", SOCWRPowerModel, nlp_solver, relax_integrality=true) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 0.1236; atol = 1e-2) diff --git a/test/modify.jl b/test/modify.jl index e08dbadd7..8554514e7 100644 --- a/test/modify.jl +++ b/test/modify.jl @@ -4,26 +4,26 @@ @testset "30-bus case file incremental" begin data = PowerModels.parse_file("../test/data/matpower/case30.m") - result = run_opf(data, ACPPowerModel, ipopt_solver) + result = run_opf(data, ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 204.96; atol = 1e-1) data["branch"]["6"]["br_status"] = 0 - result = run_opf(data, ACPPowerModel, ipopt_solver) + result = run_opf(data, ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 195.25; atol = 1e-1) data["gen"]["6"]["gen_status"] = 0 - result = run_opf(data, ACPPowerModel, ipopt_solver) + result = run_opf(data, ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 195.896; atol = 1e-1) data["load"]["4"]["pd"] = 0 data["load"]["4"]["qd"] = 0 - result = run_opf(data, ACPPowerModel, ipopt_solver) + result = run_opf(data, ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 104.428; atol = 1e-1) end @@ -53,7 +53,7 @@ PowerModels.update_data!(data, data_delta) - result = run_opf(data, ACPPowerModel, ipopt_solver) + result = run_opf(data, ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 104.428; atol = 1e-1) end diff --git a/test/multinetwork.jl b/test/multinetwork.jl index 0a6472833..2b8badc3c 100644 --- a/test/multinetwork.jl +++ b/test/multinetwork.jl @@ -111,17 +111,17 @@ TESTLOG = Memento.getlogger(PowerModels) @testset "test run_opf with multinetwork data" begin mn_data = build_mn_data("../test/data/matpower/case5.m") - @test_throws(TESTLOG, ErrorException, PowerModels.run_opf(mn_data, ACPPowerModel, ipopt_solver)) + @test_throws(TESTLOG, ErrorException, PowerModels.run_opf(mn_data, ACPPowerModel, nlp_solver)) end @testset "test run_mn_opf with single-network data" begin - @test_throws(TESTLOG, ErrorException, PowerModels.run_mn_opf("../test/data/matpower/case5.m", ACPPowerModel, ipopt_solver)) + @test_throws(TESTLOG, ErrorException, PowerModels.run_mn_opf("../test/data/matpower/case5.m", ACPPowerModel, nlp_solver)) end @testset "test multi-network solution" begin # test case where generator status is 1 but the gen_bus status is 0 mn_data = build_mn_data("../test/data/matpower/case5.m") - result = PowerModels.run_mn_opf(mn_data, ACPPowerModel, ipopt_solver) + result = PowerModels.run_mn_opf(mn_data, ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 36538.2; atol = 1e0) @@ -135,7 +135,7 @@ TESTLOG = Memento.getlogger(PowerModels) @testset "test dc polar opb" begin - result = PowerModels._solve_mn_opb(mn_data, DCPPowerModel, ipopt_solver) + result = PowerModels._solve_mn_opb(mn_data, DCPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 29620.0; atol = 1e0) @@ -153,7 +153,7 @@ TESTLOG = Memento.getlogger(PowerModels) @testset "test ac polar opf" begin - result = PowerModels.run_mn_opf(mn_data, ACPPowerModel, ipopt_solver) + result = PowerModels.run_mn_opf(mn_data, ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 35103.8; atol = 1e0) @@ -170,7 +170,7 @@ TESTLOG = Memento.getlogger(PowerModels) end @testset "test dc polar opf" begin - result = PowerModels.run_mn_opf(mn_data, DCPPowerModel, ipopt_solver) + result = PowerModels.run_mn_opf(mn_data, DCPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 34959.8; atol = 1e0) @@ -187,7 +187,7 @@ TESTLOG = Memento.getlogger(PowerModels) end @testset "test soc opf" begin - result = PowerModels.run_mn_opf(mn_data, SOCWRPowerModel, ipopt_solver) + result = PowerModels.run_mn_opf(mn_data, SOCWRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 29999.4; atol = 1e0) @@ -204,7 +204,7 @@ TESTLOG = Memento.getlogger(PowerModels) end @testset "test sdp with constraint decomposition opf" begin - result = PowerModels.run_mn_opf(mn_data, SparseSDPWRMPowerModel, scs_solver) + result = PowerModels.run_mn_opf(mn_data, SparseSDPWRMPowerModel, sdp_solver) # hits iteration limit in SCS v0.9, sense ALMOST_OPTIMAL @test result["termination_status"] == ALMOST_OPTIMAL @@ -223,7 +223,7 @@ TESTLOG = Memento.getlogger(PowerModels) end @testset "test nfa opf" begin - result = PowerModels.run_mn_opf(mn_data, NFAPowerModel, ipopt_solver) + result = PowerModels.run_mn_opf(mn_data, NFAPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 29620.0; atol = 1e0) @@ -244,7 +244,7 @@ TESTLOG = Memento.getlogger(PowerModels) mn_data = build_mn_data("../test/data/matpower/case5.m") @testset "test dc polar opf" begin - result = PowerModels.run_mn_opf(mn_data, DCPPowerModel, ipopt_solver, setting = Dict("output" => Dict("duals" => true))) + result = PowerModels.run_mn_opf(mn_data, DCPPowerModel, nlp_solver, setting = Dict("output" => Dict("duals" => true))) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 35226.4; atol = 1e0) @@ -271,21 +271,21 @@ TESTLOG = Memento.getlogger(PowerModels) mn_data = PowerModels.parse_files("../test/data/matpower/case14.m", "../test/data/matpower/case24.m") @testset "test ac polar opf" begin - result = PowerModels.run_mn_opf(mn_data, ACPPowerModel, ipopt_solver) + result = PowerModels.run_mn_opf(mn_data, ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 87886.5; atol = 1e0) end @testset "test ac rectangular opf" begin - result = PowerModels.run_mn_opf(mn_data, ACRPowerModel, ipopt_solver) + result = PowerModels.run_mn_opf(mn_data, ACRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 87886.5; atol = 1e0) end @testset "test soc opf" begin - result = PowerModels.run_mn_opf(mn_data, SOCWRPowerModel, ipopt_solver) + result = PowerModels.run_mn_opf(mn_data, SOCWRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 78765.8; atol = 1e0) @@ -297,7 +297,7 @@ TESTLOG = Memento.getlogger(PowerModels) mn_data = build_mn_data("../test/data/matpower/case5_strg.m", replicates=4) @testset "test ac polar opf" begin - result = PowerModels.run_mn_opf_strg(mn_data, PowerModels.ACPPowerModel, juniper_solver) + result = PowerModels.run_mn_opf_strg(mn_data, PowerModels.ACPPowerModel, minlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 70435.5; atol = 1e0) @@ -326,7 +326,7 @@ TESTLOG = Memento.getlogger(PowerModels) end @testset "test soc opf" begin - result = PowerModels.run_mn_opf_strg(mn_data, PowerModels.SOCWRPowerModel, juniper_solver) + result = PowerModels.run_mn_opf_strg(mn_data, PowerModels.SOCWRPowerModel, minlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 58853.5; atol = 1e0) @@ -353,7 +353,7 @@ TESTLOG = Memento.getlogger(PowerModels) end @testset "test soc bf opf" begin - result = PowerModels.run_mn_opf_bf_strg(mn_data, PowerModels.SOCBFPowerModel, juniper_solver) + result = PowerModels.run_mn_opf_bf_strg(mn_data, PowerModels.SOCBFPowerModel, minlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 58826.36; atol = 1e0) @@ -380,7 +380,7 @@ TESTLOG = Memento.getlogger(PowerModels) end @testset "test linear bf opf" begin - result = PowerModels.run_mn_opf_bf_strg(mn_data, PowerModels.BFAPowerModel, juniper_solver) + result = PowerModels.run_mn_opf_bf_strg(mn_data, PowerModels.BFAPowerModel, minlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 57980.0; atol = 1e0) @@ -421,7 +421,7 @@ TESTLOG = Memento.getlogger(PowerModels) end end - result = PowerModels.run_mn_opf_strg(mn_data, PowerModels.DCPPowerModel, juniper_solver) + result = PowerModels.run_mn_opf_strg(mn_data, PowerModels.DCPPowerModel, minlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 69703.10; atol = 1e0) @@ -438,7 +438,7 @@ TESTLOG = Memento.getlogger(PowerModels) mn_data["nw"]["1"]["storage"]["1"]["status"] = 0 # verify that storage activation does not cause error Memento.setlevel!(TESTLOG, "warn") - @test_warn(TESTLOG, "network data should specify time_elapsed, using 1.0 as a default", PowerModels.run_mn_opf_strg(mn_data, PowerModels.ACPPowerModel, juniper_solver)) + @test_warn(TESTLOG, "network data should specify time_elapsed, using 1.0 as a default", PowerModels.run_mn_opf_strg(mn_data, PowerModels.ACPPowerModel, minlp_solver)) Memento.setlevel!(TESTLOG, "error") end end @@ -447,13 +447,13 @@ TESTLOG = Memento.getlogger(PowerModels) @testset "test solution feedback" begin mn_data = build_mn_data("../test/data/matpower/case5_asym.m") - opf_result = PowerModels.run_mn_opf(mn_data, ACPPowerModel, ipopt_solver) + opf_result = PowerModels.run_mn_opf(mn_data, ACPPowerModel, nlp_solver) @test opf_result["termination_status"] == LOCALLY_SOLVED @test isapprox(opf_result["objective"], 35103.8; atol = 1e0) PowerModels.update_data!(mn_data, opf_result["solution"]) - pf_result = PowerModels._solve_mn_pf(mn_data, ACPPowerModel, ipopt_solver) + pf_result = PowerModels._solve_mn_pf(mn_data, ACPPowerModel, nlp_solver) @test pf_result["termination_status"] == LOCALLY_SOLVED @test isapprox(pf_result["objective"], 0.0; atol = 1e-3) diff --git a/test/opb.jl b/test/opb.jl index 124ecacc1..2a4f34879 100644 --- a/test/opb.jl +++ b/test/opb.jl @@ -2,67 +2,67 @@ @testset "test nfa opb" begin @testset "3-bus case" begin - result = run_opb("../test/data/matpower/case3.m", NFAPowerModel, ipopt_solver) + result = run_opb("../test/data/matpower/case3.m", NFAPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 5638.97; atol = 1e0) end @testset "5-bus tranformer swap case" begin - result = run_opb("../test/data/matpower/case5.m", NFAPowerModel, ipopt_solver) + result = run_opb("../test/data/matpower/case5.m", NFAPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 14810.0; atol = 1e0) end @testset "5-bus asymmetric case" begin - result = run_opb("../test/data/matpower/case5_asym.m", NFAPowerModel, ipopt_solver) + result = run_opb("../test/data/matpower/case5_asym.m", NFAPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 14810.0; atol = 1e0) end @testset "5-bus gap case" begin - result = run_opb("../test/data/matpower/case5_gap.m", NFAPowerModel, ipopt_solver) + result = run_opb("../test/data/matpower/case5_gap.m", NFAPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], -27410.0; atol = 1e0) end @testset "5-bus with dcline costs" begin - result = run_opb("../test/data/matpower/case5_dc.m", NFAPowerModel, ipopt_solver) + result = run_opb("../test/data/matpower/case5_dc.m", NFAPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 14810.0; atol = 1e0) end @testset "5-bus with asymmetric line charge" begin - result = run_opb("../test/data/pti/case5_alc.raw", NFAPowerModel, ipopt_solver) + result = run_opb("../test/data/pti/case5_alc.raw", NFAPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 1000.0; atol = 1e0) end @testset "5-bus with negative generators" begin - result = run_opb("../test/data/matpower/case5_npg.m", NFAPowerModel, ipopt_solver) + result = run_opb("../test/data/matpower/case5_npg.m", NFAPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], -6844.0; atol = 1e0) end @testset "5-bus with only current limit data" begin - result = run_opb("../test/data/matpower/case5_clm.m", NFAPowerModel, ipopt_solver) + result = run_opb("../test/data/matpower/case5_clm.m", NFAPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 14810.0; atol = 1e0) end @testset "5-bus with pwl costs" begin - result = run_opb("../test/data/matpower/case5_pwlc.m", NFAPowerModel, ipopt_solver) + result = run_opb("../test/data/matpower/case5_pwlc.m", NFAPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 42565.8; atol = 1e0) end @testset "6-bus case" begin - result = run_opb("../test/data/matpower/case6.m", NFAPowerModel, ipopt_solver) + result = run_opb("../test/data/matpower/case6.m", NFAPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 11277.9; atol = 1e0) end @testset "24-bus rts case" begin - result = run_opb("../test/data/matpower/case24.m", NFAPowerModel, ipopt_solver) + result = run_opb("../test/data/matpower/case24.m", NFAPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 61001.2; atol = 1e0) @@ -72,25 +72,25 @@ end @testset "test dcp opb" begin @testset "3-bus case" begin - result = run_opb("../test/data/matpower/case3.m", DCPPowerModel, ipopt_solver) + result = run_opb("../test/data/matpower/case3.m", DCPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 5638.97; atol = 1e0) end @testset "5-bus with pwl costs" begin - result = run_opb("../test/data/matpower/case5_pwlc.m", DCPPowerModel, ipopt_solver) + result = run_opb("../test/data/matpower/case5_pwlc.m", DCPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 42565.8; atol = 1e0) end @testset "6-bus case" begin - result = run_opb("../test/data/matpower/case6.m", DCPPowerModel, ipopt_solver) + result = run_opb("../test/data/matpower/case6.m", DCPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 11277.9; atol = 1e0) end @testset "24-bus case" begin - result = run_opb("../test/data/matpower/case24.m", DCPPowerModel, ipopt_solver) + result = run_opb("../test/data/matpower/case24.m", DCPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 61001.2; atol = 1e0) @@ -100,25 +100,25 @@ end @testset "test soc wr opb" begin @testset "3-bus case" begin - result = run_opb("../test/data/matpower/case3.m", SOCWRPowerModel, ipopt_solver) + result = run_opb("../test/data/matpower/case3.m", SOCWRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 5638.97; atol = 1e0) end @testset "5-bus with pwl costs" begin - result = run_opb("../test/data/matpower/case5_pwlc.m", SOCWRPowerModel, ipopt_solver) + result = run_opb("../test/data/matpower/case5_pwlc.m", SOCWRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 42565.7; atol = 1e0) end @testset "6-bus case" begin - result = run_opb("../test/data/matpower/case6.m", SOCWRPowerModel, ipopt_solver) + result = run_opb("../test/data/matpower/case6.m", SOCWRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 11277.9; atol = 1e0) end @testset "24-bus case" begin - result = run_opb("../test/data/matpower/case24.m", SOCWRPowerModel, ipopt_solver) + result = run_opb("../test/data/matpower/case24.m", SOCWRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 61001.2; atol = 1e0) diff --git a/test/opf-obj.jl b/test/opf-obj.jl index 9f56914d7..f0ba6f69a 100644 --- a/test/opf-obj.jl +++ b/test/opf-obj.jl @@ -7,21 +7,21 @@ data["gen"]["5"]["cost"] = [] @testset "nlp solver" begin - result = run_ac_opf(data, ipopt_solver) + result = run_ac_opf(data, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 5420.3; atol = 1e0) end @testset "conic solver" begin - result = run_opf(data, SOCWRConicPowerModel, scs_solver) + result = run_opf(data, SOCWRConicPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 3095.88; atol = 1e0) end @testset "lp solver" begin - result = run_dc_opf(data, cbc_solver) + result = run_dc_opf(data, milp_solver) @test result["termination_status"] == OPTIMAL # @test isapprox(result["objective"], 4679.05; atol = 1e0) # Problem upstream with JuMP.SecondOrderCone or JuMP.RotatedSecondOrderCone? @@ -35,14 +35,14 @@ end data["gen"]["5"]["cost"] = [] @testset "nlp solver" begin - result = run_ac_opf(data, ipopt_solver) + result = run_ac_opf(data, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 5420.46; atol = 1e0) end @testset "conic solver" begin - result = run_opf(data, SOCWRConicPowerModel, scs_solver) + result = run_opf(data, SOCWRConicPowerModel, sdp_solver) # ALMOST_OPTIMAL only required for Julia 1.0 on Linux @test result["termination_status"] == OPTIMAL || result["termination_status"] == ALMOST_OPTIMAL @@ -58,14 +58,14 @@ end data["gen"]["5"]["cost"] = [] @testset "opf objective" begin - result = run_ac_opf(data, ipopt_solver) + result = run_ac_opf(data, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 5452.04; atol = 1e0) end @testset "opb objective" begin - result = run_opb(data, SOCWRPowerModel, ipopt_solver) + result = run_opb(data, SOCWRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 2999.80; atol = 1e0) @@ -79,7 +79,7 @@ end @testset "nlp objective" begin data["dcline"]["1"]["cost"] = [100.0, 300.0, 4000.0, 1.0] - result = run_ac_opf(data, ipopt_solver) + result = run_ac_opf(data, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 18160.3; atol = 1e0) @@ -87,7 +87,7 @@ end @testset "qp objective" begin data["dcline"]["1"]["cost"] = [1.0, 4000.0, 1.0] - result = run_ac_opf(data, ipopt_solver) + result = run_ac_opf(data, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 18157.2; atol = 1e0) @@ -95,7 +95,7 @@ end @testset "linear objective" begin data["dcline"]["1"]["cost"] = [4000.0, 1.0] - result = run_ac_opf(data, ipopt_solver) + result = run_ac_opf(data, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 18157.2; atol = 1e0) @@ -103,7 +103,7 @@ end @testset "constant objective" begin data["dcline"]["1"]["cost"] = [1.0] - result = run_ac_opf(data, ipopt_solver) + result = run_ac_opf(data, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 17757.2; atol = 1e0) @@ -111,7 +111,7 @@ end @testset "empty objective" begin data["dcline"]["1"]["cost"] = [] - result = run_ac_opf(data, ipopt_solver) + result = run_ac_opf(data, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 17756.2; atol = 1e0) diff --git a/test/opf-ptdf.jl b/test/opf-ptdf.jl index 71dcb0d25..c72129ae1 100644 --- a/test/opf-ptdf.jl +++ b/test/opf-ptdf.jl @@ -4,8 +4,8 @@ TESTLOG = Memento.getlogger(PowerModels) @testset "test ptdf-based dc opf" begin @testset "5-bus case, LP solver" begin data = PowerModels.parse_file("../test/data/matpower/case5.m") - result_va = PowerModels.run_opf(data, DCPPowerModel, cbc_solver) - result_pg = PowerModels.run_opf_ptdf(data, DCPPowerModel, cbc_solver) + result_va = PowerModels.run_opf(data, DCPPowerModel, milp_solver) + result_pg = PowerModels.run_opf_ptdf(data, DCPPowerModel, milp_solver) @test result_va["termination_status"] == OPTIMAL @test result_pg["termination_status"] == OPTIMAL @@ -16,8 +16,8 @@ TESTLOG = Memento.getlogger(PowerModels) end @testset "5-bus gap case" begin data = PowerModels.parse_file("../test/data/matpower/case5_gap.m") - result_va = PowerModels.run_opf(data, DCPPowerModel, ipopt_solver) - result_pg = PowerModels.run_opf_ptdf(data, DCPPowerModel, ipopt_solver) + result_va = PowerModels.run_opf(data, DCPPowerModel, nlp_solver) + result_pg = PowerModels.run_opf_ptdf(data, DCPPowerModel, nlp_solver) @test result_va["termination_status"] == LOCALLY_SOLVED @test result_pg["termination_status"] == LOCALLY_SOLVED @@ -29,8 +29,8 @@ TESTLOG = Memento.getlogger(PowerModels) @testset "5-bus with pwl costs" begin data = PowerModels.parse_file("../test/data/matpower/case5_pwlc.m") data["dcline"] = Dict() - result_va = PowerModels.run_opf(data, DCPPowerModel, ipopt_solver) - result_pg = PowerModels.run_opf_ptdf(data, DCPPowerModel, ipopt_solver) + result_va = PowerModels.run_opf(data, DCPPowerModel, nlp_solver) + result_pg = PowerModels.run_opf_ptdf(data, DCPPowerModel, nlp_solver) @test result_va["termination_status"] == LOCALLY_SOLVED @test result_pg["termination_status"] == LOCALLY_SOLVED @@ -41,8 +41,8 @@ TESTLOG = Memento.getlogger(PowerModels) end @testset "14-bus case" begin data = PowerModels.parse_file("../test/data/matpower/case14.m") - result_va = PowerModels.run_opf(data, DCPPowerModel, ipopt_solver) - result_pg = PowerModels.run_opf_ptdf(data, DCPPowerModel, ipopt_solver) + result_va = PowerModels.run_opf(data, DCPPowerModel, nlp_solver) + result_pg = PowerModels.run_opf_ptdf(data, DCPPowerModel, nlp_solver) @test result_va["termination_status"] == LOCALLY_SOLVED @test result_pg["termination_status"] == LOCALLY_SOLVED @@ -55,19 +55,19 @@ TESTLOG = Memento.getlogger(PowerModels) @testset "no support for zero reference buses" begin data = PowerModels.parse_file("../test/data/matpower/case5.m") data["bus"]["4"]["bus_type"] = 2 - @test_throws(TESTLOG, ErrorException, PowerModels.run_opf_ptdf(data, DCPPowerModel, ipopt_solver)) + @test_throws(TESTLOG, ErrorException, PowerModels.run_opf_ptdf(data, DCPPowerModel, nlp_solver)) end @testset "no support for multiple connected components" begin - @test_throws(TESTLOG, ErrorException, PowerModels.run_opf_ptdf("../test/data/matpower/case6.m", DCPPowerModel, ipopt_solver)) + @test_throws(TESTLOG, ErrorException, PowerModels.run_opf_ptdf("../test/data/matpower/case6.m", DCPPowerModel, nlp_solver)) end @testset "no support for dclines" begin - @test_throws(TESTLOG, ErrorException, PowerModels.run_opf_ptdf("../test/data/matpower/case3.m", DCPPowerModel, ipopt_solver)) + @test_throws(TESTLOG, ErrorException, PowerModels.run_opf_ptdf("../test/data/matpower/case3.m", DCPPowerModel, nlp_solver)) end @testset "no support for switches" begin - @test_throws(TESTLOG, ErrorException, PowerModels.run_opf_ptdf("../test/data/matpower/case5_sw.m", DCPPowerModel, ipopt_solver)) + @test_throws(TESTLOG, ErrorException, PowerModels.run_opf_ptdf("../test/data/matpower/case5_sw.m", DCPPowerModel, nlp_solver)) end end diff --git a/test/opf-var.jl b/test/opf-var.jl index 136042b3e..17f223fe2 100644 --- a/test/opf-var.jl +++ b/test/opf-var.jl @@ -16,20 +16,20 @@ end @testset "test ac polar opf" begin @testset "5-bus case" begin data = build_current_data("../test/data/matpower/case5.m") - result = PowerModels._solve_opf_cl(data, ACPPowerModel, ipopt_solver) + result = PowerModels._solve_opf_cl(data, ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 15669.8; atol = 1e0) end @testset "5-bus current limit case" begin - result = PowerModels._solve_opf_cl("../test/data/matpower/case5_clm.m", ACPPowerModel, ipopt_solver) + result = PowerModels._solve_opf_cl("../test/data/matpower/case5_clm.m", ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 17015.5; atol = 1e0) end @testset "14-bus no limits case" begin data = build_current_data("../test/data/matpower/case14.m") - result = PowerModels._solve_opf_cl(data, ACPPowerModel, ipopt_solver) + result = PowerModels._solve_opf_cl(data, ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 8081.52; atol = 1e0) @@ -39,14 +39,14 @@ end @testset "test ac rect opf" begin @testset "5-bus case" begin data = build_current_data("../test/data/matpower/case5.m") - result = PowerModels._solve_opf_cl(data, ACRPowerModel, ipopt_solver) + result = PowerModels._solve_opf_cl(data, ACRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 15669.8; atol = 1e0) end @testset "14-bus no limits case" begin data = build_current_data("../test/data/matpower/case14.m") - result = PowerModels._solve_opf_cl(data, ACRPowerModel, ipopt_solver) + result = PowerModels._solve_opf_cl(data, ACRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 8081.52; atol = 1e0) @@ -56,14 +56,14 @@ end @testset "test ac tan opf" begin @testset "5-bus case" begin data = build_current_data("../test/data/matpower/case5.m") - result = PowerModels._solve_opf_cl(data, ACTPowerModel, ipopt_solver) + result = PowerModels._solve_opf_cl(data, ACTPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 15669.8; atol = 1e0) end @testset "14-bus no limits case" begin data = build_current_data("../test/data/matpower/case14.m") - result = PowerModels._solve_opf_cl(data, ACTPowerModel, ipopt_solver) + result = PowerModels._solve_opf_cl(data, ACTPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 8081.52; atol = 1e0) @@ -73,14 +73,14 @@ end @testset "test dc opf" begin @testset "5-bus case" begin data = build_current_data("../test/data/matpower/case5.m") - result = PowerModels._solve_opf_cl(data, DCPPowerModel, ipopt_solver) + result = PowerModels._solve_opf_cl(data, DCPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 16154.5; atol = 1e0) end @testset "14-bus no limits case" begin data = build_current_data("../test/data/matpower/case14.m") - result = PowerModels._solve_opf_cl(data, DCPPowerModel, ipopt_solver) + result = PowerModels._solve_opf_cl(data, DCPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 7642.59; atol = 1e0) @@ -90,14 +90,14 @@ end @testset "test dc+ll opf" begin @testset "5-bus case" begin data = build_current_data("../test/data/matpower/case5.m") - result = PowerModels._solve_opf_cl(data, DCPLLPowerModel, ipopt_solver) + result = PowerModels._solve_opf_cl(data, DCPLLPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 16282.6; atol = 1e0) end @testset "14-bus no limits case" begin data = build_current_data("../test/data/matpower/case14.m") - result = PowerModels._solve_opf_cl(data, DCPLLPowerModel, ipopt_solver) + result = PowerModels._solve_opf_cl(data, DCPLLPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 8110.61; atol = 1e0) @@ -107,14 +107,14 @@ end @testset "test soc (BIM) opf" begin @testset "5-bus case" begin data = build_current_data("../test/data/matpower/case5.m") - result = PowerModels._solve_opf_cl(data, SOCWRPowerModel, ipopt_solver) + result = PowerModels._solve_opf_cl(data, SOCWRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 15047.7; atol = 1e0) end @testset "14-bus no limits case" begin data = build_current_data("../test/data/matpower/case14.m") - result = PowerModels._solve_opf_cl(data, SOCWRPowerModel, ipopt_solver) + result = PowerModels._solve_opf_cl(data, SOCWRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 8075.12; atol = 1e0) @@ -124,14 +124,14 @@ end @testset "test sdp opf" begin @testset "3-bus case" begin data = build_current_data("../test/data/matpower/case3.m") - result = PowerModels._solve_opf_cl(data, SDPWRMPowerModel, scs_solver) + result = PowerModels._solve_opf_cl(data, SDPWRMPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 5728.62; atol = 1e0) end @testset "5-bus case" begin data = build_current_data("../test/data/matpower/case5.m") - result = PowerModels._solve_opf_cl(data, SDPWRMPowerModel, scs_solver) + result = PowerModels._solve_opf_cl(data, SDPWRMPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL #@test isapprox(result["objective"], 15418.4; atol = 1e0) @@ -140,7 +140,7 @@ end end @testset "14-bus case" begin data = build_current_data("../test/data/matpower/case14.m") - result = PowerModels._solve_opf_cl(data, SDPWRMPowerModel, scs_solver) + result = PowerModels._solve_opf_cl(data, SDPWRMPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 7505.33; atol = 1e0) @@ -155,7 +155,7 @@ end @testset "test ac polar mld" begin @testset "5-bus case" begin - result = PowerModels._solve_mld("../test/data/matpower/case5.m", ACPPowerModel, ipopt_solver) + result = PowerModels._solve_mld("../test/data/matpower/case5.m", ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 10.0; atol = 1e-2) @@ -165,7 +165,7 @@ end end @testset "14-bus current" begin - result = PowerModels._solve_mld("../test/data/matpower/case14.m", ACPPowerModel, ipopt_solver) + result = PowerModels._solve_mld("../test/data/matpower/case14.m", ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 3.59; atol = 1e-2) @@ -183,7 +183,7 @@ end @testset "test ac rect mld" begin @testset "5-bus case" begin - result = PowerModels._solve_mld("../test/data/matpower/case5.m", ACRPowerModel, ipopt_solver) + result = PowerModels._solve_mld("../test/data/matpower/case5.m", ACRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 10.0; atol = 1e-2) @@ -192,7 +192,7 @@ end @test all_shunts_on(result) end @testset "14-bus case" begin - result = PowerModels._solve_mld("../test/data/matpower/case14.m", ACRPowerModel, ipopt_solver) + result = PowerModels._solve_mld("../test/data/matpower/case14.m", ACRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 3.59; atol = 1e-2) @@ -204,7 +204,7 @@ end @testset "test ac tan mld" begin @testset "5-bus case" begin - result = PowerModels._solve_mld("../test/data/matpower/case5.m", ACTPowerModel, ipopt_solver) + result = PowerModels._solve_mld("../test/data/matpower/case5.m", ACTPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 10.0; atol = 1e-2) @@ -213,7 +213,7 @@ end @test all_shunts_on(result) end @testset "14-bus case" begin - result = PowerModels._solve_mld("../test/data/matpower/case14.m", ACTPowerModel, ipopt_solver) + result = PowerModels._solve_mld("../test/data/matpower/case14.m", ACTPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 3.59; atol = 1e-2) @@ -225,7 +225,7 @@ end @testset "test nfa mld" begin @testset "5-bus case" begin - result = PowerModels._solve_mld("../test/data/matpower/case5.m", NFAPowerModel, ipopt_solver) + result = PowerModels._solve_mld("../test/data/matpower/case5.m", NFAPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 10.0; atol = 1e-2) @@ -234,7 +234,7 @@ end @test all_shunts_on(result) end @testset "14-bus case" begin - result = PowerModels._solve_mld("../test/data/matpower/case14.m", NFAPowerModel, ipopt_solver) + result = PowerModels._solve_mld("../test/data/matpower/case14.m", NFAPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 3.59; atol = 1e-2) @@ -246,7 +246,7 @@ end @testset "test dc mld" begin @testset "5-bus case" begin - result = PowerModels._solve_mld("../test/data/matpower/case5.m", DCPPowerModel, ipopt_solver) + result = PowerModels._solve_mld("../test/data/matpower/case5.m", DCPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 10.0; atol = 1e-2) @@ -255,7 +255,7 @@ end @test all_shunts_on(result) end @testset "14-bus case" begin - result = PowerModels._solve_mld("../test/data/matpower/case14.m", DCPPowerModel, ipopt_solver) + result = PowerModels._solve_mld("../test/data/matpower/case14.m", DCPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 3.59; atol = 1e-2) @@ -267,7 +267,7 @@ end @testset "test soc (BIM) mld" begin @testset "5-bus case" begin - result = PowerModels._solve_mld("../test/data/matpower/case5.m", SOCWRPowerModel, ipopt_solver) + result = PowerModels._solve_mld("../test/data/matpower/case5.m", SOCWRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 10.0; atol = 1e-2) @@ -276,7 +276,7 @@ end @test all_shunts_on(result) end @testset "14-bus case" begin - result = PowerModels._solve_mld("../test/data/matpower/case14.m", SOCWRPowerModel, ipopt_solver) + result = PowerModels._solve_mld("../test/data/matpower/case14.m", SOCWRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 3.59; atol = 1e-2) @@ -288,7 +288,7 @@ end @testset "test soc conic (BIM) mld" begin @testset "5-bus case" begin - result = PowerModels._solve_mld("../test/data/matpower/case5.m", SOCWRConicPowerModel, scs_solver) + result = PowerModels._solve_mld("../test/data/matpower/case5.m", SOCWRConicPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 10.0; atol = 1e-2) @@ -297,7 +297,7 @@ end @test all_shunts_on(result) end @testset "14-bus case" begin - result = PowerModels._solve_mld("../test/data/matpower/case14.m", SOCWRConicPowerModel, scs_solver) + result = PowerModels._solve_mld("../test/data/matpower/case14.m", SOCWRConicPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 3.59; atol = 1e-2) @@ -309,7 +309,7 @@ end @testset "test sdp mld" begin @testset "5-bus case" begin - result = PowerModels._solve_mld("../test/data/matpower/case5.m", SOCWRConicPowerModel, scs_solver) + result = PowerModels._solve_mld("../test/data/matpower/case5.m", SOCWRConicPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 10.0; atol = 1e-2) @@ -318,7 +318,7 @@ end @test all_shunts_on(result) end @testset "14-bus case" begin - result = PowerModels._solve_mld("../test/data/matpower/case14.m", SOCWRConicPowerModel, scs_solver) + result = PowerModels._solve_mld("../test/data/matpower/case14.m", SOCWRConicPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 3.59; atol = 1e-2) @@ -333,7 +333,7 @@ end settings = Dict("output" => Dict("duals" => true)) @testset "ac case" begin - result = PowerModels._solve_mld("../test/data/matpower/case5.m", ACPPowerModel, ipopt_solver, setting=settings) + result = PowerModels._solve_mld("../test/data/matpower/case5.m", ACPPowerModel, nlp_solver, setting=settings) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 10.0; atol = 1e-2) @@ -346,7 +346,7 @@ end end @testset "soc case" begin - result = PowerModels._solve_mld("../test/data/matpower/case5.m", SOCWRPowerModel, ipopt_solver, setting=settings) + result = PowerModels._solve_mld("../test/data/matpower/case5.m", SOCWRPowerModel, nlp_solver, setting=settings) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 10.0; atol = 1e-2) @@ -359,7 +359,7 @@ end end @testset "dc case" begin - result = PowerModels._solve_mld("../test/data/matpower/case5.m", SOCWRPowerModel, ipopt_solver, setting=settings) + result = PowerModels._solve_mld("../test/data/matpower/case5.m", SOCWRPowerModel, nlp_solver, setting=settings) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 10.0; atol = 1e-2) @@ -378,7 +378,7 @@ end @testset "test ac opf" begin @testset "5-bus uc case" begin - result = PowerModels._solve_ucopf("../test/data/matpower/case5_uc.m", ACPPowerModel, juniper_solver) + result = PowerModels._solve_ucopf("../test/data/matpower/case5_uc.m", ACPPowerModel, minlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 18270.0; atol = 1e0) @@ -388,7 +388,7 @@ end @testset "test soc opf" begin @testset "5-bus uc case" begin - result = PowerModels._solve_ucopf("../test/data/matpower/case5_uc.m", SOCWRPowerModel, juniper_solver) + result = PowerModels._solve_ucopf("../test/data/matpower/case5_uc.m", SOCWRPowerModel, minlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 15057.09; atol = 1e0) @@ -398,7 +398,7 @@ end @testset "test dc opf" begin @testset "5-bus uc case" begin - result = PowerModels._solve_ucopf("../test/data/matpower/case5_uc.m", DCPPowerModel, cbc_solver) + result = PowerModels._solve_ucopf("../test/data/matpower/case5_uc.m", DCPPowerModel, milp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 17613.2; atol = 1e0) @@ -410,7 +410,7 @@ end for (i,load) in data["load"] load["pd"] = 0.5*load["pd"] end - result = PowerModels._solve_ucopf(data, DCPPowerModel, cbc_solver) + result = PowerModels._solve_ucopf(data, DCPPowerModel, milp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 8008.0; atol = 1e0) @@ -422,7 +422,7 @@ end @testset "test ac opf" begin @testset "5-bus uc storage case" begin - result = PowerModels._solve_ucopf("../test/data/matpower/case5_uc_strg.m", ACPPowerModel, juniper_solver) + result = PowerModels._solve_ucopf("../test/data/matpower/case5_uc_strg.m", ACPPowerModel, minlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 17740.9; atol = 1e0) @@ -435,7 +435,7 @@ end @testset "test soc opf" begin @testset "5-bus uc storage case" begin - result = PowerModels._solve_ucopf("../test/data/matpower/case5_uc_strg.m", SOCWRPowerModel, juniper_solver) + result = PowerModels._solve_ucopf("../test/data/matpower/case5_uc_strg.m", SOCWRPowerModel, minlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 14525.0; atol = 1e0) @@ -448,7 +448,7 @@ end @testset "test dc opf" begin @testset "5-bus uc storage case" begin - result = PowerModels._solve_ucopf("../test/data/matpower/case5_uc_strg.m", DCPPowerModel, cbc_solver) + result = PowerModels._solve_ucopf("../test/data/matpower/case5_uc_strg.m", DCPPowerModel, milp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 16833.2; atol = 1e0) @@ -467,7 +467,7 @@ end @testset "test ac opf" begin @testset "5-bus sw case" begin - result = PowerModels._solve_opf_sw("../test/data/matpower/case5_sw.m", ACPPowerModel, ipopt_solver) + result = PowerModels._solve_opf_sw("../test/data/matpower/case5_sw.m", ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 16641.2; atol = 1e0) @@ -480,7 +480,7 @@ end end @testset "5-bus sw nb case" begin - result = PowerModels._solve_opf_sw("../test/data/matpower/case5_sw_nb.m", ACPPowerModel, ipopt_solver) + result = PowerModels._solve_opf_sw("../test/data/matpower/case5_sw_nb.m", ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 17915.3; atol = 1e0) @@ -489,7 +489,7 @@ end @testset "test dc opf" begin @testset "5-bus sw case" begin - result = PowerModels._solve_opf_sw("../test/data/matpower/case5_sw.m", DCPPowerModel, cbc_solver) + result = PowerModels._solve_opf_sw("../test/data/matpower/case5_sw.m", DCPPowerModel, milp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 16554.7; atol = 1e0) @@ -500,7 +500,7 @@ end end @testset "5-bus sw nb case" begin - result = PowerModels._solve_opf_sw("../test/data/matpower/case5_sw_nb.m", DCPPowerModel, cbc_solver) + result = PowerModels._solve_opf_sw("../test/data/matpower/case5_sw_nb.m", DCPPowerModel, milp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 17751.3; atol = 1e0) @@ -509,7 +509,7 @@ end @testset "test soc opf" begin @testset "5-bus sw case" begin - result = PowerModels._solve_opf_sw("../test/data/matpower/case5_sw.m", SOCWRPowerModel, ipopt_solver) + result = PowerModels._solve_opf_sw("../test/data/matpower/case5_sw.m", SOCWRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 15110.0; atol = 1e0) @@ -520,7 +520,7 @@ end end @testset "5-bus sw nb case" begin - result = PowerModels._solve_opf_sw("../test/data/matpower/case5_sw_nb.m", SOCWRPowerModel, ipopt_solver) + result = PowerModels._solve_opf_sw("../test/data/matpower/case5_sw_nb.m", SOCWRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 15192.8; atol = 1e0) @@ -534,7 +534,7 @@ end @testset "test ac oswpf" begin @testset "5-bus sw case" begin - result = PowerModels._solve_oswpf("../test/data/matpower/case5_sw.m", ACPPowerModel, juniper_solver) + result = PowerModels._solve_oswpf("../test/data/matpower/case5_sw.m", ACPPowerModel, minlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 15053.6; atol = 1e0) @@ -548,7 +548,7 @@ end end @testset "5-bus sw nb case" begin - result = PowerModels._solve_oswpf("../test/data/matpower/case5_sw_nb.m", ACPPowerModel, juniper_solver) + result = PowerModels._solve_oswpf("../test/data/matpower/case5_sw_nb.m", ACPPowerModel, minlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 17915.2; atol = 1e0) @@ -560,7 +560,7 @@ end @testset "test dc oswpf" begin @testset "5-bus sw case" begin - result = PowerModels._solve_oswpf("../test/data/matpower/case5_sw.m", DCPPowerModel, cbc_solver) + result = PowerModels._solve_oswpf("../test/data/matpower/case5_sw.m", DCPPowerModel, milp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 15054.1; atol = 1e0) @@ -575,7 +575,7 @@ end end @testset "5-bus sw nb case" begin - result = PowerModels._solve_oswpf("../test/data/matpower/case5_sw_nb.m", DCPPowerModel, cbc_solver) + result = PowerModels._solve_oswpf("../test/data/matpower/case5_sw_nb.m", DCPPowerModel, milp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 15141.2; atol = 1e0) @@ -587,7 +587,7 @@ end @testset "test soc oswpf" begin @testset "5-bus sw case" begin - result = PowerModels._solve_oswpf("../test/data/matpower/case5_sw.m", SOCWRPowerModel, juniper_solver) + result = PowerModels._solve_oswpf("../test/data/matpower/case5_sw.m", SOCWRPowerModel, minlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 15053.6; atol = 1e0) @@ -601,7 +601,7 @@ end end @testset "5-bus sw nb case" begin - result = PowerModels._solve_oswpf("../test/data/matpower/case5_sw_nb.m", SOCWRPowerModel, juniper_solver) + result = PowerModels._solve_oswpf("../test/data/matpower/case5_sw_nb.m", SOCWRPowerModel, minlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 15175.7; atol = 1e0) @@ -619,7 +619,7 @@ end @testset "test ac oswpf node-breaker" begin @testset "5-bus sw nb case" begin - result = PowerModels._solve_oswpf_nb("../test/data/matpower/case5_sw_nb.m", ACPPowerModel, juniper_solver) + result = PowerModels._solve_oswpf_nb("../test/data/matpower/case5_sw_nb.m", ACPPowerModel, minlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 15350.4; atol = 1e0) @@ -634,7 +634,7 @@ end @testset "test dc oswpf node-breaker" begin @testset "5-bus sw nb case" begin - result = PowerModels._solve_oswpf_nb("../test/data/matpower/case5_sw_nb.m", DCPPowerModel, cbc_solver) + result = PowerModels._solve_oswpf_nb("../test/data/matpower/case5_sw_nb.m", DCPPowerModel, milp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 15141.2; atol = 1e0) @@ -649,7 +649,7 @@ end @testset "test soc oswpf node-breaker" begin @testset "5-bus sw nb case" begin - result = PowerModels._solve_oswpf_nb("../test/data/matpower/case5_sw_nb.m", SOCWRPowerModel, juniper_solver) + result = PowerModels._solve_oswpf_nb("../test/data/matpower/case5_sw_nb.m", SOCWRPowerModel, minlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 15175.7; atol = 1e0) @@ -669,7 +669,7 @@ end @testset "test acp polar opf" begin @testset "5-bus case" begin - result = PowerModels._solve_opf_strg("../test/data/matpower/case5_strg.m", PowerModels.ACPPowerModel, ipopt_solver) + result = PowerModels._solve_opf_strg("../test/data/matpower/case5_strg.m", PowerModels.ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 17039.7; atol = 1e0) @@ -683,7 +683,7 @@ end @testset "test mi acp polar opf" begin @testset "5-bus case" begin - result = PowerModels._solve_opf_strg_mi("../test/data/matpower/case5_strg.m", PowerModels.ACPPowerModel, juniper_solver) + result = PowerModels._solve_opf_strg_mi("../test/data/matpower/case5_strg.m", PowerModels.ACPPowerModel, minlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 17039.7; atol = 1e0) @@ -698,7 +698,7 @@ end @testset "test acr polar opf" begin @testset "5-bus case" begin - result = PowerModels._solve_opf_strg("../test/data/matpower/case5_strg.m", PowerModels.ACRPowerModel, ipopt_solver) + result = PowerModels._solve_opf_strg("../test/data/matpower/case5_strg.m", PowerModels.ACRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 17039.7; atol = 1e0) @@ -712,7 +712,7 @@ end @testset "test mi acr polar opf" begin @testset "5-bus case" begin - result = PowerModels._solve_opf_strg_mi("../test/data/matpower/case5_strg.m", PowerModels.ACRPowerModel, juniper_solver) + result = PowerModels._solve_opf_strg_mi("../test/data/matpower/case5_strg.m", PowerModels.ACRPowerModel, minlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 17039.7; atol = 1e0) @@ -727,7 +727,7 @@ end @testset "test soc polar opf" begin @testset "5-bus case" begin - result = PowerModels._solve_opf_strg("../test/data/matpower/case5_strg.m", PowerModels.SOCWRPowerModel, ipopt_solver) + result = PowerModels._solve_opf_strg("../test/data/matpower/case5_strg.m", PowerModels.SOCWRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 13799.5; atol = 1e0) @@ -741,7 +741,7 @@ end @testset "test mi soc polar opf" begin @testset "5-bus case" begin - result = PowerModels._solve_opf_strg_mi("../test/data/matpower/case5_strg.m", PowerModels.SOCWRPowerModel, juniper_solver) + result = PowerModels._solve_opf_strg_mi("../test/data/matpower/case5_strg.m", PowerModels.SOCWRPowerModel, minlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 13799.5; atol = 1e0) @@ -756,7 +756,7 @@ end @testset "test dc opf" begin @testset "5-bus case" begin - result = PowerModels._solve_opf_strg("../test/data/matpower/case5_strg.m", PowerModels.DCPPowerModel, ipopt_solver) + result = PowerModels._solve_opf_strg("../test/data/matpower/case5_strg.m", PowerModels.DCPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 16840.7; atol = 1e0) @@ -770,7 +770,7 @@ end @testset "test mi dc opf" begin @testset "5-bus case" begin - result = PowerModels._solve_opf_strg_mi("../test/data/matpower/case5_strg.m", PowerModels.DCPPowerModel, cbc_solver) + result = PowerModels._solve_opf_strg_mi("../test/data/matpower/case5_strg.m", PowerModels.DCPPowerModel, milp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 16840.7; atol = 1e0) @@ -785,7 +785,7 @@ end @testset "test dc+ll opf" begin @testset "5-bus case" begin - result = PowerModels._solve_opf_strg("../test/data/matpower/case5_strg.m", PowerModels.DCPLLPowerModel, ipopt_solver) + result = PowerModels._solve_opf_strg("../test/data/matpower/case5_strg.m", PowerModels.DCPLLPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 17048.4; atol = 1e0) @@ -801,7 +801,7 @@ end mp_data = PowerModels.parse_file("../test/data/matpower/case5_strg.m") delete!(mp_data, "time_elapsed") Memento.setlevel!(TESTLOG, "warn") - @test_warn(TESTLOG, "network data should specify time_elapsed, using 1.0 as a default", PowerModels._solve_opf_strg(mp_data, PowerModels.ACPPowerModel, ipopt_solver)) + @test_warn(TESTLOG, "network data should specify time_elapsed, using 1.0 as a default", PowerModels._solve_opf_strg(mp_data, PowerModels.ACPPowerModel, nlp_solver)) Memento.setlevel!(TESTLOG, "error") end @@ -846,31 +846,31 @@ end end @testset "3-bus case" begin - result = run_model("../test/data/matpower/case3.m", ACPPowerModel, ipopt_solver, build_opf_var) + result = run_model("../test/data/matpower/case3.m", ACPPowerModel, nlp_solver, build_opf_var) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 5907; atol = 1e0) end @testset "5-bus asymmetric case" begin - result = run_model("../test/data/matpower/case5_asym.m", ACPPowerModel, ipopt_solver, build_opf_var) + result = run_model("../test/data/matpower/case5_asym.m", ACPPowerModel, nlp_solver, build_opf_var) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 17551; atol = 1e0) end @testset "5-bus gap case" begin - result = run_model("../test/data/matpower/case5_gap.m", ACPPowerModel, ipopt_solver, build_opf_var) + result = run_model("../test/data/matpower/case5_gap.m", ACPPowerModel, nlp_solver, build_opf_var) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], -27497.7; atol = 1e0) end @testset "5-bus with dcline costs" begin - result = run_model("../test/data/matpower/case5_dc.m", ACPPowerModel, ipopt_solver, build_opf_var) + result = run_model("../test/data/matpower/case5_dc.m", ACPPowerModel, nlp_solver, build_opf_var) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 18156.2; atol = 1e0) end @testset "6-bus case" begin - result = run_model("../test/data/matpower/case6.m", ACPPowerModel, ipopt_solver, build_opf_var) + result = run_model("../test/data/matpower/case6.m", ACPPowerModel, nlp_solver, build_opf_var) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 11625.3; atol = 1e0) @@ -878,7 +878,7 @@ end @test isapprox(result["solution"]["bus"]["4"]["va"], 0.0; atol = 1e-4) end @testset "24-bus rts case" begin - result = run_model("../test/data/matpower/case24.m", ACPPowerModel, ipopt_solver, build_opf_var) + result = run_model("../test/data/matpower/case24.m", ACPPowerModel, nlp_solver, build_opf_var) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 79805; atol = 1e0) @@ -891,7 +891,7 @@ end @testset "3-bus case with fixed phase shift / tap" begin file = "../test/data/matpower/case3_oltc_pst.m" data = PowerModels.parse_file(file) - result = PowerModels.run_opf(data, ACPPowerModel, ipopt_solver) + result = PowerModels.run_opf(data, ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 5820.1; atol = 1e0) @@ -900,7 +900,7 @@ end @testset "3-bus case with optimal phase shifting / tap changing" begin file = "../test/data/matpower/case3_oltc_pst.m" data = PowerModels.parse_file(file) - result = PowerModels._solve_opf_oltc_pst(data, ACPPowerModel, ipopt_solver) + result = PowerModels._solve_opf_oltc_pst(data, ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 5738.6; atol = 1e0) @@ -926,7 +926,7 @@ end branch["tm_min"] = branch["tap"] branch["tm_max"] = branch["tap"] end - result = PowerModels._solve_opf_oltc_pst(data, ACPPowerModel, ipopt_solver) + result = PowerModels._solve_opf_oltc_pst(data, ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 5820.1; atol = 1e0) diff --git a/test/opf.jl b/test/opf.jl index 48a14184f..77282a89e 100644 --- a/test/opf.jl +++ b/test/opf.jl @@ -2,43 +2,43 @@ @testset "test ac polar opf" begin @testset "3-bus case" begin - result = run_ac_opf("../test/data/matpower/case3.m", ipopt_solver) + result = run_ac_opf("../test/data/matpower/case3.m", nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 5907; atol = 1e0) end @testset "5-bus tranformer swap case" begin - result = run_ac_opf("../test/data/matpower/case5.m", ipopt_solver) + result = run_ac_opf("../test/data/matpower/case5.m", nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 18269; atol = 1e0) end @testset "5-bus asymmetric case" begin - result = run_ac_opf("../test/data/matpower/case5_asym.m", ipopt_solver) + result = run_ac_opf("../test/data/matpower/case5_asym.m", nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 17551; atol = 1e0) end @testset "5-bus gap case" begin - result = run_ac_opf("../test/data/matpower/case5_gap.m", ipopt_solver) + result = run_ac_opf("../test/data/matpower/case5_gap.m", nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], -27497.7; atol = 1e0) end @testset "5-bus with dcline costs" begin - result = run_ac_opf("../test/data/matpower/case5_dc.m", ipopt_solver) + result = run_ac_opf("../test/data/matpower/case5_dc.m", nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 18156.2; atol = 1e0) end @testset "5-bus with asymmetric line charge" begin - result = run_ac_opf("../test/data/pti/case5_alc.raw", ipopt_solver) + result = run_ac_opf("../test/data/pti/case5_alc.raw", nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 1005.31; atol = 1e0) end @testset "5-bus with negative generators" begin - result = run_ac_opf("../test/data/matpower/case5_npg.m", ipopt_solver) + result = run_ac_opf("../test/data/matpower/case5_npg.m", nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 8190.09; atol = 1e0) @@ -46,31 +46,31 @@ @testset "5-bus with only current limit data" begin data = PowerModels.parse_file("../test/data/matpower/case5_clm.m") calc_thermal_limits!(data) - result = run_ac_opf(data, ipopt_solver) + result = run_ac_opf(data, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 16513.6; atol = 1e0) end @testset "5-bus with pwl costs" begin - result = run_ac_opf("../test/data/matpower/case5_pwlc.m", ipopt_solver) + result = run_ac_opf("../test/data/matpower/case5_pwlc.m", nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 42895; atol = 1e0) end @testset "5-bus with gen lb" begin - result = run_ac_opf("../test/data/matpower/case5_uc.m", ipopt_solver) + result = run_ac_opf("../test/data/matpower/case5_uc.m", nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 18742.2; atol = 1e0) end @testset "5-bus with dangling bus" begin - result = run_ac_opf("../test/data/matpower/case5_db.m", ipopt_solver) + result = run_ac_opf("../test/data/matpower/case5_db.m", nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 16739.1; atol = 1e0) end @testset "6-bus case" begin - result = run_ac_opf("../test/data/matpower/case6.m", ipopt_solver) + result = run_ac_opf("../test/data/matpower/case6.m", nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 11625.3; atol = 1e0) @@ -78,7 +78,7 @@ @test isapprox(result["solution"]["bus"]["4"]["va"], 0.0; atol = 1e-4) end @testset "24-bus rts case" begin - result = run_opf("../test/data/matpower/case24.m", ACPPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case24.m", ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 79805; atol = 1e0) @@ -94,7 +94,7 @@ end #= # numerical issue @testset "3-bus case" begin - result = run_opf("../test/data/matpower/case3.m", ACRPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case3.m", ACRPowerModel, nlp_solver) #@test result["termination_status"] == LOCALLY_SOLVED #@test isapprox(result["objective"], 5812; atol = 1e0) @@ -102,31 +102,31 @@ end end =# @testset "5-bus asymmetric case" begin - result = run_opf("../test/data/matpower/case5_asym.m", ACRPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case5_asym.m", ACRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 17551; atol = 1e0) end @testset "5-bus gap case" begin - result = run_opf("../test/data/matpower/case5_gap.m", ACRPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case5_gap.m", ACRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], -27497.7; atol = 1e0) end @testset "5-bus with asymmetric line charge" begin - result = run_opf("../test/data/pti/case5_alc.raw", ACRPowerModel, ipopt_solver) + result = run_opf("../test/data/pti/case5_alc.raw", ACRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 1005.31; atol = 1e0) end @testset "5-bus with pwl costs" begin - result = run_ac_opf("../test/data/matpower/case5_pwlc.m", ipopt_solver) + result = run_ac_opf("../test/data/matpower/case5_pwlc.m", nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 42895; atol = 1e0) end @testset "6-bus case" begin - result = run_opf("../test/data/matpower/case6.m", ACRPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case6.m", ACRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 11625.3; atol = 1e0) @@ -134,7 +134,7 @@ end @test isapprox(result["solution"]["bus"]["4"]["vi"], 0.0; atol = 1e-4) end @testset "24-bus rts case" begin - result = run_opf("../test/data/matpower/case24.m", ACRPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case24.m", ACRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 79805; atol = 1e0) @@ -148,37 +148,37 @@ end @testset "test ac tan opf" begin @testset "3-bus case" begin - result = run_opf("../test/data/matpower/case3.m", ACTPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case3.m", ACTPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 5907; atol = 1e0) end @testset "5-bus asymmetric case" begin - result = run_opf("../test/data/matpower/case5_asym.m", ACTPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case5_asym.m", ACTPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 17551; atol = 1e0) end @testset "5-bus gap case" begin - result = run_opf("../test/data/matpower/case5_gap.m", ACTPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case5_gap.m", ACTPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], -27438.7; atol = 1e0) end @testset "5-bus with asymmetric line charge" begin - result = run_opf("../test/data/pti/case5_alc.raw", ACTPowerModel, ipopt_solver) + result = run_opf("../test/data/pti/case5_alc.raw", ACTPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 1005.31; atol = 1e0) end @testset "5-bus with pwl costs" begin - result = run_opf("../test/data/matpower/case5_pwlc.m", ACTPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case5_pwlc.m", ACTPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 42895; atol = 1e0) end @testset "6-bus case" begin - result = run_opf("../test/data/matpower/case6.m", ACTPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case6.m", ACTPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 11625.3; atol = 1e0) @@ -186,7 +186,7 @@ end @test isapprox(result["solution"]["bus"]["4"]["va"], 0.0; atol = 1e-4) end @testset "24-bus rts case" begin - result = run_opf("../test/data/matpower/case24.m", ACTPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case24.m", ACTPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 79804; atol = 1e0) @@ -200,38 +200,38 @@ end @testset "test iv opf" begin @testset "3-bus case" begin - result = run_opf_iv("../test/data/matpower/case3.m", IVRPowerModel, ipopt_solver) + result = run_opf_iv("../test/data/matpower/case3.m", IVRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 5907; atol = 1e0) end @testset "5-bus asymmetric case" begin - result = run_opf_iv("../test/data/matpower/case5_asym.m", IVRPowerModel, ipopt_solver) + result = run_opf_iv("../test/data/matpower/case5_asym.m", IVRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 17551; atol = 1e0) end @testset "5-bus gap case" begin - result = run_opf_iv("../test/data/matpower/case5_gap.m", IVRPowerModel, ipopt_solver) + result = run_opf_iv("../test/data/matpower/case5_gap.m", IVRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], -27497.7; atol = 1e2) #numerically challenging , returns 27438.7 end @testset "5-bus with asymmetric line charge" begin - result = run_opf_iv("../test/data/pti/case5_alc.raw", IVRPowerModel, ipopt_solver) + result = run_opf_iv("../test/data/pti/case5_alc.raw", IVRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 1005.31; atol = 1e0) end @testset "5-bus with pwl costs" begin - result = run_opf_iv("../test/data/matpower/case5_pwlc.m", IVRPowerModel, ipopt_solver) + result = run_opf_iv("../test/data/matpower/case5_pwlc.m", IVRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 42895; atol = 1e0) end @testset "6-bus case" begin - result = run_opf_iv("../test/data/matpower/case6.m", IVRPowerModel, ipopt_solver) + result = run_opf_iv("../test/data/matpower/case6.m", IVRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 11625.3; atol = 1e0) @@ -239,7 +239,7 @@ end @test isapprox(result["solution"]["bus"]["4"]["vi"], 0.0; atol = 1e-4) end @testset "24-bus rts case" begin - result = run_opf_iv("../test/data/matpower/case24.m", IVRPowerModel, ipopt_solver) + result = run_opf_iv("../test/data/matpower/case24.m", IVRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 79804; atol = 1e0) @@ -254,55 +254,55 @@ end @testset "test dc opf" begin @testset "3-bus case" begin - result = run_dc_opf("../test/data/matpower/case3.m", ipopt_solver) + result = run_dc_opf("../test/data/matpower/case3.m", nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 5782; atol = 1e0) end @testset "5-bus case, LP solver" begin - result = run_dc_opf("../test/data/matpower/case5.m", cbc_solver) + result = run_dc_opf("../test/data/matpower/case5.m", milp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 17613; atol = 1e0) end @testset "5-bus asymmetric case" begin - result = run_dc_opf("../test/data/matpower/case5_asym.m", ipopt_solver) + result = run_dc_opf("../test/data/matpower/case5_asym.m", nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 17479; atol = 1e0) end @testset "5-bus gap case" begin - result = run_dc_opf("../test/data/matpower/case5_gap.m", ipopt_solver) + result = run_dc_opf("../test/data/matpower/case5_gap.m", nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], -27410.0; atol = 1e0) end @testset "5-bus with asymmetric line charge" begin - result = run_dc_opf("../test/data/pti/case5_alc.raw", ipopt_solver) + result = run_dc_opf("../test/data/pti/case5_alc.raw", nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 1000.0; atol = 1e0) end @testset "5-bus with pwl costs" begin - result = run_dc_opf("../test/data/matpower/case5_pwlc.m", ipopt_solver) + result = run_dc_opf("../test/data/matpower/case5_pwlc.m", nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 42565; atol = 1e0) end @testset "5-bus with gen lb" begin - result = run_dc_opf("../test/data/matpower/case5_uc.m", ipopt_solver) + result = run_dc_opf("../test/data/matpower/case5_uc.m", nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 18064.5; atol = 1e0) end @testset "5-bus with dangling bus" begin - result = run_dc_opf("../test/data/matpower/case5_db.m", ipopt_solver) + result = run_dc_opf("../test/data/matpower/case5_db.m", nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 16710.0; atol = 1e0) end @testset "6-bus case" begin - result = run_dc_opf("../test/data/matpower/case6.m", ipopt_solver) + result = run_dc_opf("../test/data/matpower/case6.m", nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 11391.8; atol = 1e0) @@ -311,7 +311,7 @@ end end # TODO verify this is really infeasible #@testset "24-bus rts case" begin - # result = run_opf("../test/data/matpower/case24.m", DCPPowerModel, ipopt_solver) + # result = run_opf("../test/data/matpower/case24.m", DCPPowerModel, nlp_solver) # @test result["termination_status"] == LOCALLY_SOLVED # @test isapprox(result["objective"], 79804; atol = 1e0) @@ -325,7 +325,7 @@ end @testset "test matpower dc opf" begin @testset "5-bus case with matpower DCMP model" begin - result = run_opf("../test/data/matpower/case5.m", DCMPPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case5.m", DCMPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @@ -340,43 +340,43 @@ end @testset "test nfa opf" begin @testset "3-bus case" begin - result = run_opf("../test/data/matpower/case3.m", NFAPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case3.m", NFAPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 5638.97; atol = 1e0) end @testset "5-bus asymmetric case" begin - result = run_opf("../test/data/matpower/case5_asym.m", NFAPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case5_asym.m", NFAPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 14810.0; atol = 1e0) end @testset "5-bus gap case" begin - result = run_opf("../test/data/matpower/case5_gap.m", NFAPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case5_gap.m", NFAPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], -27410.0; atol = 1e0) end @testset "5-bus with asymmetric line charge" begin - result = run_opf("../test/data/pti/case5_alc.raw", NFAPowerModel, ipopt_solver) + result = run_opf("../test/data/pti/case5_alc.raw", NFAPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 1000.0; atol = 1e0) end @testset "5-bus with pwl costs" begin - result = run_opf("../test/data/matpower/case5_pwlc.m", NFAPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case5_pwlc.m", NFAPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 42565.8; atol = 1e0) end @testset "6-bus case" begin - result = run_opf("../test/data/matpower/case6.m", NFAPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case6.m", NFAPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 11277.9; atol = 1e0) end @testset "24-bus rts case" begin - result = run_opf("../test/data/matpower/case24.m", NFAPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case24.m", NFAPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 61001.2; atol = 1e0) @@ -389,37 +389,37 @@ end @testset "test dc+ll opf" begin @testset "3-bus case" begin - result = run_opf("../test/data/matpower/case3.m", DCPLLPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case3.m", DCPLLPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 5885; atol = 1e0) end @testset "5-bus asymmetric case" begin - result = run_opf("../test/data/matpower/case5_asym.m", DCPLLPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case5_asym.m", DCPLLPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 17693; atol = 1e0) end @testset "5-bus gap case" begin - result = run_opf("../test/data/matpower/case5_gap.m", DCPLLPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case5_gap.m", DCPLLPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], -32710.0; atol = 1e0) end @testset "5-bus with asymmetric line charge" begin - result = run_opf("../test/data/pti/case5_alc.raw", DCPLLPowerModel, ipopt_solver) + result = run_opf("../test/data/pti/case5_alc.raw", DCPLLPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 1001.17; atol = 1e0) end @testset "5-bus with pwl costs" begin - result = run_opf("../test/data/matpower/case5_pwlc.m", DCPLLPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case5_pwlc.m", DCPLLPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 42937; atol = 1e0) end @testset "6-bus case" begin - result = run_opf("../test/data/matpower/case6.m", DCPLLPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case6.m", DCPLLPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 11574.3; atol = 1e0) @@ -427,7 +427,7 @@ end @test isapprox(result["solution"]["bus"]["4"]["va"], 0.0; atol = 1e-4) end @testset "24-bus rts case" begin - result = run_opf("../test/data/matpower/case24.m", DCPLLPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case24.m", DCPLLPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 82240; atol = 1e0) @@ -441,43 +441,43 @@ end @testset "test lpac-c opf" begin @testset "3-bus case" begin - result = run_opf("../test/data/matpower/case3.m", LPACCPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case3.m", LPACCPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 5908.98; atol = 1e0) end @testset "5-bus tranformer swap case" begin - result = run_opf("../test/data/matpower/case5.m", LPACCPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case5.m", LPACCPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 18288.1; atol = 1e0) end @testset "5-bus asymmetric case" begin - result = run_opf("../test/data/matpower/case5_asym.m", LPACCPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case5_asym.m", LPACCPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 17645.6; atol = 1e0) end @testset "5-bus gap case" begin - result = run_opf("../test/data/matpower/case5_gap.m", LPACCPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case5_gap.m", LPACCPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], -27554.5; atol = 1e0) end @testset "5-bus with dcline costs" begin - result = run_opf("../test/data/matpower/case5_dc.m", LPACCPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case5_dc.m", LPACCPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 18253.0; atol = 1e0) end @testset "5-bus with asymmetric line charge" begin - result = run_opf("../test/data/pti/case5_alc.raw", LPACCPowerModel, ipopt_solver) + result = run_opf("../test/data/pti/case5_alc.raw", LPACCPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 1004.58; atol = 1e0) end @testset "5-bus with negative generators" begin - result = run_opf("../test/data/matpower/case5_npg.m", LPACCPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case5_npg.m", LPACCPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 8082.54; atol = 1e0) @@ -485,19 +485,19 @@ end @testset "5-bus with only current limit data" begin data = PowerModels.parse_file("../test/data/matpower/case5_clm.m") calc_thermal_limits!(data) - result = run_opf(data, LPACCPowerModel, ipopt_solver) + result = run_opf(data, LPACCPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 16559.3; atol = 1e0) end @testset "5-bus with pwl costs" begin - result = run_opf("../test/data/matpower/case5_pwlc.m", LPACCPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case5_pwlc.m", LPACCPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 42853.4; atol = 1e0) end @testset "6-bus case" begin - result = run_opf("../test/data/matpower/case6.m", LPACCPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case6.m", LPACCPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 11615.1; atol = 1e0) @@ -506,7 +506,7 @@ end end # TODO uderstand why this is infeasible #@testset "24-bus rts case" begin - # result = run_opf("../test/data/matpower/case24.m", LPACCPowerModel, ipopt_solver) + # result = run_opf("../test/data/matpower/case24.m", LPACCPowerModel, nlp_solver) # @test result["termination_status"] == LOCALLY_SOLVED # @test isapprox(result["objective"], 79805; atol = 1e0) @@ -520,61 +520,61 @@ end @testset "test soc (BIM) opf" begin @testset "3-bus case" begin - result = run_opf("../test/data/matpower/case3.m", SOCWRPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case3.m", SOCWRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 5746.7; atol = 1e0) end @testset "5-bus transformer swap case" begin - result = run_opf("../test/data/matpower/case5.m", SOCWRPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case5.m", SOCWRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 15051; atol = 1e0) end @testset "5-bus asymmetric case" begin - result = run_opf("../test/data/matpower/case5_asym.m", SOCWRPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case5_asym.m", SOCWRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 14999; atol = 1e0) end @testset "5-bus gap case" begin - result = run_opf("../test/data/matpower/case5_gap.m", SOCWRPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case5_gap.m", SOCWRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], -28237.3; atol = 1e0) end @testset "5-bus with asymmetric line charge" begin - result = run_opf("../test/data/pti/case5_alc.raw", SOCWRPowerModel, ipopt_solver) + result = run_opf("../test/data/pti/case5_alc.raw", SOCWRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 1005.27; atol = 1e0) end @testset "5-bus with negative generators" begin - result = run_opf("../test/data/matpower/case5_npg.m", SOCWRPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case5_npg.m", SOCWRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 3603.91; atol = 1e0) end @testset "5-bus with pwl costs" begin - result = run_opf("../test/data/matpower/case5_pwlc.m", SOCWRPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case5_pwlc.m", SOCWRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 42895; atol = 1e0) end @testset "5-bus with dangling bus" begin - result = run_opf("../test/data/matpower/case5_db.m", SOCWRPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case5_db.m", SOCWRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 16739.1; atol = 1e0) end @testset "6-bus case" begin - result = run_opf("../test/data/matpower/case6.m", SOCWRPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case6.m", SOCWRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 11472.3; atol = 1e0) end @testset "24-bus rts case" begin - result = run_opf("../test/data/matpower/case24.m", SOCWRPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case24.m", SOCWRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 70690.7; atol = 1e0) @@ -587,56 +587,56 @@ end @testset "test soc conic form opf" begin @testset "3-bus case" begin - result = run_opf("../test/data/matpower/case3.m", SOCWRConicPowerModel, scs_solver) + result = run_opf("../test/data/matpower/case3.m", SOCWRConicPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 5736.94; atol = 2e0) end @testset "5-bus transformer swap case" begin - result = run_opf("../test/data/matpower/case5.m", SOCWRConicPowerModel, scs_solver) + result = run_opf("../test/data/matpower/case5.m", SOCWRConicPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 15051.4; atol = 1e1) end @testset "5-bus asymmetric case" begin - result = run_opf("../test/data/matpower/case5_asym.m", SOCWRConicPowerModel, scs_solver) + result = run_opf("../test/data/matpower/case5_asym.m", SOCWRConicPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 14999.7; atol = 1e0) end @testset "5-bus gap case" begin - result = run_opf("../test/data/matpower/case5_gap.m", SOCWRConicPowerModel, scs_solver) + result = run_opf("../test/data/matpower/case5_gap.m", SOCWRConicPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], -28236.8; atol = 1e0) end @testset "5-bus with asymmetric line charge" begin - result = run_opf("../test/data/pti/case5_alc.raw", SOCWRConicPowerModel, scs_solver) + result = run_opf("../test/data/pti/case5_alc.raw", SOCWRConicPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 1005.27; atol = 1e0) end @testset "5-bus with negative generators" begin - result = run_opf("../test/data/matpower/case5_npg.m", SOCWRConicPowerModel, scs_solver) + result = run_opf("../test/data/matpower/case5_npg.m", SOCWRConicPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 3551.71; atol = 40) end @testset "5-bus with pwl costs" begin - result = run_opf("../test/data/matpower/case5_pwlc.m", SOCWRConicPowerModel, scs_solver) + result = run_opf("../test/data/matpower/case5_pwlc.m", SOCWRConicPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 42889; atol = 1e0) end @testset "6-bus case" begin - result = run_opf("../test/data/matpower/case6.m", SOCWRConicPowerModel, scs_solver) + result = run_opf("../test/data/matpower/case6.m", SOCWRConicPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL #@test isapprox(result["objective"], 11472.2; atol = 3e0) @test isapprox(result["objective"], 11451.5; atol = 3e0) end @testset "24-bus rts case" begin - result = run_opf("../test/data/matpower/case24.m", SOCWRConicPowerModel, scs_solver) + result = run_opf("../test/data/matpower/case24.m", SOCWRConicPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL #@test isapprox(result["objective"], 70693.9; atol = 1e0) @@ -650,49 +650,49 @@ end @testset "test soc distflow opf_bf" begin @testset "3-bus case" begin - result = run_opf_bf("../test/data/matpower/case3.m", SOCBFPowerModel, ipopt_solver) + result = run_opf_bf("../test/data/matpower/case3.m", SOCBFPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 5746.7; atol = 1e0) end @testset "5-bus transformer swap case" begin - result = run_opf_bf("../test/data/matpower/case5.m", SOCBFPowerModel, ipopt_solver) + result = run_opf_bf("../test/data/matpower/case5.m", SOCBFPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 15051; atol = 1e0) end @testset "5-bus asymmetric case" begin - result = run_opf_bf("../test/data/matpower/case5_asym.m", SOCBFPowerModel, ipopt_solver) + result = run_opf_bf("../test/data/matpower/case5_asym.m", SOCBFPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 14999; atol = 1e0) end @testset "5-bus gap case" begin - result = run_opf_bf("../test/data/matpower/case5_gap.m", SOCBFPowerModel, ipopt_solver) + result = run_opf_bf("../test/data/matpower/case5_gap.m", SOCBFPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], -27659.8; atol = 1e0) end @testset "5-bus with asymmetric line charge" begin - result = run_opf_bf("../test/data/pti/case5_alc.raw", SOCBFPowerModel, ipopt_solver) + result = run_opf_bf("../test/data/pti/case5_alc.raw", SOCBFPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 1005.27; atol = 1e0) end @testset "5-bus with pwl costs" begin - result = run_opf_bf("../test/data/matpower/case5_pwlc.m", SOCBFPowerModel, ipopt_solver) + result = run_opf_bf("../test/data/matpower/case5_pwlc.m", SOCBFPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 42895; atol = 1e0) end @testset "6-bus case" begin - result = run_opf_bf("../test/data/matpower/case6.m", SOCBFPowerModel, ipopt_solver) + result = run_opf_bf("../test/data/matpower/case6.m", SOCBFPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 11472.3; atol = 1e0) end @testset "24-bus rts case" begin - result = run_opf_bf("../test/data/matpower/case24.m", SOCBFPowerModel, ipopt_solver) + result = run_opf_bf("../test/data/matpower/case24.m", SOCBFPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 70690.7; atol = 1e0) @@ -700,7 +700,7 @@ end @testset "3-bus case w/ r,x=0 on branch" begin mp_data = PowerModels.parse_file("../test/data/matpower/case3.m") mp_data["branch"]["1"]["br_r"] = mp_data["branch"]["1"]["br_x"] = 0.0 - result = run_opf_bf(mp_data, SOCBFPowerModel, ipopt_solver) + result = run_opf_bf(mp_data, SOCBFPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED end @testset "14-bus variable bounds" begin @@ -711,25 +711,25 @@ end @testset "test soc conic distflow opf_bf" begin @testset "3-bus case" begin - result = run_opf_bf("../test/data/matpower/case3.m", SOCBFConicPowerModel, scs_solver) + result = run_opf_bf("../test/data/matpower/case3.m", SOCBFConicPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 5746.7; atol = 5e1) end @testset "5-bus transformer swap case" begin - result = run_opf_bf("../test/data/matpower/case5.m", SOCBFConicPowerModel, scs_solver) + result = run_opf_bf("../test/data/matpower/case5.m", SOCBFConicPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 15051; atol = 1e1) end @testset "5-bus asymmetric case" begin - result = run_opf_bf("../test/data/matpower/case5_asym.m", SOCBFConicPowerModel, scs_solver) + result = run_opf_bf("../test/data/matpower/case5_asym.m", SOCBFConicPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 14999; atol = 1e1) end @testset "5-bus with negative generators" begin - result = run_opf_bf("../test/data/matpower/case5_npg.m", SOCBFConicPowerModel, scs_solver) + result = run_opf_bf("../test/data/matpower/case5_npg.m", SOCBFConicPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 3593.0; atol = 1e1) @@ -742,14 +742,14 @@ end @testset "test linear distflow opf_bf" begin @testset "3-bus case" begin - result = run_opf_bf("../test/data/matpower/case3.m", BFAPowerModel, ipopt_solver) + result = run_opf_bf("../test/data/matpower/case3.m", BFAPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 5638.97; atol = 1e0) end @testset "5-bus transformer swap case" begin data = PowerModels.parse_file("../test/data/matpower/case5.m") - result = run_opf_bf(data, BFAPowerModel, ipopt_solver) + result = run_opf_bf(data, BFAPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 14810; atol = 1e0) @@ -758,37 +758,37 @@ end atol = 1e-3) end @testset "5-bus asymmetric case" begin - result = run_opf_bf("../test/data/matpower/case5_asym.m", BFAPowerModel, ipopt_solver) + result = run_opf_bf("../test/data/matpower/case5_asym.m", BFAPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 14810; atol = 1e0) end @testset "5-bus gap case" begin - result = run_opf_bf("../test/data/matpower/case5_gap.m", BFAPowerModel, ipopt_solver) + result = run_opf_bf("../test/data/matpower/case5_gap.m", BFAPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], -27410.0; atol = 1e0) end @testset "5-bus with asymmetric line charge" begin - result = run_opf_bf("../test/data/pti/case5_alc.raw", BFAPowerModel, ipopt_solver) + result = run_opf_bf("../test/data/pti/case5_alc.raw", BFAPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 1002.46; atol = 1e0) end @testset "5-bus with pwl costs" begin - result = run_opf_bf("../test/data/matpower/case5_pwlc.m", BFAPowerModel, ipopt_solver) + result = run_opf_bf("../test/data/matpower/case5_pwlc.m", BFAPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 42565.8; atol = 1e0) end @testset "6-bus case" begin - result = run_opf_bf("../test/data/matpower/case6.m", BFAPowerModel, ipopt_solver) + result = run_opf_bf("../test/data/matpower/case6.m", BFAPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 11277.9; atol = 1e0) end # @testset "24-bus rts case" begin - # result = run_opf_bf("../test/data/matpower/case24.m", BFAPowerModel, ipopt_solver) + # result = run_opf_bf("../test/data/matpower/case24.m", BFAPowerModel, nlp_solver) # # @test result["termination_status"] == LOCALLY_SOLVED # @test isapprox(result["objective"], 70690.7; atol = 1e0) @@ -796,7 +796,7 @@ end @testset "3-bus case w/ r,x=0 on branch" begin mp_data = PowerModels.parse_file("../test/data/matpower/case3.m") mp_data["branch"]["1"]["br_r"] = mp_data["branch"]["1"]["br_x"] = 0.0 - result = run_opf_bf(mp_data, BFAPowerModel, ipopt_solver) + result = run_opf_bf(mp_data, BFAPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED end @testset "14-bus variable bounds" begin @@ -807,37 +807,37 @@ end @testset "test qc opf" begin @testset "3-bus case" begin - result = run_opf("../test/data/matpower/case3.m", QCRMPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case3.m", QCRMPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 5780; atol = 1e0) end @testset "5-bus asymmetric case" begin - result = run_opf("../test/data/matpower/case5_asym.m", QCRMPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case5_asym.m", QCRMPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 15921; atol = 1e0) end @testset "5-bus gap case" begin - result = run_opf("../test/data/matpower/case5_gap.m", QCRMPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case5_gap.m", QCRMPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], -27659.8; atol = 1e0) end @testset "5-bus with asymmetric line charge" begin - result = run_opf("../test/data/pti/case5_alc.raw", QCRMPowerModel, ipopt_solver) + result = run_opf("../test/data/pti/case5_alc.raw", QCRMPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 1005.27; atol = 1e0) end @testset "5-bus with pwl costs" begin - result = run_opf("../test/data/matpower/case5_pwlc.m", QCRMPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case5_pwlc.m", QCRMPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 42895; atol = 1e0) end @testset "6-bus case" begin - result = run_opf("../test/data/matpower/case6.m", QCRMPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case6.m", QCRMPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 11484.2; atol = 1e0) @@ -845,7 +845,7 @@ end @test isapprox(result["solution"]["bus"]["4"]["va"], 0.0; atol = 1e-4) end @testset "24-bus rts case" begin - result = run_opf("../test/data/matpower/case24.m", QCRMPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case24.m", QCRMPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 76599.9; atol = 1e0) @@ -858,25 +858,25 @@ end @testset "test qc opf with trilinear convexhull relaxation" begin @testset "3-bus case" begin - result = run_opf("../test/data/matpower/case3.m", QCLSPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case3.m", QCLSPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 5817.91; atol = 1e0) end @testset "5-bus asymmetric case" begin - result = run_opf("../test/data/matpower/case5_asym.m", QCLSPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case5_asym.m", QCLSPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 15929.2; atol = 1e0) end @testset "5-bus gap case" begin - result = run_opf("../test/data/matpower/case5_gap.m", QCLSPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case5_gap.m", QCLSPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], -27659.8; atol = 1e0) end @testset "6-bus case" begin - result = run_opf("../test/data/matpower/case6.m", QCLSPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case6.m", QCLSPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 11512.9; atol = 1e0) @@ -884,7 +884,7 @@ end @test isapprox(result["solution"]["bus"]["4"]["va"], 0.0; atol = 1e-4) end @testset "24-bus rts case" begin - result = run_opf("../test/data/matpower/case24.m", QCLSPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case24.m", QCLSPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 76785.4; atol = 1e0) @@ -898,31 +898,31 @@ end @testset "test sdp opf" begin @testset "3-bus case" begin - result = run_opf("../test/data/matpower/case3.m", SDPWRMPowerModel, scs_solver) + result = run_opf("../test/data/matpower/case3.m", SDPWRMPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 5818.00; atol = 1e1) end @testset "5-bus asymmetric case" begin - result = run_opf("../test/data/matpower/case5_asym.m", SDPWRMPowerModel, scs_solver) + result = run_opf("../test/data/matpower/case5_asym.m", SDPWRMPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 16662.0; atol = 1e1) end @testset "5-bus gap case" begin - result = run_opf("../test/data/matpower/case5_gap.m", SDPWRMPowerModel, scs_solver) + result = run_opf("../test/data/matpower/case5_gap.m", SDPWRMPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], -28236.1; atol = 1e1) end @testset "5-bus with asymmetric line charge" begin - result = run_opf("../test/data/pti/case5_alc.raw", SDPWRMPowerModel, scs_solver) + result = run_opf("../test/data/pti/case5_alc.raw", SDPWRMPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 1005.31; atol = 1e0) end @testset "5-bus with negative generators" begin - result = run_opf("../test/data/matpower/case5_npg.m", SDPWRMPowerModel, scs_solver) + result = run_opf("../test/data/matpower/case5_npg.m", SDPWRMPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL #@test isapprox(result["objective"], 6827.34; atol = 1e0) @@ -930,13 +930,13 @@ end end # too slow for unit tests # @testset "14-bus case" begin - # result = run_opf("../test/data/matpower/case14.m", SDPWRMPowerModel, scs_solver) + # result = run_opf("../test/data/matpower/case14.m", SDPWRMPowerModel, sdp_solver) # @test result["termination_status"] == OPTIMAL # @test isapprox(result["objective"], 8081.52; atol = 1e0) # end @testset "6-bus case" begin - result = run_opf("../test/data/matpower/case6.m", SDPWRMPowerModel, scs_solver) + result = run_opf("../test/data/matpower/case6.m", SDPWRMPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL #@test isapprox(result["objective"], 11580.8; atol = 1e1) @@ -951,33 +951,33 @@ end @testset "test sdp opf with constraint decomposition" begin @testset "3-bus case" begin - result = run_opf("../test/data/matpower/case3.m", SparseSDPWRMPowerModel, scs_solver) + result = run_opf("../test/data/matpower/case3.m", SparseSDPWRMPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL #@test isapprox(result["objective"], 5851.23; atol = 1e1) @test isapprox(result["objective"], 5818.00; atol = 1e1) end @testset "5-bus with asymmetric line charge" begin - result = run_opf("../test/data/pti/case5_alc.raw", SparseSDPWRMPowerModel, scs_solver) + result = run_opf("../test/data/pti/case5_alc.raw", SparseSDPWRMPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 1005.31; atol = 1e0) end @testset "9-bus cholesky PosDefException" begin - result = run_opf("../test/data/matpower/case9.m", SparseSDPWRMPowerModel, scs_solver) + result = run_opf("../test/data/matpower/case9.m", SparseSDPWRMPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 347.746; atol = 1e0) end @testset "14-bus case" begin - result = run_opf("../test/data/matpower/case14.m", SparseSDPWRMPowerModel, scs_solver) + result = run_opf("../test/data/matpower/case14.m", SparseSDPWRMPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 8081.5; atol = 1e0) end # multiple components are not currently supported by this form # @testset "6-bus case" begin - # result = run_opf("../test/data/matpower/case6.m", SparseSDPWRMPowerModel, scs_solver) + # result = run_opf("../test/data/matpower/case6.m", SparseSDPWRMPowerModel, sdp_solver) # @test result["termination_status"] == OPTIMAL # @test isapprox(result["objective"], 11578.8; atol = 1e0) @@ -1004,7 +1004,7 @@ end pm.ext[:SDconstraintDecomposition] = PowerModels._SDconstraintDecomposition(groups, lookup_index, sigma) PowerModels.build_opf(pm) - result = optimize_model!(pm, optimizer=scs_solver) + result = optimize_model!(pm, optimizer=sdp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 1005.31; atol = 1e0) diff --git a/test/ots.jl b/test/ots.jl index e68e078f0..1a2b5fdf9 100644 --- a/test/ots.jl +++ b/test/ots.jl @@ -13,7 +13,7 @@ end @testset "test ac ots" begin @testset "3-bus case" begin - result = run_ots("../test/data/matpower/case3.m", ACPPowerModel, juniper_solver) + result = run_ots("../test/data/matpower/case3.m", ACPPowerModel, minlp_solver) check_br_status(result["solution"], 3, 3) @@ -24,7 +24,7 @@ end #= # remove due to linux stability issue @testset "5-bus case" begin - result = run_ots("../test/data/matpower/case5.m", ACPPowerModel, juniper_solver) + result = run_ots("../test/data/matpower/case5.m", ACPPowerModel, minlp_solver) check_br_status(result["solution"], 0, 0) @@ -35,7 +35,7 @@ end end =# @testset "5-bus with negative branch reactance case" begin - result = run_ots("../test/data/matpower/case5_ext.m", ACPPowerModel, juniper_solver) + result = run_ots("../test/data/matpower/case5_ext.m", ACPPowerModel, minlp_solver) check_br_status(result["solution"], 5, 5) @@ -45,7 +45,7 @@ end #@test result["objective"] < 16600 end @testset "5-bus with asymmetric line charge" begin - result = run_ots("../test/data/pti/case5_alc.raw", ACPPowerModel, juniper_solver) + result = run_ots("../test/data/pti/case5_alc.raw", ACPPowerModel, minlp_solver) check_br_status(result["solution"], 4, 4) @@ -54,7 +54,7 @@ end end #Omitting this test, returns local infeasible #@testset "6-bus case" begin - # result = run_ots("../test/data/matpower/case6.m", ACPPowerModel, juniper_solver) + # result = run_ots("../test/data/matpower/case6.m", ACPPowerModel, minlp_solver) # check_br_status(result["solution"], 0, 0) @@ -67,7 +67,7 @@ end @testset "test dc ots" begin @testset "3-bus case" begin - result = run_ots("../test/data/matpower/case3.m", DCPPowerModel, juniper_solver) + result = run_ots("../test/data/matpower/case3.m", DCPPowerModel, minlp_solver) check_br_status(result["solution"], 3, 3) @@ -75,7 +75,7 @@ end @test isapprox(result["objective"], 5782.0; atol = 1e0) end @testset "5-bus case" begin - result = run_ots("../test/data/matpower/case5.m", DCPPowerModel, juniper_solver) + result = run_ots("../test/data/matpower/case5.m", DCPPowerModel, minlp_solver) check_br_status(result["solution"], 5, 5) @@ -83,7 +83,7 @@ end @test isapprox(result["objective"], 14991.2; atol = 1e0) end @testset "5-bus case, MIP solver" begin - result = run_ots("../test/data/matpower/case5.m", DCPPowerModel, cbc_solver) + result = run_ots("../test/data/matpower/case5.m", DCPPowerModel, milp_solver) check_br_status(result["solution"], 5, 5) @@ -91,7 +91,7 @@ end @test isapprox(result["objective"], 14991.3; atol = 1e0) end @testset "5-bus with negative branch reactance case" begin - result = run_ots("../test/data/matpower/case5_ext.m", DCPPowerModel, juniper_solver) + result = run_ots("../test/data/matpower/case5_ext.m", DCPPowerModel, minlp_solver) check_br_status(result["solution"], 5, 5) @@ -99,7 +99,7 @@ end @test isapprox(result["objective"], 14991.2; atol = 1e0) end @testset "6-bus case" begin - result = run_ots("../test/data/matpower/case6.m", DCPPowerModel, juniper_solver) + result = run_ots("../test/data/matpower/case6.m", DCPPowerModel, minlp_solver) check_br_status(result["solution"], 6, 6) @@ -111,7 +111,7 @@ end @testset "test dc+ll ots" begin @testset "3-bus case" begin - result = run_ots("../test/data/matpower/case3.m", DCPLLPowerModel, juniper_solver) + result = run_ots("../test/data/matpower/case3.m", DCPLLPowerModel, minlp_solver) check_br_status(result["solution"], 3, 3) @@ -119,7 +119,7 @@ end @test isapprox(result["objective"], 5885.2; atol = 1e0) end @testset "5-bus case" begin - result = run_ots("../test/data/matpower/case5.m", DCPLLPowerModel, juniper_solver) + result = run_ots("../test/data/matpower/case5.m", DCPLLPowerModel, minlp_solver) check_br_status(result["solution"], 5, 5) @@ -127,7 +127,7 @@ end @test isapprox(result["objective"], 15275.2; atol = 1e0) end @testset "5-bus with negative branch reactance case" begin - result = run_ots("../test/data/matpower/case5_ext.m", DCPLLPowerModel, juniper_solver) + result = run_ots("../test/data/matpower/case5_ext.m", DCPLLPowerModel, minlp_solver) check_br_status(result["solution"], 5, 5) @@ -135,7 +135,7 @@ end @test isapprox(result["objective"], 15275.2; atol = 1e0) end @testset "6-bus case" begin - result = run_ots("../test/data/matpower/case6.m", DCPLLPowerModel, juniper_solver) + result = run_ots("../test/data/matpower/case6.m", DCPLLPowerModel, minlp_solver) check_br_status(result["solution"], 6, 6) @@ -147,7 +147,7 @@ end @testset "test soc ots" begin @testset "3-bus case" begin - result = run_ots("../test/data/matpower/case3.m", SOCWRPowerModel, juniper_solver) + result = run_ots("../test/data/matpower/case3.m", SOCWRPowerModel, minlp_solver) check_br_status(result["solution"], 3, 3) @@ -155,7 +155,7 @@ end @test isapprox(result["objective"], 5746.7; atol = 1e0) end @testset "5-bus case" begin - result = run_ots("../test/data/matpower/case5.m", SOCWRPowerModel, juniper_solver) + result = run_ots("../test/data/matpower/case5.m", SOCWRPowerModel, minlp_solver) check_br_status(result["solution"], 5, 5) @@ -163,7 +163,7 @@ end @test isapprox(result["objective"], 15051.4; atol = 5e1) end @testset "5-bus with negative branch reactance case" begin - result = run_ots("../test/data/matpower/case5_ext.m", SOCWRPowerModel, juniper_solver) + result = run_ots("../test/data/matpower/case5_ext.m", SOCWRPowerModel, minlp_solver) check_br_status(result["solution"], 5, 5) @@ -171,7 +171,7 @@ end @test isapprox(result["objective"], 15009.9; atol = 1e0) end @testset "5-bus with asymmetric line charge" begin - result = run_ots("../test/data/pti/case5_alc.raw", SOCWRPowerModel, juniper_solver) + result = run_ots("../test/data/pti/case5_alc.raw", SOCWRPowerModel, minlp_solver) check_br_status(result["solution"], 4, 4) @@ -179,7 +179,7 @@ end @test isapprox(result["objective"], 1004.8; atol = 5e0) end @testset "6-bus case" begin - result = run_ots("../test/data/matpower/case6.m", SOCWRPowerModel, juniper_solver) + result = run_ots("../test/data/matpower/case6.m", SOCWRPowerModel, minlp_solver) check_br_status(result["solution"], 6, 6) @@ -191,7 +191,7 @@ end @testset "test qc ots" begin @testset "3-bus case" begin - result = run_ots("../test/data/matpower/case3.m", QCRMPowerModel, juniper_solver) + result = run_ots("../test/data/matpower/case3.m", QCRMPowerModel, minlp_solver) check_br_status(result["solution"], 3, 3) @@ -199,7 +199,7 @@ end @test isapprox(result["objective"], 5748.0; atol = 1e0) end @testset "5-bus case" begin - result = run_ots("../test/data/matpower/case5.m", QCRMPowerModel, juniper_solver) + result = run_ots("../test/data/matpower/case5.m", QCRMPowerModel, minlp_solver) check_br_status(result["solution"], 5, 5) @@ -207,7 +207,7 @@ end @test isapprox(result["objective"], 15051.4; atol = 5e1) end @testset "5-bus with negative branch reactance case" begin - result = run_ots("../test/data/matpower/case5_ext.m", QCRMPowerModel, juniper_solver) + result = run_ots("../test/data/matpower/case5_ext.m", QCRMPowerModel, minlp_solver) check_br_status(result["solution"], 5, 5) @@ -215,7 +215,7 @@ end @test isapprox(result["objective"], 15010.0; atol = 1e0) end @testset "5-bus asymmetric case" begin - result = run_ots("../test/data/matpower/case5_asym.m", QCRMPowerModel, juniper_solver) + result = run_ots("../test/data/matpower/case5_asym.m", QCRMPowerModel, minlp_solver) check_br_status(result["solution"], 6, 6) @@ -223,7 +223,7 @@ end @test isapprox(result["objective"], 15036.9; atol = 1e0) end @testset "5-bus with asymmetric line charge" begin - result = run_ots("../test/data/pti/case5_alc.raw", QCRMPowerModel, juniper_solver) + result = run_ots("../test/data/pti/case5_alc.raw", QCRMPowerModel, minlp_solver) # updated ub to 5 on 04/21/2021 to fix cross platform stability check_br_status(result["solution"], 4, 5) @@ -232,7 +232,7 @@ end @test isapprox(result["objective"], 1003.97; atol = 5e0) end @testset "6-bus case" begin - result = run_ots("../test/data/matpower/case6.m", QCRMPowerModel, juniper_solver) + result = run_ots("../test/data/matpower/case6.m", QCRMPowerModel, minlp_solver) check_br_status(result["solution"], 6, 6) @@ -244,7 +244,7 @@ end @testset "test lpac ots" begin @testset "3-bus case" begin - result = run_ots("../test/data/matpower/case3.m", LPACCPowerModel, juniper_solver) + result = run_ots("../test/data/matpower/case3.m", LPACCPowerModel, minlp_solver) check_br_status(result["solution"], 3, 3) @@ -252,7 +252,7 @@ end @test isapprox(result["objective"], 5908.98; atol = 1e0) end @testset "5-bus case" begin - result = run_ots("../test/data/matpower/case5.m", LPACCPowerModel, juniper_solver) + result = run_ots("../test/data/matpower/case5.m", LPACCPowerModel, minlp_solver) check_br_status(result["solution"], 5, 5) @@ -260,7 +260,7 @@ end @test isapprox(result["objective"], 15241.4; atol = 5e1) end @testset "5-bus with negative branch reactance case" begin - result = run_ots("../test/data/matpower/case5_ext.m", LPACCPowerModel, juniper_solver) + result = run_ots("../test/data/matpower/case5_ext.m", LPACCPowerModel, minlp_solver) check_br_status(result["solution"], 5, 5) @@ -268,7 +268,7 @@ end @test isapprox(result["objective"], 15241.4; atol = 1e0) end @testset "5-bus asymmetric case" begin - result = run_ots("../test/data/matpower/case5_asym.m", LPACCPowerModel, juniper_solver) + result = run_ots("../test/data/matpower/case5_asym.m", LPACCPowerModel, minlp_solver) check_br_status(result["solution"], 5, 5) @@ -276,7 +276,7 @@ end @test isapprox(result["objective"], 15246.9; atol = 1e0) end @testset "5-bus with asymmetric line charge" begin - result = run_ots("../test/data/pti/case5_alc.raw", LPACCPowerModel, juniper_solver) + result = run_ots("../test/data/pti/case5_alc.raw", LPACCPowerModel, minlp_solver) check_br_status(result["solution"], 4, 4) @@ -284,7 +284,7 @@ end @test isapprox(result["objective"], 998.4; atol = 5e0) end @testset "6-bus case" begin - result = run_ots("../test/data/matpower/case6.m", LPACCPowerModel, juniper_solver) + result = run_ots("../test/data/matpower/case6.m", LPACCPowerModel, minlp_solver) check_br_status(result["solution"], 6, 6) diff --git a/test/output.jl b/test/output.jl index 14be816c5..b396c10ec 100644 --- a/test/output.jl +++ b/test/output.jl @@ -1,7 +1,7 @@ @testset "test output api" begin @testset "24-bus rts case" begin - result = run_opf("../test/data/matpower/case24.m", ACPPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case24.m", ACPPowerModel, nlp_solver) @test haskey(result, "optimizer") == true @test haskey(result, "termination_status") == true @@ -29,7 +29,7 @@ gen["cost"] = gen["cost"][length(gen["cost"])-1:end] end end - result = run_opf(data, DCPPowerModel, cbc_solver) + result = run_opf(data, DCPPowerModel, milp_solver) @test haskey(result, "optimizer") @test haskey(result, "termination_status") @@ -44,7 +44,7 @@ end @testset "test branch flow output" begin @testset "24-bus rts case ac opf" begin - result = run_opf("../test/data/matpower/case24.m", ACPPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case24.m", ACPPowerModel, nlp_solver) @test haskey(result, "optimizer") == true @test haskey(result, "termination_status") == true @@ -70,7 +70,7 @@ end # A DCPPowerModel test is important because it does have variables for the reverse side of the branchs @testset "3-bus case dc opf" begin - result = run_opf("../test/data/matpower/case3.m", DCPPowerModel, ipopt_solver) + result = run_opf("../test/data/matpower/case3.m", DCPPowerModel, nlp_solver) @test haskey(result, "solution") == true @test haskey(result["solution"], "branch") == true @@ -89,7 +89,7 @@ end end @testset "24-bus rts case ac pf" begin - result = run_pf("../test/data/matpower/case24.m", ACPPowerModel, ipopt_solver) + result = run_pf("../test/data/matpower/case24.m", ACPPowerModel, nlp_solver) @test haskey(result, "optimizer") == true @test haskey(result, "termination_status") == true @@ -115,7 +115,7 @@ end # A DCPPowerModel test is important because it does have variables for the reverse side of the branchs @testset "3-bus case dc pf" begin - result = run_pf("../test/data/matpower/case3.m", DCPPowerModel, ipopt_solver) + result = run_pf("../test/data/matpower/case3.m", DCPPowerModel, nlp_solver) @test haskey(result, "solution") == true @test haskey(result["solution"], "branch") == false @@ -138,7 +138,7 @@ end settings = Dict("output" => Dict("duals" => true)) data = PowerModels.parse_file("../test/data/matpower/case14.m") calc_thermal_limits!(data) - result = run_dc_opf(data, ipopt_solver, setting = settings) + result = run_dc_opf(data, nlp_solver, setting = settings) PowerModels.make_mixed_units!(result["solution"]) @testset "14 bus - kcl duals" begin @@ -158,7 +158,7 @@ end end - result = run_dc_opf("../test/data/matpower/case5.m", ipopt_solver, setting = settings) + result = run_dc_opf("../test/data/matpower/case5.m", nlp_solver, setting = settings) PowerModels.make_mixed_units!(result["solution"]) @testset "5 bus - kcl duals" begin @@ -180,7 +180,7 @@ end end - result = run_opf("../test/data/matpower/case5.m", SOCWRPowerModel, ipopt_solver, setting = settings) + result = run_opf("../test/data/matpower/case5.m", SOCWRPowerModel, nlp_solver, setting = settings) @testset "5 bus - kcl duals soc qp" begin for (i, bus) in result["solution"]["bus"] @test bus["lam_kcl_r"] <= -2900.00 @@ -190,7 +190,7 @@ end end end - result = run_opf("../test/data/matpower/case5.m", SOCWRConicPowerModel, scs_solver, setting = settings) + result = run_opf("../test/data/matpower/case5.m", SOCWRConicPowerModel, sdp_solver, setting = settings) @testset "5 bus - kcl duals soc conic" begin for (i, bus) in result["solution"]["bus"] @test bus["lam_kcl_r"] <= -2900.00 @@ -200,7 +200,7 @@ end end end - result = run_opf("../test/data/matpower/case5.m", ACPPowerModel, ipopt_solver, setting = settings) + result = run_opf("../test/data/matpower/case5.m", ACPPowerModel, nlp_solver, setting = settings) @testset "5 bus - kcl duals acp" begin for (i, bus) in result["solution"]["bus"] @test bus["lam_kcl_r"] <= -1000.00 @@ -210,7 +210,7 @@ end end end - result = run_opf("../test/data/matpower/case5.m", ACRPowerModel, ipopt_solver, setting = settings) + result = run_opf("../test/data/matpower/case5.m", ACRPowerModel, nlp_solver, setting = settings) @testset "5 bus - kcl duals acr" begin for (i, bus) in result["solution"]["bus"] @test bus["lam_kcl_r"] <= -1000.00 @@ -220,7 +220,7 @@ end end end - result = run_opf("../test/data/matpower/case5.m", ACTPowerModel, ipopt_solver, setting = settings) + result = run_opf("../test/data/matpower/case5.m", ACTPowerModel, nlp_solver, setting = settings) @testset "5 bus - kcl duals act" begin for (i, bus) in result["solution"]["bus"] @test bus["lam_kcl_r"] <= -1000.00 @@ -237,7 +237,7 @@ end # test case where generator status is 1 but the gen_bus status is 0 data = parse_file("../test/data/matpower/case5.m") data["bus"]["4"]["bus_type"] = 4 - result = run_ac_opf(data, ipopt_solver) + result = run_ac_opf(data, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 10128.6; atol = 1e0) @@ -246,7 +246,7 @@ end @testset "test solution processors" begin @testset "sol_vr_to_vp" begin - result = run_opf("../test/data/matpower/case5.m", ACRPowerModel, ipopt_solver, solution_processors=[sol_data_model!]) + result = run_opf("../test/data/matpower/case5.m", ACRPowerModel, nlp_solver, solution_processors=[sol_data_model!]) for (i,bus) in result["solution"]["bus"] @test haskey(bus, "vm") && haskey(bus, "va") @@ -254,7 +254,7 @@ end end @testset "sol_w_to_vm" begin - result = run_opf("../test/data/matpower/case5.m", SOCWRPowerModel, ipopt_solver, solution_processors=[sol_data_model!]) + result = run_opf("../test/data/matpower/case5.m", SOCWRPowerModel, nlp_solver, solution_processors=[sol_data_model!]) for (i,bus) in result["solution"]["bus"] @test haskey(bus, "vm") @@ -263,7 +263,7 @@ end end @testset "sol_phi_to_vm" begin - result = run_opf("../test/data/matpower/case5.m", LPACCPowerModel, ipopt_solver, solution_processors=[sol_data_model!]) + result = run_opf("../test/data/matpower/case5.m", LPACCPowerModel, nlp_solver, solution_processors=[sol_data_model!]) for (i,bus) in result["solution"]["bus"] @test haskey(bus, "vm") @@ -278,13 +278,13 @@ end function solution_feedback(case, ac_opf_obj) data = PowerModels.parse_file(case) - opf_result = run_ac_opf(data, ipopt_solver) + opf_result = run_ac_opf(data, nlp_solver) @test opf_result["termination_status"] == LOCALLY_SOLVED @test isapprox(opf_result["objective"], ac_opf_obj; atol = 1e0) PowerModels.update_data!(data, opf_result["solution"]) - pf_result = run_ac_pf(data, ipopt_solver) + pf_result = run_ac_pf(data, nlp_solver) @test pf_result["termination_status"] == LOCALLY_SOLVED @test isapprox(pf_result["objective"], 0.0; atol = 1e-3) diff --git a/test/pf-native.jl b/test/pf-native.jl index 0c56a5a03..d81758a2a 100644 --- a/test/pf-native.jl +++ b/test/pf-native.jl @@ -2,7 +2,7 @@ # degenerate due to no slack bus # @testset "3-bus case" begin # data = PowerModels.parse_file("../test/data/matpower/case3.m") - # result = run_dc_pf(data, ipopt_solver) + # result = run_dc_pf(data, nlp_solver) # native = compute_dc_pf(data) # for (i,bus) in data["bus"] @@ -13,7 +13,7 @@ # end @testset "5-bus case" begin data = PowerModels.parse_file("../test/data/matpower/case5.m") - result = run_dc_pf(data, ipopt_solver) + result = run_dc_pf(data, nlp_solver) native = compute_dc_pf(data) @test length(native) >= 5 @@ -30,7 +30,7 @@ end @testset "5-bus asymmetric case" begin data = PowerModels.parse_file("../test/data/matpower/case5_asym.m") - result = run_dc_pf(data, ipopt_solver) + result = run_dc_pf(data, nlp_solver) native = compute_dc_pf(data) for (i,bus) in data["bus"] @@ -41,7 +41,7 @@ end @testset "5-bus multiple slack gens case" begin data = PowerModels.parse_file("../test/data/matpower/case5_ext.m") - result = run_dc_pf(data, ipopt_solver) + result = run_dc_pf(data, nlp_solver) native = compute_dc_pf(data) for (i,bus) in data["bus"] @@ -55,7 +55,7 @@ # compute_dc_pf does not yet support multiple slack buses # @testset "6-bus case" begin # data = PowerModels.parse_file("../test/data/matpower/case6.m") - # result = run_dc_pf(data, ipopt_solver) + # result = run_dc_pf(data, nlp_solver) # native = compute_dc_pf(data) # for (i,bus) in data["bus"] @@ -66,7 +66,7 @@ # end @testset "24-bus rts case" begin data = PowerModels.parse_file("../test/data/matpower/case24.m") - result = run_dc_pf(data, ipopt_solver) + result = run_dc_pf(data, nlp_solver) native = compute_dc_pf(data) for (i,bus) in data["bus"] @@ -84,7 +84,7 @@ end # requires dc line support in ac solver # @testset "3-bus case" begin # data = PowerModels.parse_file("../test/data/matpower/case3.m") - # result = run_dc_pf(data, ipopt_solver) + # result = run_dc_pf(data, nlp_solver) # native = compute_dc_pf(data) # @test result["termination_status"] == LOCALLY_SOLVED @@ -105,7 +105,7 @@ end # end @testset "5-bus case" begin data = PowerModels.parse_file("../test/data/matpower/case5.m") - result = run_ac_pf(data, ipopt_solver) + result = run_ac_pf(data, nlp_solver) native = compute_ac_pf(data) @test result["termination_status"] == LOCALLY_SOLVED @@ -132,7 +132,7 @@ end end @testset "5-bus asymmetric case" begin data = PowerModels.parse_file("../test/data/matpower/case5_asym.m") - result = run_ac_pf(data, ipopt_solver) + result = run_ac_pf(data, nlp_solver) native = compute_ac_pf(data) @test result["termination_status"] == LOCALLY_SOLVED @@ -154,7 +154,7 @@ end end @testset "5-bus multiple slack gens case" begin data = PowerModels.parse_file("../test/data/matpower/case5_ext.m") - result = run_ac_pf(data, ipopt_solver) + result = run_ac_pf(data, nlp_solver) native = compute_ac_pf(data) @test result["termination_status"] == LOCALLY_SOLVED @@ -180,7 +180,7 @@ end # compute_ac_pf does not yet support multiple slack buses # @testset "6-bus case" begin # data = PowerModels.parse_file("../test/data/matpower/case6.m") - # result = run_ac_pf(data, ipopt_solver) + # result = run_ac_pf(data, nlp_solver) # native = compute_ac_pf(data) # @test result["termination_status"] == LOCALLY_SOLVED @@ -201,7 +201,7 @@ end # end @testset "14-bus case, vm fixed non-1.0 value" begin data = PowerModels.parse_file("../test/data/matpower/case14.m") - result = run_ac_pf(data, ipopt_solver) + result = run_ac_pf(data, nlp_solver) native = compute_ac_pf(data) @test result["termination_status"] == LOCALLY_SOLVED @@ -223,7 +223,7 @@ end end @testset "24-bus rts case" begin data = PowerModels.parse_file("../test/data/matpower/case24.m") - result = run_ac_pf(data, ipopt_solver) + result = run_ac_pf(data, nlp_solver) native = compute_ac_pf(data) @test result["termination_status"] == LOCALLY_SOLVED @@ -370,14 +370,14 @@ end # TODO extract number of iterations and test there is a reduction # Ipopt log can be used for manual verification, for now data = PowerModels.parse_file("../test/data/matpower/case24.m") - result = run_ac_pf(data, ipopt_solver) + result = run_ac_pf(data, nlp_solver) #result = run_ac_pf(data, JuMP.optimizer_with_attributes(Ipopt.Optimizer, "tol"=>1e-6)) @test result["termination_status"] == LOCALLY_SOLVED update_data!(data, result["solution"]) set_ac_pf_start_values!(data) - result_ws = run_ac_pf(data, ipopt_solver) + result_ws = run_ac_pf(data, nlp_solver) #result_ws = run_ac_pf(data, JuMP.optimizer_with_attributes(Ipopt.Optimizer, "tol"=>1e-6)) @test result_ws["termination_status"] == LOCALLY_SOLVED @@ -431,7 +431,7 @@ end @testset "test native ac pf solver options" begin @testset "5-bus case, finite_differencing" begin data = PowerModels.parse_file("../test/data/matpower/case5.m") - result = run_ac_pf(data, ipopt_solver) + result = run_ac_pf(data, nlp_solver) native = compute_ac_pf("../test/data/matpower/case5.m", finite_differencing=true) @test result["termination_status"] == LOCALLY_SOLVED @@ -453,7 +453,7 @@ end end @testset "5-bus case, flat_start" begin data = PowerModels.parse_file("../test/data/matpower/case5.m") - result = run_ac_pf(data, ipopt_solver) + result = run_ac_pf(data, nlp_solver) native = compute_ac_pf("../test/data/matpower/case5.m", flat_start=true) @test result["termination_status"] == LOCALLY_SOLVED diff --git a/test/pf.jl b/test/pf.jl index aafa2026b..2461d64fc 100644 --- a/test/pf.jl +++ b/test/pf.jl @@ -1,6 +1,6 @@ @testset "test ac polar pf" begin @testset "3-bus case" begin - result = run_ac_pf("../test/data/matpower/case3.m", ipopt_solver) + result = run_ac_pf("../test/data/matpower/case3.m", nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 0; atol = 1e-2) @@ -21,19 +21,19 @@ #@test isapprox(result["solution"]["dcline"]["1"]["qt"], 0.0647562; atol = 1e-5) end @testset "5-bus transformer swap case" begin - result = run_pf("../test/data/matpower/case5.m", ACPPowerModel, ipopt_solver) + result = run_pf("../test/data/matpower/case5.m", ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 0; atol = 1e-2) end @testset "5-bus asymmetric case" begin - result = run_pf("../test/data/matpower/case5_asym.m", ACPPowerModel, ipopt_solver) + result = run_pf("../test/data/matpower/case5_asym.m", ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 0; atol = 1e-2) end @testset "5-bus multiple slack gens case" begin - result = run_pf("../test/data/matpower/case5_ext.m", ACPPowerModel, ipopt_solver) + result = run_pf("../test/data/matpower/case5_ext.m", ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 0; atol = 1e-2) @@ -42,7 +42,7 @@ @test isapprox(result["solution"]["gen"]["1"]["qg"], 0.30; atol = 1e-3) end @testset "5-bus case with hvdc line" begin - result = run_ac_pf("../test/data/matpower/case5_dc.m", ipopt_solver) + result = run_ac_pf("../test/data/matpower/case5_dc.m", nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 0; atol = 1e-2) @@ -61,7 +61,7 @@ @test isapprox(result["solution"]["dcline"]["1"]["pt"], -0.089; atol = 1e-5) end @testset "6-bus case" begin - result = run_pf("../test/data/matpower/case6.m", ACPPowerModel, ipopt_solver) + result = run_pf("../test/data/matpower/case6.m", ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 0; atol = 1e-2) @@ -71,7 +71,7 @@ @test isapprox(result["solution"]["bus"]["4"]["va"], 0.00000; atol = 1e-3) end @testset "24-bus rts case" begin - result = run_pf("../test/data/matpower/case24.m", ACPPowerModel, ipopt_solver) + result = run_pf("../test/data/matpower/case24.m", ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 0; atol = 1e-2) @@ -81,7 +81,7 @@ end @testset "test ac rect pf" begin @testset "5-bus asymmetric case" begin - result = run_pf("../test/data/matpower/case5_asym.m", ACRPowerModel, ipopt_solver) + result = run_pf("../test/data/matpower/case5_asym.m", ACRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 0; atol = 1e-2) @@ -89,7 +89,7 @@ end #= # numerical issues with ipopt, likely div. by zero issue in jacobian @testset "5-bus case with hvdc line" begin - result = run_pf("../test/data/matpower/case5_dc.m", ACRPowerModel, ipopt_solver, setting = Dict("output" => Dict("branch_flows" => true))) + result = run_pf("../test/data/matpower/case5_dc.m", ACRPowerModel, nlp_solver, setting = Dict("output" => Dict("branch_flows" => true))) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 0; atol = 1e-2) @@ -114,13 +114,13 @@ end @testset "test ac tan pf" begin # removed for cross platform compat (julia v1.6, linux) # @testset "5-bus asymmetric case" begin - # result = run_pf("../test/data/matpower/case5_asym.m", ACTPowerModel, ipopt_solver) + # result = run_pf("../test/data/matpower/case5_asym.m", ACTPowerModel, nlp_solver) # @test result["termination_status"] == LOCALLY_SOLVED # @test isapprox(result["objective"], 0; atol = 1e-2) # end @testset "5-bus case with hvdc line" begin - result = run_pf("../test/data/matpower/case5_dc.m", ACTPowerModel, ipopt_solver, solution_processors=[sol_data_model!]) + result = run_pf("../test/data/matpower/case5_dc.m", ACTPowerModel, nlp_solver, solution_processors=[sol_data_model!]) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 0; atol = 1e-2) @@ -144,7 +144,7 @@ end @testset "test iv pf" begin @testset "3-bus case" begin - result = run_pf_iv("../test/data/matpower/case3.m", IVRPowerModel, ipopt_solver, solution_processors=[sol_data_model!]) + result = run_pf_iv("../test/data/matpower/case3.m", IVRPowerModel, nlp_solver, solution_processors=[sol_data_model!]) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 0; atol = 1e-2) @@ -163,7 +163,7 @@ end # @test isapprox(result["solution"]["dcline"]["1"]["qt"], 0.0647562; atol = 1e-5) #no reason to expect this is unique end @testset "5-bus case with hvdc line" begin - result = run_pf_iv("../test/data/matpower/case5_dc.m", IVRPowerModel, ipopt_solver, solution_processors=[sol_data_model!]) + result = run_pf_iv("../test/data/matpower/case5_dc.m", IVRPowerModel, nlp_solver, solution_processors=[sol_data_model!]) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 0; atol = 1e-2) @@ -186,7 +186,7 @@ end @testset "test dc pf" begin @testset "3-bus case" begin - result = run_dc_pf("../test/data/matpower/case3.m", ipopt_solver) + result = run_dc_pf("../test/data/matpower/case3.m", nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 0; atol = 1e-2) @@ -198,13 +198,13 @@ end @test isapprox(result["solution"]["bus"]["3"]["va"], -0.28291891895; atol = 1e-5) end @testset "5-bus asymmetric case" begin - result = run_dc_pf("../test/data/matpower/case5_asym.m", ipopt_solver) + result = run_dc_pf("../test/data/matpower/case5_asym.m", nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 0; atol = 1e-2) end @testset "5-bus multiple slack gens case" begin - result = run_dc_pf("../test/data/matpower/case5_ext.m", ipopt_solver) + result = run_dc_pf("../test/data/matpower/case5_ext.m", nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 0; atol = 1e-2) @@ -212,7 +212,7 @@ end @test isapprox(result["solution"]["gen"]["1"]["pg"], 0.40; atol = 1e-3) end @testset "6-bus case" begin - result = run_dc_pf("../test/data/matpower/case6.m", ipopt_solver) + result = run_dc_pf("../test/data/matpower/case6.m", nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 0; atol = 1e-2) @@ -220,7 +220,7 @@ end @test isapprox(result["solution"]["bus"]["4"]["va"], 0.00000; atol = 1e-5) end @testset "24-bus rts case" begin - result = run_pf("../test/data/matpower/case24.m", DCPPowerModel, ipopt_solver) + result = run_pf("../test/data/matpower/case24.m", DCPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 0; atol = 1e-2) @@ -230,7 +230,7 @@ end @testset "test matpower dc pf" begin @testset "5-bus case with matpower DCMP model" begin - result = run_pf("../test/data/matpower/case5.m", DCMPPowerModel, ipopt_solver) + result = run_pf("../test/data/matpower/case5.m", DCMPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @@ -245,7 +245,7 @@ end @testset "test soc pf" begin # started failing 05/22/2020 when ipopt moved to jll artifacts # @testset "3-bus case" begin - # result = run_pf("../test/data/matpower/case3.m", SOCWRPowerModel, ipopt_solver, solution_processors=[sol_data_model!]) + # result = run_pf("../test/data/matpower/case3.m", SOCWRPowerModel, nlp_solver, solution_processors=[sol_data_model!]) # @test result["termination_status"] == LOCALLY_SOLVED # @test isapprox(result["objective"], 0; atol = 1e-2) @@ -264,13 +264,13 @@ end # end # started failing 05/22/2020 when ipopt moved to jll artifacts (only on travis) # @testset "5-bus asymmetric case" begin - # result = run_pf("../test/data/matpower/case5_asym.m", SOCWRPowerModel, ipopt_solver) + # result = run_pf("../test/data/matpower/case5_asym.m", SOCWRPowerModel, nlp_solver) # @test result["termination_status"] == LOCALLY_SOLVED # @test isapprox(result["objective"], 0; atol = 1e-2) # end @testset "6-bus case" begin - result = run_pf("../test/data/matpower/case6.m", SOCWRPowerModel, ipopt_solver, solution_processors=[sol_data_model!]) + result = run_pf("../test/data/matpower/case6.m", SOCWRPowerModel, nlp_solver, solution_processors=[sol_data_model!]) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 0; atol = 1e-2) @@ -278,7 +278,7 @@ end @test isapprox(result["solution"]["bus"]["4"]["vm"], 1.09999; atol = 1e-3) end @testset "24-bus rts case" begin - result = run_pf("../test/data/matpower/case24.m", SOCWRPowerModel, ipopt_solver) + result = run_pf("../test/data/matpower/case24.m", SOCWRPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 0; atol = 1e-2) @@ -289,7 +289,7 @@ end @testset "test soc distflow pf_bf" begin @testset "3-bus case" begin - result = run_pf_bf("../test/data/matpower/case3.m", SOCBFPowerModel, ipopt_solver, solution_processors=[sol_data_model!]) + result = run_pf_bf("../test/data/matpower/case3.m", SOCBFPowerModel, nlp_solver, solution_processors=[sol_data_model!]) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 0; atol = 1e-2) @@ -307,19 +307,19 @@ end @test isapprox(result["solution"]["dcline"]["1"]["pt"], -0.10; atol = 1e-4) end @testset "5-bus asymmetric case" begin - result = run_pf_bf("../test/data/matpower/case5_asym.m", SOCBFPowerModel, ipopt_solver) + result = run_pf_bf("../test/data/matpower/case5_asym.m", SOCBFPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 0; atol = 1e-2) end @testset "5-bus case with hvdc line" begin - result = run_pf_bf("../test/data/matpower/case5_dc.m", SOCBFPowerModel, ipopt_solver) + result = run_pf_bf("../test/data/matpower/case5_dc.m", SOCBFPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 0; atol = 1e-2) end @testset "6-bus case" begin - result = run_pf_bf("../test/data/matpower/case6.m", SOCBFPowerModel, ipopt_solver, solution_processors=[sol_data_model!]) + result = run_pf_bf("../test/data/matpower/case6.m", SOCBFPowerModel, nlp_solver, solution_processors=[sol_data_model!]) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 0; atol = 1e-2) @@ -327,7 +327,7 @@ end @test isapprox(result["solution"]["bus"]["4"]["vm"], 1.09999; atol = 1e-3) end @testset "24-bus rts case" begin - result = run_pf_bf("../test/data/matpower/case24.m", SOCBFPowerModel, ipopt_solver) + result = run_pf_bf("../test/data/matpower/case24.m", SOCBFPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 0; atol = 1e-2) @@ -337,7 +337,7 @@ end @testset "test linear distflow pf_bf" begin @testset "3-bus case" begin - result = run_pf_bf("../test/data/matpower/case3.m", BFAPowerModel, ipopt_solver, solution_processors=[sol_data_model!]) + result = run_pf_bf("../test/data/matpower/case3.m", BFAPowerModel, nlp_solver, solution_processors=[sol_data_model!]) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 0; atol = 1e-2) @@ -355,19 +355,19 @@ end @test isapprox(result["solution"]["dcline"]["1"]["pt"], -0.10; atol = 1e-4) end @testset "5-bus asymmetric case" begin - result = run_pf_bf("../test/data/matpower/case5_asym.m", BFAPowerModel, ipopt_solver) + result = run_pf_bf("../test/data/matpower/case5_asym.m", BFAPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 0; atol = 1e-2) end @testset "5-bus case with hvdc line" begin - result = run_pf_bf("../test/data/matpower/case5_dc.m", BFAPowerModel, ipopt_solver) + result = run_pf_bf("../test/data/matpower/case5_dc.m", BFAPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 0; atol = 1e-2) end @testset "6-bus case" begin - result = run_pf_bf("../test/data/matpower/case6.m", BFAPowerModel, ipopt_solver, solution_processors=[sol_data_model!]) + result = run_pf_bf("../test/data/matpower/case6.m", BFAPowerModel, nlp_solver, solution_processors=[sol_data_model!]) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 0; atol = 1e-2) @@ -376,7 +376,7 @@ end end @testset "24-bus rts case" begin data = parse_file("../test/data/matpower/case24.m") - result = run_pf_bf(data, BFAPowerModel, ipopt_solver) + result = run_pf_bf(data, BFAPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 0; atol = 1e-2) @@ -390,7 +390,7 @@ end @testset "test sdp pf" begin # note: may have issues on linux (04/02/18) @testset "3-bus case" begin - result = run_pf("../test/data/matpower/case3.m", SDPWRMPowerModel, scs_solver, solution_processors=[sol_data_model!]) + result = run_pf("../test/data/matpower/case3.m", SDPWRMPowerModel, sdp_solver, solution_processors=[sol_data_model!]) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 0; atol = 1e-2) @@ -409,13 +409,13 @@ end end # note: may have issues on os x (05/07/18) @testset "5-bus asymmetric case" begin - result = run_pf("../test/data/matpower/case5_asym.m", SDPWRMPowerModel, scs_solver) + result = run_pf("../test/data/matpower/case5_asym.m", SDPWRMPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 0; atol = 1e-2) end @testset "6-bus case" begin - result = run_pf("../test/data/matpower/case6.m", SDPWRMPowerModel, scs_solver, solution_processors=[sol_data_model!]) + result = run_pf("../test/data/matpower/case6.m", SDPWRMPowerModel, sdp_solver, solution_processors=[sol_data_model!]) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 0; atol = 1e-2) @@ -423,7 +423,7 @@ end @test isapprox(result["solution"]["bus"]["4"]["vm"], 1.09999; atol = 1e-3) end @testset "24-bus rts case" begin - result = run_pf("../test/data/matpower/case24.m", SDPWRMPowerModel, scs_solver) + result = run_pf("../test/data/matpower/case24.m", SDPWRMPowerModel, sdp_solver) @test result["termination_status"] == OPTIMAL @test isapprox(result["objective"], 0; atol = 1e-2) diff --git a/test/psse.jl b/test/psse.jl index 393a83d19..e355f4a04 100644 --- a/test/psse.jl +++ b/test/psse.jl @@ -28,8 +28,8 @@ end set_costs!(data_mp) - result_pti = PowerModels.run_opf(data_pti, PowerModels.ACPPowerModel, ipopt_solver) - result_mp = PowerModels.run_opf(data_mp, PowerModels.ACPPowerModel, ipopt_solver) + result_pti = PowerModels.run_opf(data_pti, PowerModels.ACPPowerModel, nlp_solver) + result_mp = PowerModels.run_opf(data_mp, PowerModels.ACPPowerModel, nlp_solver) @test result_pti["termination_status"] == LOCALLY_SOLVED @test result_mp["termination_status"] == LOCALLY_SOLVED @@ -42,8 +42,8 @@ end set_costs!(data_mp) - result_pti = PowerModels.run_opf(data_pti, PowerModels.ACPPowerModel, ipopt_solver) - result_mp = PowerModels.run_opf(data_mp, PowerModels.ACPPowerModel, ipopt_solver) + result_pti = PowerModels.run_opf(data_pti, PowerModels.ACPPowerModel, nlp_solver) + result_mp = PowerModels.run_opf(data_mp, PowerModels.ACPPowerModel, nlp_solver) @test result_pti["termination_status"] == LOCALLY_SOLVED @test result_mp["termination_status"] == LOCALLY_SOLVED @@ -58,8 +58,8 @@ end set_costs!(data_mp) - result_pti = PowerModels.run_opf(data_pti, PowerModels.ACPPowerModel, ipopt_solver) - result_mp = PowerModels.run_opf(data_mp, PowerModels.ACPPowerModel, ipopt_solver) + result_pti = PowerModels.run_opf(data_pti, PowerModels.ACPPowerModel, nlp_solver) + result_mp = PowerModels.run_opf(data_mp, PowerModels.ACPPowerModel, nlp_solver) @test result_pti["termination_status"] == LOCALLY_SOLVED @test result_mp["termination_status"] == LOCALLY_SOLVED @@ -82,12 +82,12 @@ end end end - result_opf = PowerModels.run_opf(data_pti, PowerModels.ACPPowerModel, ipopt_solver) + result_opf = PowerModels.run_opf(data_pti, PowerModels.ACPPowerModel, nlp_solver) @test result_opf["termination_status"] == LOCALLY_SOLVED @test isapprox(result_opf["objective"], 29.4043; atol=1e-4) - result_pf = PowerModels.run_pf(data_pti, PowerModels.ACPPowerModel, ipopt_solver) + result_pf = PowerModels.run_pf(data_pti, PowerModels.ACPPowerModel, nlp_solver) for (bus, vm, va) in zip(["1002", "1005", "1008", "1009"], [1.0032721, 1.0199983, 1.0203627, 1.03], @@ -105,8 +105,8 @@ end set_costs!(data_mp) - result_pti = PowerModels.run_opf(data_pti, PowerModels.ACPPowerModel, ipopt_solver) - result_mp = PowerModels.run_opf(data_mp, PowerModels.ACPPowerModel, ipopt_solver) + result_pti = PowerModels.run_opf(data_pti, PowerModels.ACPPowerModel, nlp_solver) + result_mp = PowerModels.run_opf(data_mp, PowerModels.ACPPowerModel, nlp_solver) @test result_pti["termination_status"] == LOCALLY_SOLVED @test result_mp["termination_status"] == LOCALLY_SOLVED @@ -123,8 +123,8 @@ end set_costs!(data_mp) - result_pti = PowerModels.run_opf(data_pti, PowerModels.ACPPowerModel, ipopt_solver) - result_mp = PowerModels.run_opf(data_mp, PowerModels.ACPPowerModel, ipopt_solver) + result_pti = PowerModels.run_opf(data_pti, PowerModels.ACPPowerModel, nlp_solver) + result_mp = PowerModels.run_opf(data_mp, PowerModels.ACPPowerModel, nlp_solver) @test result_pti["termination_status"] == LOCALLY_SOLVED @test result_mp["termination_status"] == LOCALLY_SOLVED @@ -143,8 +143,8 @@ end set_costs!(data_mp) - result_pti = PowerModels.run_opf(data_pti, PowerModels.ACPPowerModel, ipopt_solver) - result_mp = PowerModels.run_opf(data_mp, PowerModels.ACPPowerModel, ipopt_solver) + result_pti = PowerModels.run_opf(data_pti, PowerModels.ACPPowerModel, nlp_solver) + result_mp = PowerModels.run_opf(data_mp, PowerModels.ACPPowerModel, nlp_solver) @test result_pti["termination_status"] == LOCALLY_SOLVED @test result_mp["termination_status"] == LOCALLY_SOLVED @@ -161,8 +161,8 @@ end set_costs!(data_mp) - result_pti = PowerModels.run_opf(data_pti, PowerModels.ACPPowerModel, ipopt_solver) - result_mp = PowerModels.run_opf(data_mp, PowerModels.ACPPowerModel, ipopt_solver) + result_pti = PowerModels.run_opf(data_pti, PowerModels.ACPPowerModel, nlp_solver) + result_mp = PowerModels.run_opf(data_mp, PowerModels.ACPPowerModel, nlp_solver) @test result_pti["termination_status"] == LOCALLY_SOLVED @test result_mp["termination_status"] == LOCALLY_SOLVED @@ -178,8 +178,8 @@ end set_costs!(data_mp) - result_pti = PowerModels.run_opf(data_pti, PowerModels.ACPPowerModel, ipopt_solver) - result_mp = PowerModels.run_opf(data_mp, PowerModels.ACPPowerModel, ipopt_solver) + result_pti = PowerModels.run_opf(data_pti, PowerModels.ACPPowerModel, nlp_solver) + result_mp = PowerModels.run_opf(data_mp, PowerModels.ACPPowerModel, nlp_solver) @test result_pti["termination_status"] == LOCALLY_SOLVED @test result_mp["termination_status"] == LOCALLY_SOLVED @@ -196,8 +196,8 @@ end set_costs!(data_mp) - result_pti = PowerModels.run_opf(data_pti, PowerModels.ACPPowerModel, ipopt_solver) - result_mp = PowerModels.run_opf(data_mp, PowerModels.ACPPowerModel, ipopt_solver) + result_pti = PowerModels.run_opf(data_pti, PowerModels.ACPPowerModel, nlp_solver) + result_mp = PowerModels.run_opf(data_mp, PowerModels.ACPPowerModel, nlp_solver) @test result_pti["termination_status"] == LOCALLY_SOLVED @test result_mp["termination_status"] == LOCALLY_SOLVED @@ -259,12 +259,12 @@ end @test length(data_pti["bus"]) == 8 @test length(data_pti["branch"]) == 15 - result_opf = PowerModels.run_opf(data_pti, PowerModels.ACPPowerModel, ipopt_solver) + result_opf = PowerModels.run_opf(data_pti, PowerModels.ACPPowerModel, nlp_solver) @test result_opf["termination_status"] == LOCALLY_SOLVED @test isapprox(result_opf["objective"], 10.00027; atol=1e-5) - result_pf = PowerModels.run_pf(data_pti, PowerModels.ACPPowerModel, ipopt_solver) + result_pf = PowerModels.run_pf(data_pti, PowerModels.ACPPowerModel, nlp_solver) bus_data = zip( ["1001", "1002", "1003", "10002", "10003", "10004"], @@ -303,12 +303,12 @@ end end - result_opf = PowerModels.run_opf(data_pti, PowerModels.ACPPowerModel, ipopt_solver) + result_opf = PowerModels.run_opf(data_pti, PowerModels.ACPPowerModel, nlp_solver) @test result_opf["termination_status"] == LOCALLY_SOLVED @test isapprox(result_opf["objective"], 10.0; atol=1e-5) - result_pf = PowerModels.run_pf(data_pti, PowerModels.ACPPowerModel, ipopt_solver) + result_pf = PowerModels.run_pf(data_pti, PowerModels.ACPPowerModel, nlp_solver) for (bus, vm, va) in zip(["1001", "1002", "1003", "10001"], [1.09, 1.0, 1.0, 0.997], [2.304, 0., 6.042244, 2.5901]) @test isapprox(result_pf["solution"]["bus"][bus]["vm"], vm; atol=1e-1) @@ -327,12 +327,12 @@ end @test isapprox(data_pti["branch"]["1"]["g_fr"], 5e-3; atol=1e-4) @test isapprox(data_pti["branch"]["1"]["b_fr"], 6.74e-3; atol=1e-4) - result_opf = PowerModels.run_opf(data_pti, PowerModels.ACPPowerModel, ipopt_solver) + result_opf = PowerModels.run_opf(data_pti, PowerModels.ACPPowerModel, nlp_solver) @test result_opf["termination_status"] == LOCALLY_SOLVED @test isapprox(result_opf["objective"], 701.637157; atol=1e-5) - result_pf = PowerModels.run_pf(data_pti, PowerModels.ACPPowerModel, ipopt_solver) + result_pf = PowerModels.run_pf(data_pti, PowerModels.ACPPowerModel, nlp_solver) @test result_pf["termination_status"] == LOCALLY_SOLVED @test result_pf["objective"] == 0.0 @@ -351,12 +351,12 @@ end @test isapprox(data_pti["branch"]["1"]["g_fr"], 5e-3; atol=1e-4) @test isapprox(data_pti["branch"]["1"]["b_fr"], 6.74e-3; atol=1e-4) - result_opf = PowerModels.run_opf(data_pti, PowerModels.ACPPowerModel, ipopt_solver) + result_opf = PowerModels.run_opf(data_pti, PowerModels.ACPPowerModel, nlp_solver) @test result_opf["termination_status"] == LOCALLY_SOLVED @test isapprox(result_opf["objective"], 10.4001; atol=1e-2) - result_pf = PowerModels.run_pf(data_pti, PowerModels.ACPPowerModel, ipopt_solver) + result_pf = PowerModels.run_pf(data_pti, PowerModels.ACPPowerModel, nlp_solver) @test result_pf["termination_status"] == LOCALLY_SOLVED @test result_pf["objective"] == 0.0 @@ -412,7 +412,7 @@ end end end - result = PowerModels.run_opf(data, PowerModels.ACPPowerModel, ipopt_solver) + result = PowerModels.run_opf(data, PowerModels.ACPPowerModel, nlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["objective"], 297.878089; atol=1e-4) @@ -479,11 +479,11 @@ end @test length(data["dcline"]) == 1 @test length(data["dcline"]["1"]) == 26 - opf = PowerModels.run_opf(data, PowerModels.ACPPowerModel, ipopt_solver) + opf = PowerModels.run_opf(data, PowerModels.ACPPowerModel, nlp_solver) @test opf["termination_status"] == LOCALLY_SOLVED @test isapprox(opf["objective"], 10.5; atol=1e-3) - pf = PowerModels.run_pf(data, PowerModels.ACPPowerModel, ipopt_solver) + pf = PowerModels.run_pf(data, PowerModels.ACPPowerModel, nlp_solver) @test pf["termination_status"] == LOCALLY_SOLVED end @@ -493,11 +493,11 @@ end @test length(data["dcline"]) == 1 @test length(data["dcline"]["1"]) == 26 - opf = PowerModels.run_opf(data, PowerModels.ACPPowerModel, ipopt_solver) + opf = PowerModels.run_opf(data, PowerModels.ACPPowerModel, nlp_solver) @test opf["termination_status"] == LOCALLY_SOLVED @test isapprox(opf["objective"], 21.8842; atol=1e-3) - pf = PowerModels.run_pf(data, PowerModels.ACPPowerModel, ipopt_solver) + pf = PowerModels.run_pf(data, PowerModels.ACPPowerModel, nlp_solver) @test pf["termination_status"] == LOCALLY_SOLVED end end diff --git a/test/pti.jl b/test/pti.jl index 071c14fec..242ef285e 100644 --- a/test/pti.jl +++ b/test/pti.jl @@ -344,12 +344,12 @@ end function compare_pf(filename::String) source_data = parse_file(filename) - source_solution = PowerModels.run_pf(source_data, ACPPowerModel, ipopt_solver)["solution"] + source_solution = PowerModels.run_pf(source_data, ACPPowerModel, nlp_solver)["solution"] file_tmp = "../test/data/tmp.raw" PowerModels.export_pti(file_tmp, source_data) destination_data = parse_file(file_tmp) - destination_solution = PowerModels.run_pf(destination_data, ACPPowerModel, ipopt_solver)["solution"] + destination_solution = PowerModels.run_pf(destination_data, ACPPowerModel, nlp_solver)["solution"] rm(file_tmp) @test InfrastructureModels.compare_dict(source_solution, destination_solution) diff --git a/test/runtests.jl b/test/runtests.jl index f18113fca..49c7902f0 100644 --- a/test/runtests.jl +++ b/test/runtests.jl @@ -29,12 +29,12 @@ end # default setup for solvers -ipopt_solver = JuMP.optimizer_with_attributes(Ipopt.Optimizer, "tol"=>1e-6, "print_level"=>0) -ipopt_ws_solver = JuMP.optimizer_with_attributes(Ipopt.Optimizer, "tol"=>1e-6, "mu_init"=>1e-4, "print_level"=>0) +nlp_solver = JuMP.optimizer_with_attributes(Ipopt.Optimizer, "tol"=>1e-6, "print_level"=>0) +nlp_ws_solver = JuMP.optimizer_with_attributes(Ipopt.Optimizer, "tol"=>1e-6, "mu_init"=>1e-4, "print_level"=>0) -cbc_solver = JuMP.optimizer_with_attributes(HiGHS.Optimizer, "output_flag"=>false) -juniper_solver = JuMP.optimizer_with_attributes(Juniper.Optimizer, "nl_solver"=>JuMP.optimizer_with_attributes(Ipopt.Optimizer, "tol"=>1e-4, "print_level"=>0), "log_levels"=>[]) -scs_solver = JuMP.optimizer_with_attributes(SCS.Optimizer, "verbose"=>false) +milp_solver = JuMP.optimizer_with_attributes(HiGHS.Optimizer, "output_flag"=>false) +minlp_solver = JuMP.optimizer_with_attributes(Juniper.Optimizer, "nl_solver"=>JuMP.optimizer_with_attributes(Ipopt.Optimizer, "tol"=>1e-4, "print_level"=>0), "log_levels"=>[]) +sdp_solver = JuMP.optimizer_with_attributes(SCS.Optimizer, "verbose"=>false) include("common.jl") diff --git a/test/tnep.jl b/test/tnep.jl index a00c01cec..34642ef88 100644 --- a/test/tnep.jl +++ b/test/tnep.jl @@ -11,7 +11,7 @@ end @testset "3-bus case" begin data = PowerModels.parse_file("../test/data/matpower/case3_tnep.m") calc_thermal_limits!(data) - result = run_tnep(data, ACPPowerModel, juniper_solver) + result = run_tnep(data, ACPPowerModel, minlp_solver) check_tnep_status(result["solution"]) @@ -20,7 +20,7 @@ end end @testset "5-bus case" begin - result = run_tnep("../test/data/matpower/case5_tnep.m", ACPPowerModel, juniper_solver) + result = run_tnep("../test/data/matpower/case5_tnep.m", ACPPowerModel, minlp_solver) check_tnep_status(result["solution"]) @@ -34,7 +34,7 @@ end @testset "3-bus case" begin data = PowerModels.parse_file("../test/data/matpower/case3_tnep.m") calc_thermal_limits!(data) - result = run_tnep(data, SOCWRPowerModel, juniper_solver) + result = run_tnep(data, SOCWRPowerModel, minlp_solver) check_tnep_status(result["solution"]) @@ -43,7 +43,7 @@ end end @testset "5-bus rts case" begin - result = run_tnep("../test/data/matpower/case5_tnep.m", SOCWRPowerModel, juniper_solver) + result = run_tnep("../test/data/matpower/case5_tnep.m", SOCWRPowerModel, minlp_solver) check_tnep_status(result["solution"]) @@ -57,7 +57,7 @@ end @testset "3-bus case" begin data = PowerModels.parse_file("../test/data/matpower/case3_tnep.m") calc_thermal_limits!(data) - result = run_tnep(data, QCRMPowerModel, juniper_solver) + result = run_tnep(data, QCRMPowerModel, minlp_solver) check_tnep_status(result["solution"]) @@ -66,7 +66,7 @@ end end @testset "5-bus rts case" begin - result = run_tnep("../test/data/matpower/case5_tnep.m", QCRMPowerModel, juniper_solver) + result = run_tnep("../test/data/matpower/case5_tnep.m", QCRMPowerModel, minlp_solver) check_tnep_status(result["solution"]) @@ -80,7 +80,7 @@ end @testset "3-bus case" begin data = PowerModels.parse_file("../test/data/matpower/case3_tnep.m") calc_thermal_limits!(data) - result = run_tnep(data, DCPPowerModel, juniper_solver) + result = run_tnep(data, DCPPowerModel, minlp_solver) check_tnep_status(result["solution"]) @@ -89,7 +89,7 @@ end end @testset "5-bus case" begin - result = run_tnep("../test/data/matpower/case5_tnep.m", DCPPowerModel, juniper_solver) + result = run_tnep("../test/data/matpower/case5_tnep.m", DCPPowerModel, minlp_solver) check_tnep_status(result["solution"]) @@ -101,7 +101,7 @@ end @testset "test matpower dc tnep" begin @testset "5-bus case with matpower DCMP model and TNEP" begin - result = run_tnep("../test/data/matpower/case5_tnep.m", DCMPPowerModel, juniper_solver) + result = run_tnep("../test/data/matpower/case5_tnep.m", DCMPPowerModel, minlp_solver) @test result["termination_status"] == LOCALLY_SOLVED @test isapprox(result["solution"]["ne_branch"]["1"]["built"], 1.0; atol = 1e-5) @@ -116,7 +116,7 @@ end @testset "3-bus case" begin data = PowerModels.parse_file("../test/data/matpower/case3_tnep.m") calc_thermal_limits!(data) - result = run_tnep(data, DCPLLPowerModel, juniper_solver) + result = run_tnep(data, DCPLLPowerModel, minlp_solver) check_tnep_status(result["solution"]) @@ -126,7 +126,7 @@ end =# @testset "5-bus case" begin - result = run_tnep("../test/data/matpower/case5_tnep.m", DCPLLPowerModel, juniper_solver) + result = run_tnep("../test/data/matpower/case5_tnep.m", DCPLLPowerModel, minlp_solver) check_tnep_status(result["solution"]) @@ -139,7 +139,7 @@ end @testset "3-bus case" begin data = PowerModels.parse_file("../test/data/matpower/case3_tnep.m") calc_thermal_limits!(data) - result = run_tnep(data, LPACCPowerModel, juniper_solver) + result = run_tnep(data, LPACCPowerModel, minlp_solver) check_tnep_status(result["solution"]) @@ -148,7 +148,7 @@ end end @testset "5-bus case" begin - result = run_tnep("../test/data/matpower/case5_tnep.m", LPACCPowerModel, juniper_solver) + result = run_tnep("../test/data/matpower/case5_tnep.m", LPACCPowerModel, minlp_solver) check_tnep_status(result["solution"]) @@ -161,7 +161,7 @@ end @testset "3-bus case" begin data = PowerModels.parse_file("../test/data/matpower/case3_tnep.m") calc_thermal_limits!(data) - result = run_tnep(data, SOCWRPowerModel, juniper_solver) + result = run_tnep(data, SOCWRPowerModel, minlp_solver) check_tnep_status(result["solution"]) diff --git a/test/util.jl b/test/util.jl index 85a858382..a9d916301 100644 --- a/test/util.jl +++ b/test/util.jl @@ -1,15 +1,15 @@ @testset "obbt with trilinear convex hull relaxation" begin @testset "3-bus case" begin - result_ac = run_ac_opf("../test/data/matpower/case3.m", ipopt_solver); + result_ac = run_ac_opf("../test/data/matpower/case3.m", nlp_solver); upper_bound = result_ac["objective"] - data, stats = run_obbt_opf!("../test/data/matpower/case3.m", ipopt_solver, model_type=QCLSPowerModel); + data, stats = run_obbt_opf!("../test/data/matpower/case3.m", nlp_solver, model_type=QCLSPowerModel); @test isapprox(stats["final_relaxation_objective"], 5901.96; atol=1e0) @test isnan(stats["final_rel_gap_from_ub"]) @test stats["iteration_count"] == 5 - data, stats = run_obbt_opf!("../test/data/matpower/case3.m", ipopt_solver, + data, stats = run_obbt_opf!("../test/data/matpower/case3.m", nlp_solver, model_type = QCLSPowerModel, upper_bound = upper_bound, upper_bound_constraint = true, @@ -23,15 +23,15 @@ end @testset "obbt with qc relaxation" begin @testset "3-bus case" begin - result_ac = run_ac_opf("../test/data/matpower/case3.m", ipopt_solver); + result_ac = run_ac_opf("../test/data/matpower/case3.m", nlp_solver); upper_bound = result_ac["objective"] - data, stats = run_obbt_opf!("../test/data/matpower/case3.m", ipopt_solver, model_type=QCRMPowerModel); + data, stats = run_obbt_opf!("../test/data/matpower/case3.m", nlp_solver, model_type=QCRMPowerModel); @test isapprox(stats["final_relaxation_objective"], 5900.04; atol=1e0) @test isnan(stats["final_rel_gap_from_ub"]) @test stats["iteration_count"] == 5 - data, stats = run_obbt_opf!("../test/data/matpower/case3.m", ipopt_solver, + data, stats = run_obbt_opf!("../test/data/matpower/case3.m", nlp_solver, model_type = QCRMPowerModel, upper_bound = upper_bound, upper_bound_constraint = true, @@ -47,10 +47,10 @@ end data["gen"]["1"]["cost"] = data["gen"]["1"]["cost"][2:3] data["gen"]["2"]["cost"] = data["gen"]["2"]["cost"][2:3] - result_ac = run_ac_opf(data, ipopt_solver); + result_ac = run_ac_opf(data, nlp_solver); upper_bound = result_ac["objective"] - data, stats = run_obbt_opf!(data, ipopt_solver, + data, stats = run_obbt_opf!(data, nlp_solver, model_type=QCRMPowerModel, upper_bound = upper_bound, upper_bound_constraint = true); @@ -63,8 +63,8 @@ end @testset "opf with flow cuts" begin @testset "ac 5-bus case" begin - result_base = run_opf("../test/data/matpower/case5.m", ACPPowerModel, ipopt_solver) - result_cuts = run_opf_branch_power_cuts("../test/data/matpower/case5.m", ACPPowerModel, ipopt_solver) + result_base = run_opf("../test/data/matpower/case5.m", ACPPowerModel, nlp_solver) + result_cuts = run_opf_branch_power_cuts("../test/data/matpower/case5.m", ACPPowerModel, nlp_solver) @test result_base["termination_status"] == LOCALLY_SOLVED @test result_cuts["termination_status"] == LOCALLY_SOLVED @@ -75,8 +75,8 @@ end end end @testset "ac 14-bus case" begin - result_base = run_opf("../test/data/matpower/case14.m", ACPPowerModel, ipopt_solver) - result_cuts = run_opf_branch_power_cuts("../test/data/matpower/case14.m", ACPPowerModel, ipopt_solver) + result_base = run_opf("../test/data/matpower/case14.m", ACPPowerModel, nlp_solver) + result_cuts = run_opf_branch_power_cuts("../test/data/matpower/case14.m", ACPPowerModel, nlp_solver) @test result_base["termination_status"] == LOCALLY_SOLVED @test result_cuts["termination_status"] == LOCALLY_SOLVED @@ -88,8 +88,8 @@ end end @testset "soc 5-bus case" begin - result_base = run_opf("../test/data/matpower/case5.m", SOCWRPowerModel, ipopt_solver) - result_cuts = run_opf_branch_power_cuts("../test/data/matpower/case5.m", SOCWRPowerModel, ipopt_solver) + result_base = run_opf("../test/data/matpower/case5.m", SOCWRPowerModel, nlp_solver) + result_cuts = run_opf_branch_power_cuts("../test/data/matpower/case5.m", SOCWRPowerModel, nlp_solver) @test result_base["termination_status"] == LOCALLY_SOLVED @test result_cuts["termination_status"] == LOCALLY_SOLVED @@ -99,8 +99,8 @@ end end end @testset "soc 14-bus case" begin - result_base = run_opf("../test/data/matpower/case14.m", SOCWRPowerModel, ipopt_solver) - result_cuts = run_opf_branch_power_cuts("../test/data/matpower/case14.m", SOCWRPowerModel, ipopt_solver) + result_base = run_opf("../test/data/matpower/case14.m", SOCWRPowerModel, nlp_solver) + result_cuts = run_opf_branch_power_cuts("../test/data/matpower/case14.m", SOCWRPowerModel, nlp_solver) @test result_base["termination_status"] == LOCALLY_SOLVED @test result_cuts["termination_status"] == LOCALLY_SOLVED @@ -111,8 +111,8 @@ end end @testset "dc 5-bus case" begin - result_base = run_opf("../test/data/matpower/case5.m", DCPPowerModel, cbc_solver) - result_cuts = run_opf_branch_power_cuts("../test/data/matpower/case5.m", DCPPowerModel, cbc_solver) + result_base = run_opf("../test/data/matpower/case5.m", DCPPowerModel, milp_solver) + result_cuts = run_opf_branch_power_cuts("../test/data/matpower/case5.m", DCPPowerModel, milp_solver) @test result_base["termination_status"] == OPTIMAL @test result_cuts["termination_status"] == OPTIMAL @@ -122,8 +122,8 @@ end end end @testset "dc 14-bus case" begin - result_base = run_opf("../test/data/matpower/case14.m", DCPPowerModel, ipopt_solver) - result_cuts = run_opf_branch_power_cuts("../test/data/matpower/case14.m", DCPPowerModel, ipopt_solver) + result_base = run_opf("../test/data/matpower/case14.m", DCPPowerModel, nlp_solver) + result_cuts = run_opf_branch_power_cuts("../test/data/matpower/case14.m", DCPPowerModel, nlp_solver) @test result_base["termination_status"] == LOCALLY_SOLVED @test result_cuts["termination_status"] == LOCALLY_SOLVED @@ -137,8 +137,8 @@ end @testset "ptdf opf with flow cuts" begin @testset "dc 5-bus case" begin - result_base = run_opf("../test/data/matpower/case5.m", DCPPowerModel, cbc_solver) - result_cuts = run_opf_ptdf_branch_power_cuts("../test/data/matpower/case5.m", cbc_solver) + result_base = run_opf("../test/data/matpower/case5.m", DCPPowerModel, milp_solver) + result_cuts = run_opf_ptdf_branch_power_cuts("../test/data/matpower/case5.m", milp_solver) @test result_base["termination_status"] == OPTIMAL @test result_cuts["termination_status"] == OPTIMAL @@ -148,8 +148,8 @@ end end end @testset "dc 5-bus ext case" begin - result_base = run_opf("../test/data/matpower/case5_ext.m", DCPPowerModel, cbc_solver) - result_cuts = run_opf_ptdf_branch_power_cuts("../test/data/matpower/case5_ext.m", cbc_solver) + result_base = run_opf("../test/data/matpower/case5_ext.m", DCPPowerModel, milp_solver) + result_cuts = run_opf_ptdf_branch_power_cuts("../test/data/matpower/case5_ext.m", milp_solver) @test result_base["termination_status"] == OPTIMAL @test result_cuts["termination_status"] == OPTIMAL @@ -159,8 +159,8 @@ end end end @testset "dc 14-bus case" begin - result_base = run_opf("../test/data/matpower/case14.m", DCPPowerModel, ipopt_solver) - result_cuts = run_opf_ptdf_branch_power_cuts("../test/data/matpower/case14.m", ipopt_solver) + result_base = run_opf("../test/data/matpower/case14.m", DCPPowerModel, nlp_solver) + result_cuts = run_opf_ptdf_branch_power_cuts("../test/data/matpower/case14.m", nlp_solver) @test result_base["termination_status"] == LOCALLY_SOLVED @test result_cuts["termination_status"] == LOCALLY_SOLVED diff --git a/test/warmstart.jl b/test/warmstart.jl index 18a00d458..e956b357e 100644 --- a/test/warmstart.jl +++ b/test/warmstart.jl @@ -40,32 +40,32 @@ end @testset "dc warm starts" begin @testset "5 bus case" begin data = PowerModels.parse_file("../test/data/matpower/case5.m") - result = run_dc_opf(data, ipopt_solver; setting = Dict("output" => Dict("branch_flows" => true))); + result = run_dc_opf(data, nlp_solver; setting = Dict("output" => Dict("branch_flows" => true))); PowerModels.update_data!(data, result["solution"]) # 14 iterations - result = run_dc_opf(data, ipopt_solver); + result = run_dc_opf(data, nlp_solver); set_dc_start!(data) # 6 iterations - result = run_dc_opf(data, ipopt_ws_solver); + result = run_dc_opf(data, nlp_ws_solver); end @testset "5 bus pwl case" begin data = PowerModels.parse_file("../test/data/matpower/case5_pwlc.m") - result = run_dc_opf(data, ipopt_solver; setting = Dict("output" => Dict("branch_flows" => true))); + result = run_dc_opf(data, nlp_solver; setting = Dict("output" => Dict("branch_flows" => true))); PowerModels.update_data!(data, result["solution"]) # 35 iterations - result = run_dc_opf(data, ipopt_solver); + result = run_dc_opf(data, nlp_solver); set_dc_start!(data) # 6 iterations - result = run_dc_opf(data, ipopt_ws_solver); + result = run_dc_opf(data, nlp_ws_solver); end end @@ -73,32 +73,32 @@ end @testset "ac warm starts" begin @testset "5 bus case" begin data = PowerModels.parse_file("../test/data/matpower/case5.m") - result = run_ac_opf(data, ipopt_solver; setting = Dict("output" => Dict("branch_flows" => true))); + result = run_ac_opf(data, nlp_solver; setting = Dict("output" => Dict("branch_flows" => true))); PowerModels.update_data!(data, result["solution"]) # 22 iterations - result = run_ac_opf(data, ipopt_solver); + result = run_ac_opf(data, nlp_solver); set_ac_start!(data) # 19 iterations - result = run_ac_opf(data, ipopt_ws_solver); + result = run_ac_opf(data, nlp_ws_solver); end @testset "5 bus pwl case" begin data = PowerModels.parse_file("../test/data/matpower/case5_pwlc.m") - result = run_ac_opf(data, ipopt_solver; setting = Dict("output" => Dict("branch_flows" => true))); + result = run_ac_opf(data, nlp_solver; setting = Dict("output" => Dict("branch_flows" => true))); PowerModels.update_data!(data, result["solution"]) # 40 iterations - result = run_ac_opf(data, ipopt_solver); + result = run_ac_opf(data, nlp_solver); set_ac_start!(data) # 12 iterations - result = run_ac_opf(data, ipopt_ws_solver); + result = run_ac_opf(data, nlp_ws_solver); end end