pub fn qr_ndarray<T>(a: &Array2<T>) -> LapackResult<QrResult<T>>where T: Zeroable + Real + Field + Clone,
Computes the QR decomposition of a matrix.
A = Q * R
a
QR decomposition with Q and R