cannon矩陣乘法mpi實(shí)現(xiàn)及性能的分析

《cannon矩陣乘法mpi實(shí)現(xiàn)及性能的分析》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在應(yīng)用文檔-天天文庫。

1、...頁眉#include#include#include#include#include#includeMPI_Statusstatus;double**A,**B,**C;//C=A*Bdouble*a,*b,*c;//各個(gè)進(jìn)程的緩沖區(qū)intn;//矩陣的行列數(shù)intnp;//每個(gè)進(jìn)程控制的小矩陣的行列數(shù)intp,rank;//進(jìn)程個(gè)個(gè)數(shù)、當(dāng)前進(jìn)程的編號,笛卡爾進(jìn)程編號double*tempa,*tempb;voidProduceABC();//在根處理器中生成矩陣AB,初始化矩陣Cvo

2、idPrintABC();//輸出結(jié)果voidScatterAB();//分發(fā)矩陣AB中的元素到各個(gè)進(jìn)程中voidMainProcess();//cannon算法的主過程voidcollectC();//收集結(jié)果矩陣CvoidMutiply();//矩陣相乘voidPrintab();voidPrintc();intmain(intargc,char*argv[]){inti;doublestarttime,endtime;MPI_Init(&argc,&argv);MPI_Comm_size(MPI_COMM_WORLD,&p);MPI_Comm_rank(MPI_COMM_WORLD,&r

3、ank);if(rank==0){....頁腳...頁眉printf("請輸入矩陣的行列數(shù)n=");fflush(stdout);scanf_s("%d",&n);printf("");}MPI_Bcast(&n,1,MPI_DOUBLE,0,MPI_COMM_WORLD);//n=atoi(argv[1]);np=n/(int)sqrt(p);a=(double*)malloc(np*np*sizeof(double));b=(double*)malloc(np*np*sizeof(double));c=(double*)malloc(np*np*sizeof(double));mems

4、et(c,0,np*np*sizeof(double));tempa=(double*)malloc(np*np*sizeof(double));tempb=(double*)malloc(np*np*sizeof(double));if(rank==0){//在根處理器中為矩陣ABC分配空間A=(double**)malloc(n*sizeof(double*));B=(double**)malloc(n*sizeof(double*));C=(double**)malloc(n*sizeof(double*));for(i=0;i

5、*sizeof(double));B[i]=(double*)malloc(n*sizeof(double));C[i]=(double*)malloc(n*sizeof(double));}ProduceABC();//在根處理器中隨機(jī)生成矩陣AB，初始化矩陣CScatterAB();//分發(fā)矩陣AB中的元素到各個(gè)進(jìn)程中}....頁腳...頁眉else{MPI_Recv(a,np*np,MPI_DOUBLE,0,1,MPI_COMM_WORLD,&status);MPI_Recv(b,np*np,MPI_DOUBLE,0,2,MPI_COMM_WORLD,&status);}starttim

6、e=MPI_Wtime();//開始時(shí)間MainProcess();//cannon算法的主過程if(rank==0){collectC();//收集結(jié)果矩陣CPrintABC();//輸出結(jié)果endtime=MPI_Wtime();printf("timeused:%lf",endtime-starttime);for(i=0;i

7、ee(b);free(c);....頁腳...頁眉free(tempa);free(tempb);MPI_Finalize();return0;}voidProduceABC()//在根處理器中生成矩陣AB{inti,j;for(i=0;i