MPI计算π

Posted zhangyazhou

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了MPI计算π相关的知识,希望对你有一定的参考价值。

MPI计算(pi)
  • 利用公式

[int_0^1 frac{4}{1+x^2}dx = pi ]

技术图片
技术图片
#include<stdio.h>
#include<mpi.h>
#include<stdlib.h>
#include<time.h>
int main(int argc, char** argv)
{
	int rank, size;
	int n = 100;
	double sum,width,local,mypi,pi;
	MPI_Init(&argc, &argv);
	MPI_Comm_rank(MPI_COMM_WORLD, &rank);
	MPI_Comm_size(MPI_COMM_WORLD, &size);
	MPI_Bcast(&n, 1, MPI_INT, 0, MPI_COMM_WORLD);
	sum = 0.0;
	width = 1.0 / n;
	/*
	* 把0~1分成n份,把计算任务分成size份。
	*/
	for (int i = rank; i < n; i += size)
	{
		local = width * ((double)i + 0.5);
		sum += 4.0 / (1.0 + local * local);
	}
	mypi = width * sum;
	MPI_Reduce(&mypi, &pi, 1, MPI_DOUBLE, MPI_SUM, 0, MPI_COMM_WORLD);
	if (rank == 0)
	{
		printf("pi is %.20f
", pi);
		fflush(stdout);
	}
	MPI_Finalize();
	return 0;
}

以上是关于MPI计算π的主要内容,如果未能解决你的问题,请参考以下文章

MPI和OpenMP混合编程计算pi π值

MPI从双数组计算差异符号

在python的IDLE中计算Π值的代码错误在哪?

使用 PBLAS 算法进行矩阵计算 - 使用 MPI 的多线程

java运算中,如何引入π进行计算(希望能附上例子)?

如何在Python上实现用文本进度条体现π的计算过程