def _run_ac_pf_without_qlims_enforced(ppci, options):
if options["numba"]:
makeYbus = makeYbus_numba
pfsoln = pfsoln_numba
else:
makeYbus = makeYbus_pypower
pfsoln = pfsoln_pypower
baseMVA, bus, gen, branch, ref, pv, pq, _, _, V0 = _get_pf_variables_from_ppci(
ppci)
ppci, Ybus, Yf, Yt = _get_Y_bus(ppci, options, makeYbus, baseMVA, bus,
branch)
## compute complex bus power injections [generation - load]
Sbus = makeSbus(baseMVA, bus, gen)
## run the newton power flow
V, success, _, ppci["internal"]["J"] = newtonpf(Ybus, Sbus, V0, pv, pq,
ppci, options)
## update data matrices with solution
bus, gen, branch = pfsoln(baseMVA, bus, gen, branch, Ybus, Yf, Yt, V, ref)
return ppci, success, bus, gen, branch
def _run_ac_pf_without_qlims_enforced(ppci, options):
if options["numba"]:
try:
makeYbus = makeYbus_numba
except:
makeYbus = makeYbus_pypower
else:
makeYbus = makeYbus_pypower
baseMVA, bus, gen, branch, ref, pv, pq, _, gbus, V0 = _get_pf_variables_from_ppci(
ppci)
ppci, Ybus, Yf, Yt = _get_Y_bus(ppci, options, makeYbus, baseMVA, bus,
branch)
## compute complex bus power injections [generation - load]
Sbus = makeSbus(baseMVA, bus, gen)
## compute complex bus current injections from constant current loads
Ibus = _get_ibus(ppci)
## run the newton power flow
V, success, _ = newtonpf(Ybus, Sbus, V0, pv, pq, options, Ibus=Ibus)
## update data matrices with solution
bus, gen, branch = pfsoln(baseMVA,
bus,
gen,
branch,
Ybus,
Yf,
Yt,
V,
ref,
pv,
pq,
Ibus=Ibus)
return ppci, success, bus, gen, branch