[論文レビュー] SWIFT: Fast algorithms for multi-resolution SPH on multi-core architectures

Pedro Gonnet, Matthieu Schaller+2 more|arXiv (Cornell University)|Sep 15, 2013

Algorithms and Data Compression参考文献 4被引用数 8

ひとこと要約

本稿では、階層的セル分割とソート済み相互作用を用いることで、マルチコアシステム上のマルチスケールSPHを高速化する新しいタスクベースのSPHシミュレーションフレームワークSWIFTを提案する。SWIFTは32コアで75%を超える並列効率を達成し、宇宙論的シミュレーションにおいてGadget-2の8倍以上も高速である。これは、定数時間の近接粒子探索、対称性の活用、および明示的なSIMDベクタ化なしのキャッシュ効率の向上によるものである。

ABSTRACT

This paper describes a novel approach to neighbour-finding in Smoothed Particle Hydrodynamics (SPH) simulations with large dynamic range in smoothing length. This approach is based on hierarchical cell decompositions, sorted interactions, and a task-based formulation. It is shown to be faster than traditional tree-based codes, and to scale better than domain decomposition-based approaches on shared-memory parallel architectures such as multi-cores.

研究の動機と目的

• マルチコアアーキテクチャ上でスムージング長の動的レンジが大きいSPHシミュレーションにおける性能ボトルネックを解消すること。
• ツリーに基づく近接粒子探索手法の限界、すなわち低いキャッシュ効率、高い通信オーバーヘッド、対称性の活用不足を克服すること。
• ハードウェアに依存しない最適化に依存せずに、共有メモリ型SPHシミュレーションにおける並列スケーリングと計算効率を向上させること。
• 細粒度の負荷バランスと低遅延同期を実現する、タスクベースで非同期な並列化スキームを開発すること。
• 将来のハイブリッド共有/分散メモリ型およびGPUアクセcelerated SPHシミュレーションへの拡張を可能にすること。

提案手法

• 最大スムージング長と同等以上に大きなセルサイズを用いたシミュレーション領域の階層的セル分割を採用し、各粒子に対して定数時間で近接粒子を探索可能にする。
• 空間的近接度に基づいて事前にソートされた順序で粒子ペアを処理することで、重複する距離計算を最小限に抑えるための、ソート済み粒子相互作用を導入する。
• 各粒子ペアを1回だけ計算することで対称的相互作用を活用し、密度計算および力計算における二重作業を排除する。
• 非同期計算を管理するためのタスクベースの並列化モデルを採用し、細粒度の負荷バランスを実現するとともに、コア間の同期ポイントを削減する。
• ツリー走査を完全に回避し、体系的に列挙・処理されるセルペア相互作用に置き換える。
• 実装は標準的なコンパイラ最適化を用いるが、明示的なSIMDベクタ化は行わず、アルゴリズム的設計による性能向上を示している。

実験結果

リサーチクエスチョン

• RQ1 セルベースの近接粒子探索アプローチは、マルチスケールSPHにおいて、計算効率とキャッシュ利用効率の面でツリー基盤の手法を上回ることができるか？
• RQ2 従来のドメイン分割やツリー基盤の並列化と比較して、タスクベースの並列化は、共有メモリ型マルチコアシステムにおけるスケーリングと効率をどの程度向上させることができるか？
• RQ3 ソート済み相互作用と対称性の活用を組み合わせることで、SPHの近接粒子探索における重複計算はどの程度削減されるか？
• RQ4 明示的なSIMDベクタ化に依存せず、アルゴリズム的およびデータ構造的革新のみで顕著な性能向上を達成できるか？
• RQ5 提案手法は32コアでどの程度スケーリングするか？また、Gadget-2のような既存のコードと比較して、並列効率はどの程度か？

主な発見

• SWIFTは、全テストシミュレーションにおいて32コアで75%を超える並列効率を達成し、マルチコアアーキテクチャにおける強スケーリングを示した。
• 宇宙論的ボリュームシミュレーションにおいて、SWIFTは同じハードウェアと同等の解像度を用いても、Gadget-2の8倍以上も高速であった。
• タスクの取得および管理のオーバーヘッドは、合計実行時間の3.5%未満であり、タスクベースのランタイムシステムの実行時コストが最小限であることが示された。
• 主な性能損失は、自己およびペア相互作用の計算に起因しており、コア数の増加に伴い有効メモリ帯域幅が低下するためと考えられる。
• Sedov爆発およびSodショックテストのシミュレーション結果は、解析解とよく一致しており、正しさと精度が確認された。
• 性能向上の要因は、定数時間の近接粒子探索、ソートによる重複距離計算の低減、対称性の活用であり、いずれも明示的なSIMDベクタ化なしに達成された。