并行计算划分和分治策略.ppt

* * * * * * * * * * * * * * * * * * * * * * * * 采用第二种方法，即计算每个小区间的梯形面积的计算公式： d (f(a)+f(a+d)) 2 d (f(a+d)+f(a+2d)) 2 d (f(a+(n-1)d)+f(b)) 2 面积= + +…+ = d + f(a+d) + f(a+2d) +…+ f(a+(n-1)d) + f(b) 2 f(a) 2 算法只需做少量修改： part_sum = 0.5 *(f(start)+f(end)); for (x = start + d; xend; x = x + d) { //calcs partial sums part_sum += f(x); } part_sum = d * part_sum; 除以上计算数值积分方法，我们还可以用递归的分治法来解决问题。如用C-Linda 实现的计算 p 的程序。 C-Linda 含有几个在“元组空间”上的操作。 元组空间(tuple space)：并行编程环境中的虚拟存储器，由一组有序的元组构成。并行执行的进程通过共享数据元组进行交互。 见 ComputePi . doc 文件。 元组空间 进程 进程 进程 元组 元组 元组 work(1, 0, 400000 , 1.0/400000 ) work(1, 0, 400000 , 0.0000025 ) work(2, 0, 200000 , 0.0000025 ) work(3, 200000, 400000 , 0.0000025 ) work(1, 0, 400000 , 0.0000025 ) work(2, 0, 200000 , 0.0000025 ) work(3, 200000, 400000 , 0.0000025 ) work(4, 0, 100000 , 0.0000025 ) work(5, 100000, 200000 , 0.0000025 ) work(6, 200000, 300000 , 0.0000025 ) work(7, 300000, 400000 , 0.0000025 ) ”worker”, 1, work(2, 0, 200000, 0.0000025 ) ”worker”, 1, work(3, 200000, 400000, 0.0000025 ) eval ”worker”, 2, work(4, 0, 100000, 0.0000025 ) ”worker”, 2, work(5, 100000, 200000, 0.0000025 ) ”worker”, 3, work(6, 200000, 300000, 0.0000025 ) ”worker”, 3, work(7, 300000, 400000, 0.0000025 ) ”worker”, 2, 0.9799 ”worker”, 2, 0.8747 ”worker”, 3, 0.7194 ”worker”, 3, 0.5676 work(1, 0, 400000 , 0.0000025 ) work(2, 0, 200000 , 0.0000025 ) work(3, 200000, 400000 , 0.0000025 ) in ”worker”, 1, 1.8546 ”worker”, 1, 1.2870 work(1, 0, 400000 , 1.0/400000 ) in eval : : 问题的描述： N体问题是利用牛顿的万有引力定律和运动定律模拟太空中星体的运动轨迹。 万有引力定律：质量为 ma 和 mb 的两个物体间的引力为： 应用3 ---- N体问题 其中：G (6.67259*10-11米3/(千克.秒2)) 是引力常数，r 为两物体间的距离。 G ma mb r 2 F = ————— 应用3 ---- N体问题 一个物体在受力的情况下，将根据牛顿第二定律产生加速度： F = ma m 是物体的质量，F 是物体所受的力，a 是物体获得的加速度 。 实现细节： 设模拟天体运动的时间间隔为Δt，物体的质量为m，物体所受的力为： m(vt+1 - vt) Δt F =