并行计算矩阵分块乘法.docxVIP

一、题目及要求1、题目简单并行分块乘法：（1）情形1: 超立方连接；（2）情形2：二维环绕网孔连接已知求。2、要求（1）题目分析、查阅与题目相关的资料；（2）设计算法；（3）算法实现、代码编写；（4）结果分析和性能分析与改进；（5）论文撰写、答辩；二、设计算法、算法原理要考虑的计算问题是C=AB,其中A与B分别是矩阵。①A、B和C分成的方块阵,和,大小均为 ，p个处理器编号为, 存放,和。②通讯:每行处理器进行A矩阵块的多到多播送(得到, k=0~ )每列处理器进行B矩阵块的多到多播送(得到, k=0~ )③乘-加运算: 做三、算法描述、设计流程3.1算法描述超立方情形下矩阵的简单并行分块算法输入：待选路的信包在源处理器中输出：将原处理器中的信包送至其目的地Beginfori=1 to n doendfor (2) (3) whiledo (3.1)ifthen从当前节点V选路到节点为V1 （3.2）endwhileEnd二维网孔情形下矩阵的简单并行分块算法输入：待选路的信包处于源处理器中输出：将各信包送至各自的目的地中Begin沿x维将信包向左或向右选路至目的地的处理器所在的列沿y维将信包向上或向下选路至目的地的处理器所在的行分块乘法算法 //输入: , ; 子快大小均为 输出: nBegin (1)fori=0 to dofor all par-doifik thenendifif jk thenB(i+1)mod , jendifendforendforfori=0 to dofor allpar-do=+endforEndforEnd3.2设计流程以下是二维网孔与超立方连接设计流程。 如图3-1二维网孔步骤：(1)先进行行播送;(2)再同时进行列播送;图3-1 二维网孔示意图超立方步骤：依次从低维到高维播送, d-立方, d=0,1,2,3,4…;算法流程如图所示：图3-2 算法流程四、源程序代码及运行结果1、超立方1.1超立方的源程序代码#include stdio.h#include stdlib.h#include mpi.h#define intsizesizeof(int)#define floatsizesizeof(float)#define charsizesizeof(char)#define A(x,y) A[x*K+y]#define B(x,y) B[x*N+y]#define C(x,y) C[x*N+y]#define a(x,y) a[x*K+y]#define b(x,y) b[x*n+y]#define buffer(x,y) buffer[x*n+y] #define c(l,x,y) c[x*N+y+l*n]float *a,*b,*c,*buffer;int s;float *A,*B,*C; int M,N,K,P ;intm,n;intmyid;int p; FILE *dataFile; MPI_Status status;double time1;doublestarttime,endtime;voidreadData(){inti,j;starttime = MPI_Wtime();dataFile=fopen(yin.txt,r);fscanf(dataFile,%d%d, M, K); A=(float *)malloc(floatsize*M*K); for(i = 0; i M; i++) {for(j = 0; j K; j++) {fscanf(dataFile,%f, A+i*K+j); } }fscanf(dataFile,%d%d, P, N); if (K!=P) {printf(the input is wrong\n);exit(1); } B=(float *)malloc(floatsize*K*N); for(i = 0; i K; i++) {for(j = 0; j N; j++) {fscanf(dataFile,%f, B+i*N+j); } }fclose(dataFile);printf(Input of