Function Matrix_cannon_scatter

code/Chapter11/mympi.h:46–62

}

void Matrix_cannon_scatter(int N,int id,int num,int part,data_t *source_matrix,data_t *target_matrix,MPI_Comm MPI_COMM_CART){
    data_t *matrixpool;//用于分发矩阵的一维向量
    int i,j,w,h,count=0;
    if(id==0){//(0,0)进程
		matrixpool=malloc(sizeof(data_t)*N*N);
        //将矩阵source_matrix重新排序放入matrixpool
        for(w=0;w<part;w++)
            for(h=0;h<part;h++)
                for(i=w*num;i<(w+1)*num;i++)
                    for(j=h*num;j<(h+1)*num;j++)
                        matrixpool[count++] = source_matrix[i*N+j]; 
    }
    MPI_Scatter(matrixpool,num*num,MPI_FLOAT,target_matrix,num*num,MPI_FLOAT,0,MPI_COMM_CART);
    if(id==0){
        free(matrixpool);
	}
}
void Matrix_cannon_gather(int N,int id,int num,int part,data_t *source_c,data_t *target_c,MPI_Comm MPI_COMM_CART){
    data_t *matrixpool;//用于收集C矩阵的一维向量
    int i,j,w,h,posi;