Function lu_transpose 

pub fn lu_transpose<T: RlstScalar, Array2MutImpl: MutableArrayImpl<T, 2>>(
    mat: &mut Array<Array2MutImpl, 2>,
) -> Vec<usize>

Compute the LU decomposition of the transpose of a square matrix