程序 程序5.6.docVIP

程序5.6 /*文件名：mtpi.c*/ #includempi.h #include stdio.h #include stdlib.h int main(int argc,char **argv) { int myid, numprocs; int namelen,source; long count=1000000; char processor_name[MPI_MAX_PROCESSOR_NAME]; MPI_Status status; MPI_Init(argc,argv); MPI_Comm_rank(MPI_COMM_WORLD,myid);//得到当前进程的进程号 MPI_Comm_size(MPI_COMM_WORLD,numprocs);//得到通信域中的总进程数 MPI_Get_processor_name(processor_name,namelen);//得到节点主机名称 srand((int)time(0));//设置随机种子 double y; double x; long m=0,m1=0,i=0,p=0; double pi=0.0,n=0.0; for(i=0;icount;i++) { x=(double)rand()/(double)RAND_MAX;//得到0～1之间的随机数，x坐标 y=(double)rand()/(double)RAND_MAX;//得到0～1之间的随机数，y坐标 if((x-0.5)*(x-0.5)+(y-0.5)*(y-0.5)0.25)//判断产生的随机点坐标是否在圆内 m++; } n=4.0*m/1000000; printf(Process %d of %d on %s pi= %f\n,myid,numprocs,processor_name,n); if(myid!=0) //判断是否是主节点 { MPI_Send(m,1,MPI_DOUBLE,0,1,MPI_COMM_WORLD);//子节点向主节点传送结果 } else { p=m; /*分别接收来自于不同子节点的数据*/ for(source=1;sourcenumprocs;source++) { MPI_Recv(m1,1,MPI_DOUBLE,source,1,MPI_COMM_WORLD,status);//主节点接收数据 p+=m1; } printf(pi= %f\n,4.0*p/(count* numprocs));//汇总计算pi值 } MPI_Finalize(); } 《云计算与大数据技术》 并行计算实验室