Function rusfun::utils::LU_decomp

pub fn LU_decomp<T>(A: &Array2<T>) -> (Array2<T>, Array2<T>, Array2<T>) where
    T: NdFloat, 

Performs partial pivoted LU decomposition of A such that P A = LU with L a lower triangular matrix and U an upper triangular matrix A needs to be a square matrix