pub fn qr<T>(input: &Array<T>) -> (Array<T>, Array<T>, Array<T>) where
T: HasAfEnum + FloatingPoint,
Perform QR decomposition
inputis the input matrix
A triplet of Arrays.
The first Array is the orthogonal matrix from QR decomposition
The second Array is the upper triangular matrix from QR decomposition
The third Array will contain additional information needed for solving a least squares problem using q and r