快速球谐函数展开的并行算法设计及实现  

Design and Implementation of Parallel Algorithm for Fast Spherical Harmonic Expansions

作　　者：王翔[1] 宋君强[1] 卢风顺[1] 杨锦辉[1] 

机构地区：[1]国防科学技术大学计算机学院,湖南长沙410073

出　　处：《微电子学与计算机》2011年第8期26-31,共6页Microelectronics & Computer

基　　金：国家自然科学基金项目(40505023)

摘　　要：球谐函数展开是一种广泛应用的数学方法.在研究Rokhlin-Tygert球谐函数展开快速算法(RT算法)的基础上,扩充并完善了RT算法,建立了所有阶连带Legendre多项式展开系数的计算过程,进而研究了该算法的MPI并行策略并给出了RT并行算法.数值实验表明,RT算法的计算量随三角截断波数的增大而急剧增长;MPI并行技术能够有效提高其运算速度,较好地提升算法的整体性能,但并行效率会受多核处理器cache大小和访存带宽的限制.Spherical harmonic expansions are a widely used tool of applied mathematics.Based on the research of Rokhlin-Tygert＇s fast algorithm for spherical harmonic expansions（RT algorithm）,the procedure for computing expansion coefficients of all orders of associated Legendre functions is established,while the RT algorithm is extended and consummated.Furthermore,the strategy of MPI parallelization of RT algorithm is developed,and its parallel algorithm is designed and implemented.Numerical results show that,the computation of RT algorithm grows violently with the lift of the highest Fourier wave number of triangular truncations;MPI parallelization can efficiently bring up the computing speed and improve the performance of the algorithm.However,the efficiency of parallelization could be limited by the cache size of multi-cores processors and the bandwidth of memory access.

关 键 词：球谐函数 连带Legendre多项式 MPI 并行算法 

分 类 号：TP312[自动化与计算机技术—计算机软件与理论]