Function lu 

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

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.