Particle-Mesh Ewald(PME)算法的GPU加速

计算物理 ›› 2010, Vol. 27 ›› Issue (4): 548-554.

Particle-Mesh Ewald(PME)算法的GPU加速

徐骥^1,2, 葛蔚¹, 任瑛¹, 李静海¹

1. 中国科学院过程工程研究所多相反应实验室, 北京 100190;
2. 中国科学院研究生院, 北京 100039

收稿日期:2009-04-30 修回日期:2009-09-04 出版日期:2010-07-25 发布日期:2010-07-25
作者简介:徐骥(1984-),男,江苏扬州,博士生,从事复杂分子体系的微观模拟方面的研究.
基金资助:
国家自然科学基金(20221603);国家科技支撑计划(2008BAF33B01);中国科学院知识创新工程(KGCX2-YW-124)资助项目

Implementation of Particle-Mesh Ewald(PME) on Graphics Processing Units

XU Ji^1,2, GE Wei¹, REN Ying¹, LI Jinghai¹

1. State Key Laboratory of Multiphase Complex System, Institute of Process Engineering, Chinese Academy of Sciences, Beijing 100190, China;
2. Graduate University of the Chinese Academy of Sciences, Beijing 100039, China

Received:2009-04-30 Revised:2009-09-04 Online:2010-07-25 Published:2010-07-25

摘要/Abstract

摘要： 讨论在NVIDIACUDA开发环境下,用GPU加速分子动力学模拟中静电作用的长程受力计算部分.采用Particle-Mesh Ewald(PME)方法,将其分解为参数确定、点电荷网格离散、离散网格的傅立叶变换、静电热能求解与静电力求解5个部分,并分别分析各部分的GPU实现.此方法已成功用于7个不同大小的生物分子体系的模拟计算,达到了7倍左右的加速.该程序可耦合到现有分子动力学模拟软件中,或作为进一步开发的GPU分子动力学程序的一部分,显著加速传统分子动力学程序.

关键词: PME(Particle-Mesh Ewald)加速, GPU(Graphic Processing Unit图形处理单元), CUDA(Compute Unified Device Architecture计算统一设备架构)

Abstract: In this article,long range part of electrostatic interaction is accelerated using GPUs in NIVIDIA CUDA programming environment.Particle Mesh Ewald algorithm is adopted and split into 5 procedures:parameters designation,discretization of point charges into grids,Fourier transformation of grids,potential of electrostatic interaction,and force calculation of electrostatic interaction.The codes are tested by seven biomolecular systems with different sizes.Roughly 7-fold speedup over one core of mainstream CPUs is obtained.The codes can be integrated into molecular dynamic simulation software packages that already exst or used as part of GPU codes developed in the future to further speedup traditional MD simulations.

Key words: PME, GPU, CUDA, speedup

中图分类号:

O241

徐骥, 葛蔚, 任瑛, 李静海. Particle-Mesh Ewald(PME)算法的GPU加速[J]. 计算物理, 2010, 27(4): 548-554.

XU Ji, GE Wei, REN Ying, LI Jinghai. Implementation of Particle-Mesh Ewald(PME) on Graphics Processing Units[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2010, 27(4): 548-554.

[1]	魏雪丹, 戴厚平, 李梦军, 郑洲顺. 一维空间Riesz分数阶对流扩散方程的格子Boltzmann方法[J]. 计算物理, 2021, 38(6): 683-692.
[2]	房尧立, 王一. 一种求解二维辐射流体力学方程组的显隐式格式[J]. 计算物理, 2021, 38(4): 401-417.
[3]	徐潜龙, 李晔. 有限深水域中三维频域自由面格林函数及其导数的计算方法[J]. 计算物理, 2021, 38(1): 35-46.
[4]	赵菲, 盛志强, 袁光伟. 基于二阶格式的有限体积保正格式[J]. 计算物理, 2020, 37(4): 379-392.
[5]	徐骁, 高志明, 戴自换. 三维拉氏理想磁流体数值模拟方法[J]. 计算物理, 2020, 37(4): 403-412.
[6]	张林, 葛永斌. 二维半线性扩散反应方程的高精度全隐格式及其多重网格方法[J]. 计算物理, 2020, 37(3): 307-319.
[7]	刘妍, 茅德康. 相容拉氏方法中的ADRS近似Riemann解算器[J]. 计算物理, 2020, 37(2): 140-152.
[8]	丁琪, 尚月强. 非定常Navier-Stokes方程基于两重网格离散的有限元并行算法[J]. 计算物理, 2020, 37(1): 10-18.
[9]	刘全, 倪国喜, 牛霄, 胡军. Vlasov-Poisson方程半拉格朗日守恒算法[J]. 计算物理, 2020, 37(1): 19-25.
[10]	赵国忠, 蔚喜军, 郭虹平, 董自明. 求解含有高阶导数偏微分方程的局部间断Petrov-Galerkin方法[J]. 计算物理, 2019, 36(5): 517-532.
[11]	孔令海, 孔令波, 许海波, 贾清刚. 含Laplace-Gauss型混合噪声图像二阶正则化重建方法[J]. 计算物理, 2019, 36(3): 280-290.
[12]	郭子滔, 冯仁忠. 一种高精度的修正Hermite-ENO格式[J]. 计算物理, 2019, 36(2): 141-152.
[13]	戴自换. ALE计算中基于声速分布的网格优化技术[J]. 计算物理, 2019, 36(1): 15-24.
[14]	徐丞君, 徐胜利. SPH二阶粒子近似光滑函数的充分条件及数值验证[J]. 计算物理, 2019, 36(1): 25-38.
[15]	符芳芳, 孔令华, 王兰, 徐远, 曾展宽. 一维Gross-Pitaevskii方程的高阶紧致分裂步多辛格式[J]. 计算物理, 2018, 35(6): 657-667.