Function lu

src/lapack/mod.rs:150–167

Perform LU decomposition # Parameters - `input` is the input matrix # Return Values A triplet of Arrays. The first Array will contain the lower triangular matrix of the LU decomposition. The second Array will contain the lower triangular matrix of the LU decomposition. The third Array will contain the permutation indices to map the input to the decomposition.

(input: &Array<T>)

///
/// The third Array will contain the permutation indices to map the input to the decomposition.
pub fn lu<T>(input: &Array<T>) -> (Array<T>, Array<T>, Array<i32>)
where
    T: HasAfEnum + FloatingPoint,
{
    unsafe {
        let mut lower: af_array = std::ptr::null_mut();
        let mut upper: af_array = std::ptr::null_mut();
        let mut pivot: af_array = std::ptr::null_mut();
        let err_val = af_lu(
            &mut lower as *mut af_array,
            &mut upper as *mut af_array,
            &mut pivot as *mut af_array,
            input.get(),
        );
        HANDLE_ERROR(AfError::from(err_val));
        (lower.into(), upper.into(), pivot.into())
    }
}

/// Perform inplace LU decomposition
///