【高性能并行计算-MPI】求取圆周率PI的值

DAY

原理说明:
我们通过求解定积分|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;
}

编译和运行: