矩陣相乘的并行算法的設(shè)計(jì)與實(shí)現(xiàn).doc

《矩陣相乘的并行算法的設(shè)計(jì)與實(shí)現(xiàn).doc》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫(kù)。

1、仲愷農(nóng)業(yè)工程學(xué)院實(shí)驗(yàn)報(bào)告紙計(jì)算機(jī)科學(xué)與工程學(xué)院（院、系）網(wǎng)絡(luò)工程專業(yè)083班組并行計(jì)算應(yīng)用試驗(yàn)課學(xué)號(hào)：200810224311姓名：李志冬實(shí)驗(yàn)日期：2011-05-19教師評(píng)定實(shí)驗(yàn)三矩陣相乘的并行算法的設(shè)計(jì)與實(shí)現(xiàn)一、實(shí)驗(yàn)?zāi)康睦斫夂驼莆站仃囅喑说牟⑿兴惴ǖ脑O(shè)計(jì)思想以及實(shí)現(xiàn)原理二、實(shí)驗(yàn)內(nèi)容編譯和運(yùn)行一個(gè)兩矩陣相乘算法的并行程序三、實(shí)驗(yàn)步驟1使用vi編輯器輸入并行計(jì)算的代碼，保存在multi.c中#include#include"mpi.h"#defineNRA62#defineNCA15#defineNCB7#defineMASTER0#defineFROM_MA

2、STER1#defineFROM_WORKER2MPI_Statusstatus;intmain(intargc,char*argv[]){intnumtasks,taskid,numworkers,source,dest,nbytes,mtype,intsize,dbsize,rows,averow,extra,offset,i,j,k,count;doublea[NRA][NCA],b[NCA][NCB],c[NRA][NCB];intsize=sizeof(int);第7頁共7頁dbsize=sizeof(double);MPI_Init(&argc,&argv);MPI

3、_Comm_rank(MPI_COMM_WORLD,&taskid);MPI_Comm_size(MPI_COMM_WORLD,&numtasks);numworkers=numtasks-1;if(taskid==MASTER){printf("Numberofworkertasks=%d",numworkers);for(i=0;i

4、mworkers;offset=0;mtype=FROM_MASTER;for(dest=1;dest<=numworkers;dest++){rows=(dest<=extra)?averow+1:averow;printf("sending%drowstotask%d",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][

5、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

6、,MPI_COMM_WORLD,&status);count=rows*NCB;MPI_Recv(&c[offset][0],count,MPI_DOUBLE,source,mtype,MPI_COMM_WORLD,&status);}printf("Hereistheresultmatrix");for(i=0;iMASTER){mtype=FROM_MASTER;source=MA

7、STER;printf("Master=%d,mtype=%d",source,mtype);MPI_Recv(&offset,1,MPI_INT,source,mtype,MPI_COMM_WORLD,&status);printf("offset=%d",offset);MPI_Recv(&rows,1,MPI_INT,source,mtype,MPI_COMM_WORLD,&status);printf("rows=%d",rows);count=rows*NCA;MPI_