Function gramschmidt::cgs2

pub fn cgs2<S>(a: &ArrayBase<S, Ix2>) -> Result<(Array2<f64>, Array2<f64>)> where
    S: Data<Elem = f64>, 

Convenience function that calculates a Reorthogonalized Gram Schmmidt QR factorization (see Giraud et al. for details), returning a tuple (Q,R).

If you want to repeatedly calculate QR factorizations, then prefer constructing a Reorthogonalized struct and calling its GramSchmidt::compute method implemented through the GramSchmidt trait.