并行计算MPI(三):计算π

最新推荐文章于 2024-08-06 08:11:08 发布

sj_lxd

最新推荐文章于 2024-08-06 08:11:08 发布

阅读量1.7w

点赞数 10

分类专栏：并行计算文章标签：并行计算数学

本文链接：https://blog.csdn.net/sj_lxd/article/details/44197863

版权

学习了之前的一些并行计算MPI基本知识之后，其实可以尝试应用起来，今天这篇文章就是用并行计算求解π的。首先需要知道π的求解方式：

上面是高等数学里的基本知识，所以π的求解就可以变为求解y(x)=4/(1+x^2)在(0,1)区间的面积。而且求解面积这一块可以采用将(0,1)平均分n个块，将n个矩形面积相加即可得到近似解。而且n取得越大，近似解越精确。

比如我们取n=100,开4个进程来处理，那么相应的每个进程平均分到计算25个矩形面积，如果开5个进程，那么平均分下来每个进程计算20个矩形面积。下面列出代码：

//返回函数值，因为调用次数很多，这里写成inline提高效率
inline double f(double x) {
	return 4/(1+x*x);
}

int main(int argc, char *argv[]) {
	double pi, h, sum, x, startime, endtime;
	int size, myid;
	long long n;
	MPI_Init(&argc, &argv);
	MPI_Comm_rank(MPI_COMM_WORLD, &myid);
	MPI_Comm_size(MPI_COMM_WORLD, &size);

	n=0;
	if(0 == myid) {
		cout <