多线程并行快速求解e值的六种方法.pdfVIP

文章编号：1007—1423(2013)17—0015—06 DOI：10．3969／j．issn．1007—1423．2013．17．004 多线程并行快速求解 e值的六种方法 朱建伟 ．刘 荣 (中国石油大学 (华东)理学院，青 岛 266580) 摘 要 ：欧拉常数 e是数学最重要的符号之一，通常利用叠加无穷级数有限项来进行计算。采用多 线程技术，将 串行算法改为并行算法 ，通过使用VC-inAPI、OpenMP、MPI三种并行模式得到 六种计算方法对 e值进行计算。根据数值试验结果，对各种并行方法优缺点进行 比较 ，发现 MPI并行模式计算速度最快，最终获得结论并行计算模式可推广到各种数值 问题 。 关键词 ：并行计算 ；多线程 ；加速 比；XX-／inAPI；OpenMP；MPI 0 引 言 进行改造 ，得出优化后算法。进而考虑将并行计算中高 效的计算方式移植到其他数值 问题 。 自然常数 e最先 由欧拉在 1727年使用 ．具体值为 e：2．718281828…。e是个超越数．是 自然对数的底 ．并 1 串行算法描述 且在Euler公式中e也十分重要 ．国际上对 e值计算方 国际上计算 e值的算法有很多．比较常用有极限 法研究从未间断。文献 3『1与4『1介绍了快速精确计算e 逼近法、无穷级数逼近法、连分数法、Pippenger积法、 值的多种方式 ．计算精度高。但是所有算法均为串行算 Monte—Carlo模拟法等。由于本文 目的在于将某串行算 法 ．没有发挥出多核 CPU在并行计算发面的强大优 法并行实现 ．级数逼近法计算机实现起来比较方便，所 势 ．在高性能计算领域对 e值并行计算却未能有较为 以选取级数逼近法进行改造 ，下面进行简要描述 。 详细的实现 为弥补并行计算领域对 e值求解不够系 利用泰勒展开法将指数函数e．令x=l展开成麦 统的缺憾 ．也为充分发挥多CPU在计算时的强大功 克劳林级数 能 ．依托于级数求解算法 。利用三种并行模式改造成六 种并行算法并进行程序实现。 e= + }+++} 。 WinAPI方式就是编写线程 函数 ．利用 Windows系 作近似计算时．采用前有限项的和 作为 e的近似 统中的API接 口进行线程调度实现并行 ：OpenMP方式 值 是使用 OpenMP语言．利用封装好 的函数进行并行计 算 ．不用考虑底层实现 ：而 MPI方式则是利用 MPI并 + }+++} 行语言．以消息通信的模式 ．利用多节点或多CPU模 进行近似那么误差为 可以解得为En=e—e，显然 拟多节点进行并行计算。然后对于各种实现方式。获得 级数计算 e的收敛速度是 0( )阶的。 其计算速度、加速 比及加速效果。依托于计算机数值计 串行实现作为最基本的实现方式，其实现的核心 算大数计算 比大数乘法计算更加快速的特性 。对算法 ★基金项 目：国家大学生创新训练计划 (No．111042566) 收稿 日期 ：2013—05—07 修稿 日期 ：2013—05—27 作者简介 ：朱建伟(1991一)，男，山东 日照人 ，本科 ，研究方向为高性能并