Function gramschmidt::cgs2

pub fn cgs2<S>(
    a: &ArrayBase<S, Ix2>
) -> Result<(Array<f64, Ix2>, Array<f64, Ix2>)> 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.