Return perm matrix with LU factorization
Return perm matrix with #factorize_lu

renaming

refactor

refactor