Jacobi iterative C++ implementation