Fixed a bug in LDU solvers; factorization is symmetric for a symmetric input matrix but need to check that it actually works