pub fn invert_spd(matrix: &ArrayView2<'_, f64>) -> Result<Array2<f64>, Error>