[論文レビュー] Accelerating glassy dynamics using graphics processing units

Peter H. Colberg, Felix Höfling|arXiv (Cornell University)|Dec 20, 2009

Theoretical and Computational Physics参考文献 1被引用数 2

ひとこと要約

この論文は、CUDAを用いてGPUアクセラレーションされた分子動力学（MD）実装を提示し、最大100万粒子の系をシミュレートし、シリアルCPUコードと比較して最大80倍の高速化を達成した。単一精度浮動小数点演算による数値不安定性を解消するため、二重単一精度のエミュレーションを採用し、10⁸ MDステップにわたり厳密なエネルギー保存を確保し、過冷状態のLennard-Jones混合系において物理的に誤った結果を防いだ。

ABSTRACT

Modern graphics processing units (GPUs) provide impressive computing resources, which can be accessed conveniently through the CUDA programming interface. We describe how GPUs can be used to considerably speed up molecular dynamics (MD) simulations for system sizes ranging up to about 1 million particles. Particular emphasis is put on the numerical long-time stability in terms of energy and momentum conservation, and caveats on limited floating-point precision are issued. Strict energy conservation over 10^8 MD steps is obtained by double-single emulation of the floating-point arithmetic in accuracy-critical parts of the algorithm. For the slow dynamics of a supercooled binary Lennard-Jones mixture, we demonstrate that the use of single-floating point precision may result in quantitatively and even physically wrong results. For simulations of a Lennard-Jones fluid, the described implementation shows speedup factors of up to 80 compared to a serial implementation for the CPU, and a single GPU was found to compare with a parallelised MD simulation using 64 distributed cores.

研究の動機と目的

• 最新のGPUを用いてガラス状態の系の分子動力学シミュレーションを高速化すること。
• GPUの浮動小数点精度が限られているにもかかわらず、エネルギーおよび運動量保存の長期的数値安定性を確保すること。
• 過冷却液体の遅い動的挙動において、単一精度算術が物理的に誤った結果を生じることを示すこと。
• 1つのGPUで64コアの分散CPUクラスタと同等の性能を達成すること。

提案手法

• CUDAプログラミングインターフェースを活用し、MDシミュレーションをGPUアーキテクチャに移植すること。
• 高精度浮動小数点演算をエミュレートするため、精度が重要なセクションに二重単一精度算術を実装すること。
• 長時間の安定性および物理的正確性を評価するため、テスト用に二成分Lennard-Jones混合系を用いること。
• シリアルCPUコードおよび64コアに並列化されたCPU実装と比較して、パフォーマンスをベンチマークすること。
• GPUのメモリ階層と並列性を最適化した力の計算を伴う標準的なMD積分スキームを適用すること。

実験結果

リサーチクエスチョン

• RQ1 GPUアクセラレーションは、ガラス形成系の分子動力学シミュレーションにおいて顕著な高速化を達成できるか？
• RQ2 単一精度浮動小数点算術は、過冷却液体の長時間MDシミュレーションにおいて物理的に誤った挙動を引き起こすか？
• RQ3 二重単一精度エミュレーションは、GPU上で10⁸ MDステップにわたりエネルギー保存を厳密に確保できるか？
• RQ4 MDシミュレーションにおいて、GPUのパフォーマンスは高度に並列化されたCPUクラスタと比べてどうか？

主な発見

• GPU実装は、シリアルCPU実装と比較して最大80倍の高速化を達成した。
• 1つのGPUが64コアの分散CPU並列化と同等の性能を発揮した。
• 二重単一精度エミュレーションにより、10⁸ MDステップにわたり厳密なエネルギー保存が可能になり、数値的安定性が確保された。
• 単一精度算術は、過冷却二成分Lennard-Jones混合系のシミュレーションにおいて、定量的および物理的に誤った結果を生じさせた。
• 本手法は運動量保存を維持し、長時間スケールのシミュレーションにおいても安定であった。