2

Jacobi or Gauss-Seidel are not really state of the art for solving systems of linear equations. It is more done by preconditioned conjugated gradient (for symmetric positive semi-definite matrices) and preconditioned (F)GMRES (or other Krylov subspace methods) for arbitrary matrices. The crucial part here is the preconditioner. There was recently (=21th ...


Only top voted, non community-wiki answers of a minimum length are eligible