Perm matrix with LU fact

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

renaming

refactor

refactor

docs