[論文レビュー] GALAXY package for N-body simulation

J. A. Sellwood|arXiv (Cornell University)|Jun 25, 2014

Astronomy and Astrophysical Research被引用数 20

ひとこと要約

GALAXYパッケージは、粒子・メッシュ法を用い、複数の力計算手法を備えた、孤立した衝突なし星系（例：銀河）をモデル化する高度に最適化されたN体シミュレーションコードである。特定の応用、特に大規模な球状およびディスク型系において、広く使われているGadget-2ツリー法コードと比較して50～200倍の高速性能を達成しており、グリッドベースの重力ソルバーとリアルタイム分析ツールにより高い精度を維持している。

ABSTRACT

This posting announces public availability of the GALAXY software package developed by the author over the past 40 years. It is a highly efficient code for the evolution of (almost) isolated, collisionless stellar systems, both disk-like and ellipsoidal. In addition to the N-body code galaxy, which offers eleven different methods to compute the gravitational accelerations, the package also includes sophisticated set-up and analysis software. This paper gives an outline of the contents of the package and provides links to the source code and a comprehensive on-line manual. While not as versatile as tree codes, the particle-mesh methods in this package are shown, for certain restricted applications, to be between 50 and 200 times faster than a widely-used tree code.

研究の動機と目的

• 孤立星系、特に銀河を対象とした高性能で公開可能なN体シミュレーションパッケージの提供。
• 衝突なしシミュレーションにおける計算ボトルネックを、最適化された粒子・メッシュ法による効率最大化によって克服すること。
• ディスクにバルジやハローを含む複雑な多成分系に対して、正確な初期条件生成とリアルタイム分析を支援すること。
• コードの高い学習曲線にもかかわらず、研究者がより高解像度のシミュレーションを、数値的整合性を保ちながら実行できるようにすること。
• 長期間にわたる保守性と拡張性を備えた、GNU GPLの下での無料かつオープンソースの代替手段を提供すること。

提案手法

• コアN体コードであるgalaxyは、球状および円柱状グリッド構成と調和分解を含む、11種類の異なる粒子・メッシュ法を用いて重力加速度を計算する。
• FFTPAKからのFFTベースのソルバーを採用し、2次元および3次元グリッド幾何学をサポート。力の精度を高めるために、変動する調和モード（例：l ≤ 4, m ≤ 8）を用いる。
• リアルタイムでの分析を実装し、各タイムステップで主要な動力学的量（例：エネルギー、角運動量、密度プロファイル）を記録することで、ストレージを最小限に抑え、コンactなデータ出力を可能にする。
• ディスク、バルジ、ハローを含む複数の質量成分を持つ平衡初期条件を生成するための特別な初期化ツール（例：begin, start）をパッケージに含む。
• 単一プロセッサおよびMPI並列実行を別々の実行可能ファイル（例：galaxy_mpi）でサポートし、マルチコアシステムにおけるスケーラブルなパフォーマンスを実現する。
• 外部ライブラリと統合：FFTsにはFFTPAK、プロットにはPGPLOT、任意の高度機能にはNAGをサポート。オープンソース代替品へのフォールバックも実装済み。

実験結果

リサーチクエスチョン

• RQ1 粒子・メッシュ法をどのように最適化すれば、孤立銀河の衝突なしN体シミュレーションにおいて極めて高い性能を達成できるか？
• RQ2 特定の銀河モデルにおいて、Gadget-2のような既存のツリー法コードと比較して、GALAXYパッケージは計算効率でどの程度優れているか？
• RQ3 N体シミュレーションにおけるリアルタイム分析は、科学的有用性を損なわせることなく、データストレージ要件を顕著に削減できるか？
• RQ4 高パフォーマンスであるが複雑なコードGALAXYを、天文学的研究における新規ユーザーが使用する際の実用的制限およびリスクは何か？
• RQ5 粒子数およびプロセッサ数の増加に伴い、GALAXYのパフォーマンス優位性は、異なる系の幾何形状においてどのようにスケーリングするか？

主な発見

• 球状プラマー・モデルでは、テストされた全粒子数およびプロセッサ数において、Gadget-2と比較してGALAXYコードが約50倍高速に動作する。
• 厚みをもたせたクズミン＝クツェフォ・ディスク・モデルでは、粒子数の増加に伴いGALAXYの性能優位性が増大し、N ≈ 10^8粒子でGadget-2と比較して200倍以上高速に達する。
• GALAXYにおける粒子・メッシュ法は、固定グリッドのオーバーヘッドが一定であるため、Nの増加に伴うパフォーマンス劣化が最小限に抑えられ、大規模シミュレーションに効率的である。
• リアルタイム分析により、後処理を要せず、各分析ステップで結果をコンactに保存できるため、データストレージの必要量が顕著に削減される。
• 単一プロセッサでも高いパフォーマンスを発揮し、MPI対応実行可能ファイルを用いることで、並列処理により巨大なスピードアップが可能である。
• 高速である一方で、ユーザーの熟練度が求められる。研究者が出版可能な結果を得る前に、数週間を費やしてコードを習得することが推奨される。