MPI 实现梯形积分法

最新推荐文章于 2021-10-14 15:06:16 发布

白菜白菜白

最新推荐文章于 2021-10-14 15:06:16 发布

阅读量6.3k

点赞数 1

分类专栏： parallel computing

本文链接：https://blog.csdn.net/u012417189/article/details/25717263

版权

本文介绍如何使用MPI实现梯形积分法的并行计算，通过将任务分解、识别通信信道、任务聚合及核上的任务分配，有效提高计算效率。内容包括每个小梯形面积的计算及所有梯形面积的总和计算，确保n是MPI进程数（comm_sz）的整数倍。

摘要由CSDN通过智能技术生成

梯形面积和 = h * ( f(x0)/2 + f(x1) + ....... + f(xn)/2) h = (b-a)/n;

// 正常版伪代码
h= (b-a)/n;
approx = (f(a)+f(b))/2.0;
for (int i=0;i<n;i++)
{
 	x_i = a + i*h;
	approx += f(x_i);
}
approx = h*approx;

并行化需要考虑这几个步骤 1.将问题的解决方案划分为多个任务 2.在任务间识别出需要的通信信道 3.将任务聚合成符合任务 4.在核上分配复合任务

计算每一个小梯形的面积和计算这些梯形的面积和

Get a,b,n;
h = (b-a)/n;
local_n = n/comm_sz;
local_a = a + my_rank*local_n*h;
local_b = local_a +  local_n*h;
local_integral = Trap(local_a,local_b,local_n,h);  // 面积的串行实现
if(my_rank !=0){
send local_intergral to process 0;
}
else
 total_intergral = local_intergral;
for(proc = 1; proc<comm_sz;proc++){
receive local_intergal from proc;
total_