I am inverting a sparse, symmetric, ill-conditioned matrix. I have used both SVD and the LDL decomposition. I find that my results are better with the latter. Why?
I understand that LDL decomposition works only on symmetric matrices. But that doesn't necessarily say why it works better than a "more general" algorithm.