MPI程序 对等模式的 Jacobi迭代

/*
   对等模式的MPI程序：test_8_1_2.c
   Jacobi迭代：迭代数据按列进行分割，
   并假设一共有4个进程同时并行计算。



时间:15.7.27
Jason Zhou   
热爱你所写下的程序，他是你的伙伴，而不是工具.
*/
#include "mpi.h"
#include <stdio.h>

#define totalsize 8
#define mysize totalsize/4    //分成四块，每块大小
#define steps 10

int  main(int argc,char* argv[])
{
        int myid,numprocs,n,i,j,rc;
        float a[totalsize][mysize+2],b[totalsize][mysize+2],tmp[mysize][totalsize],c[totalsize][totalsize];  //除分块大小外，还包括左右两边各一列
        float temp[totalsize];  /* 临时数组 */
        int begin_col,end_col,ierr;
        MPI_Status status;


        //初始化c数组
        for(i=0;i<totalsize;i++)
        {
                for(j=0;j<=totalsize;j++)
                        c[i][j]=0;
        }

        for(int j=0;j<mysize;j++)
                for(int i=0;i<totalsize;i++)
                        tmp[j][i]=8;


        MPI_Init(&argc,&argv);

        /* 得到当前进程标识和总的进程个数 */
        MPI_Comm_rank(MPI_COMM_WORLD,&am