MPI編程環(huán)境配置與示例

1、.目錄一.系統(tǒng)安裝21.下載地點22.安裝步驟23.VisualStdio設置2二.實驗程序31.簡單的MPI編程示例32.信息傳達MPI編程示例143.信息傳達MPI編程示例254.MonteCarlo方法計算圓周率75.計算積分9三.心得意會10.一.系統(tǒng)安裝下載地點FTP匿名登岸，在pub/mpi/nt文件夾中安裝步驟在安裝有MPI的計算機上要成立一個有管理員權限的賬戶，不能夠夠沒有密碼；雙擊exe文件，按默認設置安裝MPI；注冊MPI賬戶，調用MPIRegister.exe，用戶名和密碼即為第一步中的賬戶。VisualStdio設置為預防每新建一個項目都要設置一次，能夠對它進行通用設置

2、。打開視圖-其他窗口-屬性管理器，點擊Debug|Win32目錄下的Microsoft.Cpp.Win32.user，在VC+目錄下的包含目錄中增添MPICH的Include路徑，庫目錄中增添MPICH的Lib路徑；在鏈接器-輸入目錄下的附加依靠項中增添mpich.lib、mpe.lib、mped.lib、mpichd.lib。.二.實驗程序簡單的MPI編程示例源代碼#include#includeintmain(intargc,char*argv)intnum,rk;MPI_Init(&argc,&argv);MPI_Comm_size(MPI_COMM_WORLD,&num);MPI_Co

3、mm_rank(MPI_COMM_WORLD,&rk);printf(HelloworldfromProcess%dof%dn,rk,num);MPI_Finalize();return0;運行截圖.信息傳達MPI編程示例1源代碼#include#includeintmain(intargc,charargv)intmyid,numprocs,source;MPI_Statusstatus;charmessages100;MPI_Init(&argc,&argv);MPI_Comm_size(MPI_COMM_WORLD,&numprocs);MPI_Comm_rank(MPI_COMM_WO

4、RLD,&myid);if(myid!=0)strcpy(messages,HelloWorld!);MPI_Send(messages,strlen(messages)+1,MPI_CHAR,0,99,MPI_COMM_WORLD);elsefor(source=1;sourcenumprocs;source+).MPI_Recv(messages,100,MPI_CHAR,source,99,MPI_COMM_WORLD,&status);printf(Iamprocess%dIrecvstring%sfromprocess%dn,myid,messages,source);MPI_Fin

5、alize();return0;運行截圖信息傳達MPI編程示例2源代碼#include#include#include#defineN1002intmain(intargc,charargv).intmyid,numprocs,C=0;doubledataN,SqrtSum=0.0;_doubled;for(inti=0;iN;i+)datai=i;MPI_Statusstatus;MPI_Init(&argc,&argv);_MPI_Comm_rank(MPI_COMM_WORLD,&myid);MPI_Comm_size(MPI_COMM_WORLD,&numprocs);-numproc

6、s;if(myid=0)for(inti=0;iN;+i)MPI_Send(&datai,1,MPI_DOUBLE,i%numprocs+1,1,MPI_COMM_WORLD);for(intsource=1;source=numprocs;+source)MPI_Recv(&d,1,MPI_DOUBLE,source,99,MPI_COMM_WORLD,&status);SqrtSum+=d;elsefor(inti=myid;iN;i=i+numprocs)MPI_Recv(&d,1,MPI_DOUBLE,0,1,MPI_COMM_WORLD,&status);SqrtSum+=sqrt(

7、d);+C;MPI_Send(&SqrtSum,1,MPI_DOUBLE,0,99,MPI_COMM_WORLD);printf(Iamprocess%d.Irecvtotal%dfromprocess0,andSqrtSum=%f.n,myid,C,SqrtSum);_MPI_Finalize();_return0;.運行截圖MonteCarlo方法計算圓周率源代碼#includempi.h#include#includeintmain(intargc,charargv)intmyid,numprocs;intsource;longcount=1000000;MPI_Statusstatus

8、;MPI_Init(&argc,&argv);MPI_Comm_rank(MPI_COMM_WORLD,&myid);MPI_Comm_size(MPI_COMM_WORLD,&numprocs);srand(myid);doubley,x,pi=0.0,n=count;longm=0,m1=0,i=0,p=0;for(i=0;icount;i+).x=(double)rand()/(double)RAND_MAX;y=(double)rand()/(double)RAND_MAX;if(x-0.5)*(x-0.5)+(y-0.5)*(y-0.5)0.25)+m;pi=4.0*m/n;prin

9、tf(Process%dof%pi=%fn,myid,numprocs,pi);if(myid!=0)MPI_Send(&m,1,MPI_DOUBLE,0,1,MPI_COMM_WORLD);elsep=m;for(source=1;sourcenumprocs;source+)MPI_Recv(&m1,1,MPI_DOUBLE,source,1,MPI_COMM_WORLD,&status);p=p+m1;printf(pi=%fn,4.0*p/(n*numprocs);MPI_Finalize();return0;運行截圖.計算積分源代碼#defineN100000000#definea0

10、#defineb10#include#include#include#includempi.hintmain(intargc,charargv)intmyid,numprocs;inti;_doublelocal=0.0,dx=(double)b-(double)a)/(double)N;doubleinte,x;MPI_Init(&argc,&argv);_MPI_Comm_rank(MPI_COMM_WORLD,&myid);_MPI_Comm_size(MPI_COMM_WORLD,&numprocs);for(i=myid;iN;i=i+numprocs).x=a+i*dx+dx/2;local+=x*x*dx;MPI_Reduce(&local,&inte,1,MPI_DOUBLE,MPI_SUM,0,MPI_COMM_WORLD);_if(myid=0)printf(Theintegalofx*xinregion%d,%d=%16.15fn,a,b,inte);MPI_F