并行计算MPI(三):计算π |
您所在的位置：网站首页 › 级数公式求圆周率c语言程序 › 并行计算MPI(三):计算π |
并行计算MPI(三):计算π
|
学习了之前的一些并行计算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 |
【本文地址】
今日新闻 |
推荐新闻 |