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

[DOC] -矩阵相乘的并行算法的设计与实现 矩阵相乘的并行算法的设计与实现 仲恺农业工程学院实验报告纸 计算机科学与工程学院(院、系)网络工程专业083班组并行计算 应用试验课 实验三矩阵相乘的并行算法的设计与实现 一、 实验目的 理解和掌握矩阵相乘的并行算法的设讣思想以及实现原理 二、 实验内容 编译和运行一个两矩阵相乘算法的并行程序 三、 实验步骤 1使用vi编辑器输入并行计算的代码，保存在multi, c中 #include stdio. h #include mpi. h ttdefine NRA 62 #define NCA 15 #define NCB 7 ^define MASTER 0 itdefine FROM_MASTER 1 itdefine 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): 第1页共7页 MPI_Init(Aargc, Aargv); MPI_Comm_rank(MPI_C0MM_W0RLD, Ataskid); MPI_Comm_size(MPI_COMM_WORLD, numtasks): numworkers = numtasks-1; 辻(taskid二二MASTER) { printf (Number of worker tasks =%d\n，z, 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 (destHT dest〈Hnunworkers-dest++)宀 rows(dest〈Hextra) 7averow+ldverow- printf (sending %d rows Co task %dxn、「rows-dest)- MPIISend (冷offset. MPIIINT? dest? nltypo MPIICOMMIWORLD)- MPIISend(rowyrMPIIINL dest》mtype》MPIICOMMIWORLD) - countrows*NC MPIISend(々a〔offset〕o〕9 counrMPIIDOUBL0 desrIntypo MPIICOMMIWORLD)八 countNCA*NCB- MPIISendKFcounty MPIIDOUBLEy desr-nltypo MPIICOMMIWORLD)- offsetuoffset+rows- ) IntypeHFROMlwoRKE?0for (iHl 二〈Hnumworkers- i++) (source H T MPIIRecvKoffserrMPIlINTy source》mtypo MPIICOMMIWORLP^status)- MPIIRecvQrowyrMPIIINL source》mtypp MPIICOMMIWORLP々status)- counturows*NCB- MPIIRecv (冷 c〔off set〕o〕9 count》MPIID0UBLE9 source》Intype》MPIICOMMIWORLP □status: ) printf(、xere is the result matrixxn、、： for (ipi〈NRA=++) (printf (、、\n、、?rfor(ll-oL.〈NCBL.++) printf (、、％6? 2f、「c「r-.c?〕)【) printf(、、\n、、1 )黑 2 刃艸 7 刃 if(taskid〉MASTER)( IntypeHFROMIMASTER - sourceHMASTER八 printf (vasterld mtype=%dxn、、9 sourc 卩 mtype)- MPIIRecvQoffserrMPIIINT9 source》mtypp MPIICOM