Jacobi and Gauss-Seidel solvers