MPI—并行求和

素颜马尾好姑娘i 2022-07-15 10:49 256阅读 0赞

// MPI1.cpp : 定义控制台应用程序的入口点。
    //MPI求和
    
    #include "stdafx.h"
    #include "mpi.h"
    #include <stdio.h>
    #include <math.h>
    
    double f(double a)  //定义函数f(x)
    {
    	return (4.0/(1.0+a*a));
    }
    
    void main(int argc, char * argv[])
    {
    	int done=0,n,myid,numprocs,i;
    	double mypi,pi,sum;
    	double startwtime,endwtime;
    	int namelen;
    	char processor_name[MPI_MAX_PROCESSOR_NAME];
    
    	MPI_Init(&argc,&argv);
    	MPI_Comm_size(MPI_COMM_WORLD,&numprocs);
    	MPI_Comm_rank(MPI_COMM_WORLD,&myid);
    	MPI_Get_processor_name(processor_name,&namelen);
    
    	fprintf(stderr,"Process %d on %s\n",myid,processor_name);
    	fflush(stderr);
    
    	n=0;
    	while(!done)
    	{
    		if(myid==0)
    		{
    			printf("输入一个数字不超过：（0退出）");
    		    fflush(stdout);
    			scanf_s("%d",&n);
    			startwtime=MPI_Wtime();
    		}
    		MPI_Bcast(&n,1,MPI_INT,0,MPI_COMM_WORLD);   //将n值广播出去
    		if(n==0)
    			done=1;
    		else
    		{
    			sum=0.0;
    			for(i=myid+1;i<=n;i+=numprocs)
    				sum+=i;
    			mypi=sum;   //各个进程并行计算得到的部分和
    			MPI_Reduce(&mypi,&pi,1,MPI_DOUBLE,MPI_SUM,0,MPI_COMM_WORLD);
    			if(myid==0)       //执行累加的号进程，将近似值打印出来
    			{
    				printf("结果%.16f\n",pi);
    				endwtime=MPI_Wtime();
    				printf("时间=%f\n",endwtime-startwtime);
    			}
    		}
    	}
    	MPI_Finalize();
    }

一个处理器：

![Center][]

运行结果：

![Center 1][]

两个处理器（双核）：

![Center 2][]

运行结果：

![Center 3][]

[Center]: /images/20220715/2aec191ba4fb43b295c131a22102b7df.png
[Center 1]: /images/20220715/c423a915dfe5415c9d785613f9b4a2c2.png
[Center 2]: /images/20220715/dabe87c73a5e40ae93644782dbeab7f7.png
[Center 3]: /images/20220715/128fa6dc2eed489abf8c84b116dc3d58.png