原理不解释,直接上代码
代码中被注释的源程序可用于打印中间结果,检查运算是否正确。
#include "mpi.h"#include#include #include #include void scatter_matrix(int* fstream,int n1,int n2,int*Q,int root,int tag){ /*每个矩阵块的大小*/ int rows=(n1+root-1)/root; int cols=(n2+root-1)/root; int* tmp_matrix=(int*)malloc(rows*cols*sizeof(int)); int i,j; memset(Q,0,rows*cols*sizeof(int)); for(i=0;i