Function qr 

pub fn qr<T>(input: &Array<T>) -> (Array<T>, Array<T>, Array<T>)
where
    T: HasAfEnum + FloatingPoint,

Perform QR decomposition

Parameters

• input is the input matrix

Return Values

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