Function spd_inverse 

fn spd_inverse(a: &[f64], n: usize) -> Vec<f64>

SPD inverse via Cholesky (needed for the explicit nf · M⁻¹ term).