MPI—统计素数个数

爱被打了一巴掌 2022-07-14 13:55 176阅读 0赞

## MPI求素数个数 ##

// MPI1.cpp : 定义控制台应用程序的入口点。
    //MPI求素数个数
    
    #include "stdafx.h"
    #include "mpi.h"
    #include <stdio.h>
    #include <math.h>
    
    int isPrime(int lyy_num)   //判断是否为素数  
    {  
        int lyy_flag=1;  
        int lyy_s=sqrt(lyy_num*1.0);  
        for(int j=2;j<=lyy_s;j++)  
            {  
                if(lyy_num%j==0)     
                {  
                    lyy_flag=0;  
                    break;  
                }  
            }  
        return lyy_flag;  
    }  
    
    void main(int argc, char * argv[])
    {
    	int n=0,myid,numprocs,i,pi,sum,mypi;
    	double startwtime,endwtime;
    
    	MPI_Init(&argc,&argv);
    	MPI_Comm_size(MPI_COMM_WORLD,&numprocs);
    	MPI_Comm_rank(MPI_COMM_WORLD,&myid);
    
    	if(myid==0)
    	{
    		printf("输入一个数字：");
    		fflush(stdout);
    		scanf_s("%d",&n);
    		startwtime=MPI_Wtime();
    	}
    	MPI_Bcast(&n,1,MPI_INT,0,MPI_COMM_WORLD);   //将n值广播出去
    	sum=0;
    	for(i=myid*2+1;i<=n;i+=numprocs*2)
    		sum+=isPrime(i);
    	mypi=sum;  
    	MPI_Reduce(&mypi,&pi,1,MPI_INT,MPI_SUM,0,MPI_COMM_WORLD);  //规约
    	if(myid==0)      
    	{
    		printf("结果=%d\n",pi);
    		endwtime=MPI_Wtime();
    		printf("并行时间=%f\n",endwtime-startwtime);
    	}
    	//串行程序
    	sum=0;
    	double startwtime2=MPI_Wtime();
    	if(myid==0)
    	{
    		for(i=1;i<=n;i+=2)
    			sum+=isPrime(i);
    		double endwtime2=MPI_Wtime();
    		printf("结果=%d\n",sum);
    		printf("串行时间=%f\n",endwtime2-startwtime2);
    	}
    
    	MPI_Finalize();
    }

## 小数据运行结果： ##

![Center][]

## 大数据运行结果： ##

## ![Center 1][]   实验加速比：       ##

加速比计算：0.209725/0.105496=1.988

[Center]: /images/20220714/6f131733feaa4d43916e9ddc55048edc.png
[Center 1]: /images/20220714/ba591fe6e4154ab1940cf3938bb0c41c.png