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