bugfix shape penalty factors

NREL-Sienna · Jan 17, 2025 · 2106c0c · 2106c0c
1 parent d2fb7b6
commit 2106c0c
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/src/newton_ac_powerflow.jl b/src/newton_ac_powerflow.jl
@@ -724,7 +724,7 @@ function _newton_powerflow(
     else
         solver_data.J = J
         solver_data.dSbus_dV_ref =
-            [vec(real.(dSbus_dVa[ref, :][:, pvpq])); vec(real.(dSbus_dVm[ref, :][:, pvpq]))]
+            [vec(real.(dSbus_dVa[ref, :][:, pvpq])); vec(real.(dSbus_dVm[ref, :][:, pq]))]
         @info("The powerflow solver with KLU converged after $i iterations")
     end
 

diff --git a/test/test_utils/legacy_pf.jl b/test/test_utils/legacy_pf.jl
@@ -118,7 +118,7 @@ function _newton_powerflow(
         Sbus_result = V .* conj(Ybus * V)
         solver_data.J = J
         solver_data.dSbus_dV_ref =
-            [vec(real.(dSbus_dVa[ref, :][:, pvpq])); vec(real.(dSbus_dVm[ref, :][:, pvpq]))]
+            [vec(real.(dSbus_dVa[ref, :][:, pvpq])); vec(real.(dSbus_dVm[ref, :][:, pq]))]
         @info("The powerflow solver with KLU converged after $i iterations")
     end
     return (converged, V, Sbus_result)