程序 程序5.8.docVIP

程序5.8 /*文件名 myreduce.c*/ #define N 100000000 #include stdio.h #include stdlib.h #include time.h #include mpi.h void Myreduce(double *sendbuf, double *recvbuf,int count,int root);//定义自己的reduce函数 int main(int argc, char** argv) { int myid,numprocs; int i; double local=0.0; double inte,tmp=0.0,x; MPI_Init(argc, argv); MPI_Comm_rank(MPI_COMM_WORLD, myid); MPI_Comm_size(MPI_COMM_WORLD,numprocs); /*采用归约对y=x*x在[1,10]区间求积分*/ srand((int)time(0)); for(i=myid;iN;i=i+numprocs) { x=10.0*rand()/(RAND_MAX+1.0); tmp=x*x/N; local=tmp+local; } Myreduce(local,inte,1,0);//调用自定义的规约函数 if(myid==0) { printf(The integal of x*x=%16.15f\n,inte); } MPI_Finalize(); } /*自定义的归约函数，sendbuf为发送缓冲区，recvbuf为接收缓冲区，count为数据个数，root为指定根节点*/ /*该函数实现归约求和的功能*/ void Myreduce(double *sendbuf,double *recvbuf,int count,int root) { MPI_Status status; int i; int myid,numprocs; *recvbuf=0.0; MPI_Comm_rank(MPI_COMM_WORLD, myid); MPI_Comm_size(MPI_COMM_WORLD,numprocs); double *tmp; //非root节点向root节点发送数据 if(myid!=root) { MPI_Send(sendbuf,count,MPI_DOUBLE,root,99,MPI_COMM_WORLD); } //root节点接收数据并对数据求和，完成规约操作 if(myid==root) { *recvbuf=*sendbuf; for(i=0;inumprocs;i++) { if(i!=root) { MPI_Recv(tmp,count,MPI_DOUBLE,i,99,MPI_COMM_WORLD,status); *recvbuf=*recvbuf+*tmp; } } } } 《云计算与大数据技术》 并行计算实验室