矩阵相乘的并行算法的设计与实现.docVIP

仲恺农业工程学院实验报告纸 计算机科学与工程学院（院、系） 网络工程 专业 083 班 组 并行计算应用试验 课 学号：200810224311 姓名：李志冬 实验日期：2011-05-19 教师评定 实验三 矩阵相乘的并行算法的设计与实现 实验目的 理解和掌握矩阵相乘的并行算法的设计思想以及实现原理 实验内容 编译和运行一个两矩阵相乘算法的并行程序 实验步骤 1 使用vi编辑器输入并行计算的代码，保存在multi.c中 #include stdio.h #include mpi.h #define NRA 62 #define NCA 15 #define NCB 7 #define MASTER 0 #define FROM_MASTER 1 #define FROM_WORKER 2 MPI_Status status; int main(int argc, char *argv[]) { int numtasks, taskid, numworkers, source, dest, nbytes, mtype, intsize, dbsize, rows, averow,extra,offset, i,j,k, count; double a[NRA][NCA],b[NCA][NCB],c[NRA][NCB]; intsize = sizeof(int); dbsize = sizeof(double); MPI_Init(argc,argv); MPI_Comm_rank(MPI_COMM_WORLD,taskid); MPI_Comm_size(MPI_COMM_WORLD,numtasks); numworkers = numtasks-1; if(taskid==MASTER) { printf(Number of worker tasks = %d\n,numworkers); for(i=0;iNRA;i++) for(j=0;jNCA;j++) a[i][j]=i+j; for(i=0;iNCA;i++) for(j=0;jNCB;j++) b[i][j]=i*j; averow=NRA/numworkers; extra=NRA%numworkers; offset=0; mtype=FROM_MASTER; for(dest=1;dest=numworkers;dest++) { rows=(dest=extra)?averow+1:averow; printf(sending %d rows to task %d\n,rows,dest); MPI_Send(offset,1,MPI_INT,dest,mtype,MPI_COMM_WORLD); MPI_Send(rows,1,MPI_INT,dest,mtype,MPI_COMM_WORLD); count=rows*NCA; MPI_Send(a[offset][0],count,MPI_DOUBLE,dest,mtype,MPI_COMM_WORLD); count=NCA*NCB; MPI_Send(b,count,MPI_DOUBLE,dest,mtype,MPI_COMM_WORLD); offset=offset+rows; } mtype=FROM_WORKER; for(i=1;i=numworkers;i++) { source = i; MPI_Recv(offset,1,MPI_INT,source,mtype,MPI_COMM_WORLD,status); MPI_Recv(rows,1,MPI_INT,source,mtype,MPI_COMM_WORLD,status); count=rows*NCB; MPI_Recv(c[offset][0],count,MPI_DOUBLE,source,mtype,MPI_COMM_WORLD,status); } printf(Here is the result matrix\n); for(i=0;iNRA;i++) { printf(\n); for(j=0;jNCB;j++) printf(%6.2f ,c[i][j]); } printf(\n); } if(taskidMASTER) { mtype=FROM_MASTER;