Function qr

src/lapack/mod.rs:208–225

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

(input: &Array<T>)

/// The third Array will contain additional information needed for solving a least squares problem
/// using q and r
pub fn qr<T>(input: &Array<T>) -> (Array<T>, Array<T>, Array<T>)
where
    T: HasAfEnum + FloatingPoint,
{
    unsafe {
        let mut q: af_array = std::ptr::null_mut();
        let mut r: af_array = std::ptr::null_mut();
        let mut tau: af_array = std::ptr::null_mut();
        let err_val = af_qr(
            &mut q as *mut af_array,
            &mut r as *mut af_array,
            &mut tau as *mut af_array,
            input.get(),
        );
        HANDLE_ERROR(AfError::from(err_val));
        (q.into(), r.into(), tau.into())
    }
}

/// Perform inplace QR decomposition
///