【高性能并行计算-MPI】求取圆周率PI的值
本帖最后由 DAY 于 2018-12-4 09:00 编辑原理说明:
我们通过求解定积分|1,0（4/1+x^2）（定积分不清楚的可以看看后面手工画的图片）来确定PI的值，首先我们把区间【0，1】分解为K个长子区间，用这些高度可以构造k个矩行，这些矩形的面积逼近曲线下图形的面积，当k增加时，估算精度随之增加。
代码分析:
#include<stdio.h>
#include<mpi.h>
const int INTERVALS = 100000000;
int main(int argc,char*argv[])
{
int rank,numProcs;//定义进程号和进程数目
double area;//记录每个进程计算的面积
double PI;//最后的PI
MPI_Init(&argc,&argv);
MPI_Comm_rank(MPI_COMM_WORLD,&rank);
MPI_Comm_size(MPI_COMM_WORLD,&numProcs);
//进行进程分配
for(int i = rank; i < INTERVALS;i += numProcs ){
area += 4.0 /(1.0+(1.0/INTERVALS * i)*(1.0/INTERVALS *i))*(1.0/INTERVALS);
}
//规约面积和
MPI_Reduce(&area,&PI,1,MPI_DOUBLE,MPI_SUM,0,MPI_COMM_WORLD);
//阻塞
MPI_Barrier(MPI_COMM_WORLD);
//打印输出
if(0 == rank){
printf("PI是%.12lf\n",PI);
}
MPI_Finalize();
return 0;
}
编译和运行:
页:
[1]