Fix element type of Dict in pf.jl

CHANGELOG.md

@@ -2,7 +2,7 @@ PowerModels.jl Change Log
 =========================
 
 ### Staged
-- nothing
+- Fix InexactError in `compute_ac_pf` (#939)
 
 ### v0.21.3
 - Fix no-buses bug in `calc_connected_components` (#933)

src/prob/pf.jl

@@ -312,7 +312,7 @@ function compute_ac_pf(pf_data::PowerFlowData; kwargs...)
             if bus["bus_type"] != 4
                 bus_idx = am.bus_to_idx[bus["index"]]
 
-                bus_assignment[i] = Dict(
+                bus_assignment[i] = Dict{String,Float64}(
                     "vm" => pf_data.vm_idx[bus_idx],
                     "va" => pf_data.va_idx[bus_idx]
                 )
@@ -322,7 +322,7 @@ function compute_ac_pf(pf_data::PowerFlowData; kwargs...)
         gen_assignment= Dict{String,Any}()
         for (i,gen) in data["gen"]
             if gen["gen_status"] != 0
-                gen_assignment[i] = Dict(
+                gen_assignment[i] = Dict{String,Float64}(
                     "pg" => gen["pg"],
                     "qg" => gen["qg"]
                 )

test/data/json/issue_938.json

@@ -0,0 +1,118 @@
+{
+    "baseMVA": 1000,
+    "per_unit": true,
+    "shunt": {},
+    "storage": {},
+    "switch": {},
+    "dcline": {},
+    "load": {
+        "1": {
+            "load_bus": 1,
+            "pd": 0.1,
+            "qd": 0,
+            "index": 1,
+            "status": 1
+        }
+    },
+    "gen": {
+        "2": {
+            "gen_bus": 2,
+            "pg": 0.1,
+            "qg": 0,
+            "pmax": 1,
+            "pmin": 0,
+            "qmax": 9999,
+            "qmin": -9999,
+            "index": 2,
+            "vg": 1,
+            "gen_status": 1
+        },
+        "3": {
+            "gen_bus": 3,
+            "pg": 0,
+            "qg": 0,
+            "pmax": 9999,
+            "pmin": 0,
+            "qmax": 9999,
+            "qmin": -9999,
+            "index": 3,
+            "vg": 1,
+            "gen_status": 1
+        }
+    },
+    "branch": {
+        "5": {
+            "index": 5,
+            "f_bus": 1,
+            "t_bus": 2,
+            "br_r": 0.06815,
+            "br_x": 1.096,
+            "b_fr": 0.03092,
+            "b_to": 0.03092,
+            "br_status": 1,
+            "base_kv": 220,
+            "c_rating_a": 1.959,
+            "angmin": -0.17453292519943295,
+            "angmax": 0.17453292519943295,
+            "transformer": false,
+            "tap": 1,
+            "shift": 0,
+            "g_fr": 0,
+            "g_to": 0
+        },
+        "6": {
+            "index": 6,
+            "f_bus": 2,
+            "t_bus": 3,
+            "br_r": 0.06744,
+            "br_x": 1.084,
+            "b_fr": 0.03059,
+            "b_to": 0.03059,
+            "br_status": 1,
+            "base_kv": 220,
+            "c_rating_a": 1.959,
+            "angmin": -0.17453292519943295,
+            "angmax": 0.17453292519943295,
+            "transformer": false,
+            "tap": 1,
+            "shift": 0,
+            "g_fr": 0,
+            "g_to": 0
+        }
+    },
+    "bus": {
+        "1": {
+            "bus_type": 1,
+            "vm": 1,
+            "va": 0,
+            "vmin": 0.5,
+            "vmax": 1.5,
+            "index": 1,
+            "gen_status": 1,
+            "base_kv": 400.0,
+            "bus_i": 1
+        },
+        "2": {
+            "bus_type": 2,
+            "vm": 1,
+            "va": 0,
+            "vmin": 0.5,
+            "vmax": 1.5,
+            "index": 2,
+            "gen_status": 1,
+            "base_kv": 400.0,
+            "bus_i": 2
+        },
+        "3": {
+            "bus_type": 3,
+            "vm": 1,
+            "va": 0,
+            "vmin": 0.5,
+            "vmax": 1.5,
+            "index": 3,
+            "gen_status": 1,
+            "base_kv": 220.0,
+            "bus_i": 3
+        }
+    }
+}

test/pf-native.jl

@@ -501,5 +501,11 @@ end
             @test isapprox(data["gen"][i]["qg"], gen["qg"]; atol = 1e-6)
         end
     end
+    @testset "test_issue_938" begin
+        filename = joinpath(@__DIR__, "data/json/issue_938.json")
+        data = PowerModels.parse_file(filename)
+        native = PowerModels.compute_ac_pf(data)
+        @test native["termination_status"]
+    end
 end