add LU-decompose, LDU-decompose