Scalable and portable implementation of the fast multipole method on parallel computers

详细信息    查看全文

• 作者：Ogata ; Shuji ; Campbell ; Timothy J. ; Kalia ; Rajiv K. ; Nakano ; Aiichiro ; Vashishta ; Priya ; Vemparala ; Satyavani
• 关键词：02.60.-x ; 02.70.-c ; Fast multipole method ; Parallel computation ; Stress calculation ; Periodic boundary conditions ; Coulomb interaction
• 刊名：Computer Physics Communications
• 出版年：2003
• 出版时间：July 1, 2003
• 年：2003
• 卷：153
• 期：3
• 页码：445-461
• 全文大小：366 K

文摘

A scalable and portable Fortran code is developed to calculate Coulomb interaction potentials of charged particles on parallel computers, based on the fast multipole method. The code has a unique feature to calculate microscopic stress tensors due to the Coulomb interactions, which is useful in constant-pressure simulations and local stress analyses. The code is applicable to various boundary conditions, including periodic boundary conditions in two and three dimensions, corresponding to slab and bulk systems, respectively. Numerical accuracy of the code is tested through comparison of its results with those obtained by the Ewald summation method and by direct calculations. Scalability tests show the parallel efficiency of 0.98 for 512 million charged particles on 512 IBM SP3 processors. The timing results on IBM SP3 are also compared with those on IBM SP4.