基于Cannon方法的矩阵相乘mpi编程.docxVIP

基于Cannon方法的矩阵相乘mpi编程 问题描述： 矩阵A, B 均为N*N的方阵，试计算矩阵C=A*B；使用P个进程并行计算（N可以被P整除）；矩阵A，B及C均采用分布式存储，采用Cannon、Fox、Systolic或DNS算法；编写计算C矩阵的MPI程序，并进行计算。实际计算时， 矩阵A, B请采用如下值： 要求：计算出C矩阵后，请计算，并由根节点打印出来。 将S值与串行程序的结果进行对比，校验程序的正确性； 分别N值设为100,200,500，并分别用1,2,4,10个进程进行计算，并利用MPI_Wtime( )函数计算程序的运行时间，考核不同N值和不同进程数下的程序可观测加速比。 Cannon方法原理说明 图 SEQ 图 \* ARABIC 1 分块与进程Cannon方法要求进程数p是完全平方数。如下图，先将A、B、C分成的方阵块，，，对应储存在进程（i，j=0，1，……）中。然后按下面的算法原理写程序即可。 图 SEQ 图 \* ARABIC 1 分块与进程 Cannon算法原理很简单，主要分四步： 所有块Aij向左循环移动i步；所有块Bij向上循环移动j步； 所有处理器Pij执行Aij和Bij的乘加运算； 所有Aij向左循环移动1步；所有Bij向上循环移动1步； 转②执行 次。 图 SEQ 图 \* ARABIC 图 SEQ 图 \* ARABIC 2 执行步骤① 图 图 SEQ 图 \* ARABIC 3 执行步骤③和④ 图 SEQ 图 \* ARABIC 4 执行步骤② 详细算法如下： 输入： , 输出： Begin 对所有处理器myid(myid=0,……p-1)同时执行如下算法： 计算子块行号i=myid/m，j=myid%m;（其中m=sqrt(p)）； for k=0 to m-1 do if(ik)then LeftMoveOneStep(); if(jk)then UpMoveOneStep(); end for for k=0 to L-1 do //(L=N/m为子块维度) for l=0 to L-1 do C[k*L+j]=0 end for end for for s=0 to m-1 do for k=0 to L-1 do for l=0 to L-1 do for t=0 to L-1 do C[k*L+l]+=A[k*L+t]*B[t*L+l]; end for end for end for LeftMoveOneStep(); UpMoveOneStep(); end for End 函数LeftMoveOneStep()算法： Begin ① 将所存的A子块发给Id_Move(i,j-1,m)进程的bufferA； ② 从Id_Move(I,j+1,m)进程接收bufferA，复制到本进程的A。 End 函数UpMoveOneStep()算法： Begin ①将所存的B子块发给Id_Move(i-1,j,m)进程的bufferB； ②从Id_Move(i+1,j,m)进程接收bufferB，复制到本进程的B。 End 函数Id_Move(i,j,m)算法： return [(i+m)%m]*m+(j+m)%m即可。 程序源代码 先给出串行代码，增强对矩阵相乘的理解，从而降低编程难度；不仅可用于检验并行程序结果，还可比较并行效率。 //运行模式为./name 500.500为矩阵大小。 #includestdio.h #includestdlib.h #includemath.h #includetime.h int main(int argc,char *argv[]) { int i,j,k,N; clock_t tstart,tend; double *x,*y,**A,**B,**C,C2_sum=0,ttotal,*tmp1,*tmp2,*tmp3,*tmp4;//定义tmp1-4是因为初始化A和B时，tmp们只需计算一次。 tstart=clock(); N=atoi(argv[1]); x=(double *)malloc(N*sizeof(double)); y=(double *)malloc(N*sizeof(double)); A=(double **)malloc(N*sizeof(double **)*N); B=(double **)malloc(N*sizeof(double **)*N); C=(double **)malloc(N*sizeof(double **)*N); tmp1=(double *)malloc(N*sizeof(double));