[Porygon]

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

[slavik81]

I've never worked with multi-grid, but I assume the preconditioner is also based on the Cholesky factorization? Incomplete Cholesky was pretty effective as the preconditioner for the pressure solve in my toy fluid sim.