Function calc_susceptance_diagonal 

pub fn calc_susceptance_diagonal(b: &[f64]) -> CsMat<f64>

B = diag(b), shape m × m.