def calc_power(V, Ybus, Ibus):
"""
:param V:
:param Ybus:
:param Ibus:
:return:
"""
# S = V * np.conj(Ybus * V - Ibus)
# Y = Ybus.tocsr()
S = calc_power_csr_numba(n=V.shape[0], Yp=Ybus.indptr, Yj=Ybus.indices, Yx=Ybus.data, V=V, I=Ibus, n_par=500)
return S
def calc_power(V, Ybus, Ibus):
"""
Compute the power from a voltage solution
:param V: Voltages Vector
:param Ybus: Admittance Matrix
:param Ibus: Currents vector
:return: Power injections
"""
# S = V * np.conj(Ybus * V - Ibus)
# Y = Ybus.tocsr()
# return S
return calc_power_csr_numba(n=V.shape[0],
Yp=Ybus.indptr,
Yj=Ybus.indices,
Yx=Ybus.data,
V=V,
I=Ibus,
n_par=500)