MPI—统计素数个数

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();
}

小数据运行结果：

大数据运行结果：

实验加速比：

加速比计算：0.209725/0.105496=1.988