Bug in the weeds of _reactive_power_redistribution_pv

[GabrielKS]

The helper function _reactive_power_redistribution_pv:

PowerFlows.jl/src/post_processing.jl

Line 283 in 558a337

function _reactive_power_redistribution_pv(sys::PSY.System, Q_gen::Float64, bus::PSY.Bus)

takes a system, a reactive power, and a bus; and distributes that reactive power among the devices that contribute to the bus. We would therefore expect that after this function is run, the sum of the reactive power among the relevant devices is equal to the Q_gen provided to the function. However:

it isn't.

[GabrielKS]

@jd-lara could you verify that this is in fact a bug? It should definitely be possible to distribute that amount of reactive power because that amount of reactive power is already distributed.

[GabrielKS]

I have a test here that would fail if I did not manually exclude reactive power from the comparison. The fact that it works for active power (and voltage magnitude and voltage angle) suggests that this is indeed a bug.