Comment by Porygon

Conjugate gradient descent with a multigrid preconditioner also works quite well in my experience, especially for larger systems.