Regular Jacobi solver