2009并计算-5.pptVIP

并行计算 ——mpi程序设计 程序2 π的计算 程序2 π的计算 #include “mpi.h” #include stdio.h #include math.h #define N 1000 double f(double x);/*定义函数f(x) */ {return(4.0/(1.0+x*x)); } int main (int argc,char * argv[]) {int done =0,n,myid,numprocs,i;double PI25DT=3.62643;double mypi,pi,h,sum,x;double startwtime=0.0,endwtime;int namelen;char processor_name[MPI_MAXPROCESSOR_NAME];MPI_Status status;MPI_Init(argc,argv);MPI_Comm_size(MPI_COMM_WORLD,numprocs);MPI_Comm_rank(MPI_COMM_WORLD,myid);printf(”Process %d of %d\n”,myid,numprocs);n=N;if (myid==0){startwtime=MPI_Wtime();for (j=1;jnumprocs;j++){MPI_Send(n,1,MPI_INT,j,99,MPI_COMM_WORLD);}}else { MPI_Recv(n,1,MPI_INT,MPI_ANY_SOURCE,99,MPI_COMM_WORLD,status); }h=1.0/(double) n;sum=0.0;for(i=myid+1;i=n;i+=numprocs) /* 每一个进程计算一部分矩形的面积，若进程总数numprocs为4，将 0-1区间划分为100个矩形，则各个进程分别计算矩形块0进程 1，5，9，13，……，971进程 2，6，10，14，……，982进程 3，7，11，15，……，993进程 4，8，12，16，……，100 */{x=h*((double)i-0.5);sum+=f(x);}mypi=h*sum; /*各进程并行计算得到的部分和*//*将部分和累加得到所有矩形的面积，该面积和即为近似PI值*/if (myid != 0)MPI_Send(mypi,1,MPI_DOUBLE,0,myid,MPI_COMM_WORLD) ；else{pi=0.0; pi=pi+mypi;for (j=1;jnumprocs;j++){MPI_Recv(mypi,1,MPI_DOUBLE,MPI_ANY_SOURCE,MPI_ANY_TAG,MPI_COMM_WORLD,status);pi=pi+mypi;}printf(“pi is approximately %.16f,Error is %.16f\n”,pi,fabs(pi-PI25DT));endwtime=MPI_Wtime();printf(“wall clock time=% f\n”,endwtime-startwtime);fflush(stdout);}MPI_Finalize(); } 5. MPI群体通信 群体通信：在一个通信体的所有进程间同时进行地通信。 群体通信包括： 广播（MPI_Bcast） 数据收集（MPI_Gather） 数据散发（MPI_Scatter） 数据转置（MPI_Alltoall） 归约（MPI_Reduce） a）广播 (MPI_Bcast) ： MPI_Bcast ( void *buffer, int count, MPI_Datatype datatype, int root, MPI_Comm comm ) 所有进程都要执行此函数，只是指定的root进程进行广播操作。 将一条消息广播发送到进程组内的所有进程，包括root自己。 此函数没有对应的receive函数，不需要进行receive。 所有进程的缓冲元素数目必须一样大。 b）收集(MPI_Gather) ： int MPI_Gather (void *sendbuf,int sendcount,MPI_Datatype sendtype,void *recvbuf,int recvcount,MPI_Datatype recvtype,int root,MPI_Comm comm) c）散发 (MPI_Scatter) ： MPI_Scatter( void* sendbuf, int sendcount, MPI_Datatype sendtype,