Jacobi Iteration: C MPI Program Sketch
/* sizes of data and data files */
int xpoints = nx/npx; int ypoints = ny/npy;
double phi[ypoints+2][xpoints+2], oldphi[ypoints+2][xpoints+2];/* communication variables */int rank; int rankx, ranky; int coords[NDIM];int reorder = 0; int dims[NDIM], periods[NDIM];MPI_Comm comm2d;MPI_Datatype contig, strided;