MPI求素数个数

最新推荐文章于 2020-08-15 23:06:24 发布

fxy流年无悔

最新推荐文章于 2020-08-15 23:06:24 发布

阅读量3.1k

点赞数

文章标签：并行计算

本文链接：https://blog.csdn.net/a716121/article/details/53819234

版权

博客探讨了在MPI并行计算中求素数个数时遇到的问题及解决方案。错误分析指出MPI_Bcast的用法不当，应分别广播输入数值，避免广播数组多个元素。还提到负载不平衡问题，因偶数素数较少导致进程间工作不均衡，影响并行效率。建议改进方法是更好地分配任务以实现真正的并行计算。

摘要由CSDN通过智能技术生成

附上另一个调用函数后针对一个输入正确代码：

#include "stdafx.h"
#include "mpi.h"
#include "stdio.h"
#include "string.h"
#include "math.h"
int isPrime(int n)
{
        int  flag=1;
		for(int m=2;m<=sqrt(n*1.0);m++)
		{
			if(n%m==0)
			{
				flag=0;
				break;
			}
		}
		return flag;
}
int main(int argc, char *argv[])
{

    int a,mypid,sum=0,local=0,size;
	int flag;
    double starttime1,endtime1,starttime2,endtime2,t1,t2;
    MPI_Init(&argc,&argv);
    MPI_Comm_rank(MPI_COMM_WORLD,&mypid);
    MPI_Comm_size(MPI_COMM_WORLD,&size);
    if(mypid==0)
    {
        scanf_s("%d",&a);
    }
    MPI_Bcast(&a,1,MPI_INT,0,MPI_COMM_WORLD);
	starttime1=MPI_Wtime();
    for(int i=mypid*2+1; i<=a; i+=size*2) //每个进程算自己的任务，保存在local中
    {
		local+=isPrime(i);
 
    }
    MPI_Reduce(&local,&sum,1