\bm{v}^{(m+1)} = R \bm{v}^{(m)} - B \bm{f}$$ We also saw that the convergence of the iterative calculation is dominated by the spectral radius $${ \rho(R) }$$ of the iteration matrix $${ R }$$. In ...
double* matrix_vec_prod(double* A, double* vec, int n, double* result_vec, int rank, int size, MPI_Comm comm); void gauss_jacobi(double* A, double* rhs_vec, double* initial_guess, int n, int max_iter, ...