Basic LUP factorization