复化梯形求定积分的并行计算方法.pdfVIP

理论探讨/环球市场 复化梯形求定积分的并行计算方法 杨菁蓓 中国石油大学(华东)理学院 摘要：使用复化梯形法求积分，需要将区间 [0,1] 进行分割，分 int threadID=*((int*)arg); 割数越大求得的积分近似值就越精确。但是当分割数 N 较大时，算 double a=0.0,b=1.0; 法的计算量会比较大，算法的运行效率比较低。可以通过并行计算 double s1=(f1(a)+f1(b))/2; 的方法提高运算的速度。这样在达到较高计算精度的同时提高了算 double h=(b-a)/numSteps; 法运行的速度。 for (int i=threadID;inumSteps;i+=numThreads) 关键词：复化梯形法；并行计算；定积分 { s1+=f1(a+i*h); 一、某问题的串行算法描述 } 1.1 实现梯形积分法的基本步骤 s1*=h; (1) 输入区间 [a,b] 的端点 a,b 值以及分割数 N； EnterCriticalSection(cs); (2) 将 区 间 [a,b] 等 分 成 N 个 小 区 间， 每 一 个 小 区 间 的 长 度 s+=s1; h=(b-a)/N; LeaveCriticalSection(cs); (3) 计算每一个等分点的函数值 y f a ih i N return 0; i= ( + )( =0,1,…, )； (4) 计算 } 3.3 基于OpenMP 的多核并行算法的实现 二、该问题的并行算法描述 OpenMp 提供了对并行算法的高层的抽象描述，通过在源代码中 2.1 基于API 的多核并行算法的设计 加入专用 pragma 来指明自己的意图，由此编译器可以自动将程序进 Windows 系统提供很多 API 接口，可以利用 WINAPI 定义线程 行并行化，并在必要之处加入同步互斥以及通信。归约操作的方式 函数，函数内部设计好改线程所要进行的工作，然后主函数里创建 是拷贝复制执行 线程，将线程函数导入创建好的线程中运行，计算根据创建的线程 3.4 基于MPI 的并行算法实现 数目，调用相应数目的 CPU 进行计算，保证每个 CPU 运行一个线 MPI 是基于消息传递的并行计算模式，建立消息传递标准的主 程，最终计算结束后将结果合并可得到正确结果。 要优点是可移植性和易于使用。以低级消息