Function connection_power_map 

pub fn connection_power_map(
    incidence: &[Vec<f64>],
    reference: &[Complex64],
) -> Result<ConnectionPowerMap>

Construct H = diag(vbar) D' diag(D vbar)^-1.

A channel power s_channel maps to terminal powers as s_terminal = H s_channel; the split is frozen at the reference phasors.

Errors

The incidence matrix is empty, ragged, non-finite, has the wrong terminal arity, or produces a zero reference winding voltage.