LDU factorization