[論文レビュー] Solving the Klein-Gordon equation using Fourier spectral methods: A benchmark test for computer performance

Samar Aseeri, Oleg Batrasev+13 more|arXiv (Cornell University)|Jan 19, 2015

Numerical Methods and Algorithms参考文献 27被引用数 6

ひとこと要約

この論文は、13種類の多様なスーパーコンピューターシステムを対象に、2DECOMP&FFTライブラリを用いて3次元立方体クライン＝ゴードン方程式を解くためのフーリエスペクトル法の性能をベンチマークしている。強スケーリング特性を示し、相互接続遅延およびプロセッサアーキテクチャの違いに起因する性能の顕著な変動を特定しており、HPCシステムおよびPDEソルバーにおける数値手法の実用的かつ現実世界のベンチマークを提供する。

ABSTRACT

The cubic Klein-Gordon equation is a simple but non-trivial partial differential equation whose numerical solution has the main building blocks required for the solution of many other partial differential equations. In this study, the library 2DECOMP&FFT is used in a Fourier spectral scheme to solve the Klein-Gordon equation and strong scaling of the code is examined on thirteen different machines for a problem size of 512^3. The results are useful in assessing likely performance of other parallel fast Fourier transform based programs for solving partial differential equations. The problem is chosen to be large enough to solve on a workstation, yet also of interest to solve quickly on a supercomputer, in particular for parametric studies. Unlike other high performance computing benchmarks, for this problem size, the time to solution will not be improved by simply building a bigger supercomputer.

研究の動機と目的

• 幅広いハイパフォーマンスコンピューティング（HPC）システム上で、3次元立方体クライン＝ゴードン方程式を解くためのフーリエスペクトル法の強スケーリング性能を評価すること。
• 異なるコンピュータアーキテクチャおよび相互接続が、特に512³のような中程度の問題サイズにおけるFFTベースのPDEソルバーの性能に与える影響を評価すること。
• 合成的または合成的と似たワークロードを超える、複雑なPDEシミュレーションの計算的要件を反映する、実用的かつ現実世界のHPCシステムベンチマークを確立すること。
• 非線形PDEシミュレーションの将来のスケーリングに向けた最適なハードウェアおよびソフトウェアスタックの選定を支援すること、特にパrametricスタディで高いスループットを必要とする分野を想定する。
• このベンチマークが、Linpack や NAS Parallel Benchmarks といった既存のHPCランク付けの補完的または更新版として機能できるかどうかを検討すること。

提案手法

• 立方体クライン＝ゴードン方程式は、512³の直交座標グリッド上で、擬似スペクトル法を用いたフーリエスペクトル空間離散化により解かれる。
• 並列FFT処理には2DECOMP&FFTライブラリが用いられ、2次元ブロックサイクリック分割に基づきMPI_AlltoallvまたはMPI_Alltoallを用いたプロセッサ間通信が実行される。
• 13種類の異種HPCシステム（リーダーシップクラスのスーパーコンピュータおよびワークステーションを含む）において強スケーリング実験が実施され、ソリューション時間と効率が測定された。
• 浮動小数点演算回数、通信帯域幅、相互接続遅延の観点から性能を分析し、異なるプロセッサタイプおよびネットワークトポロジー間で比較された。
• スーパーコンピュータにとって十分に大きな問題サイズでありながら、ワークステーションでローカルに後処理が可能なサイズであるため、パrametricスタディが可能である。
• ベンチマークは、時間ステップ処理、非線形項の評価、FFTベースのスペクトル変換を含む、現実世界のPDEシミュレーションワークロードを反映するように設計されている。

実験結果

リサーチクエスチョン

• RQ1 異なる相互接続およびプロセッサタイプを有する多様なHPCアーキテクチャにおいて、3次元クライン＝ゴードン方程式に対するフーリエスペクトル法の強スケーリング性能はどのように変動するか？
• RQ2 現代のスーパーコンピュータにおいて、相互接続遅延およびメモリ帯域幅の違いが、FFTベースのPDEソルバーの性能に及ぼす影響の程度はどの程度か？
• RQ3 このベンチマークは、Linpack や NAS Parallel Benchmarks といった既存のHPCベンチマークの意味のある代替または補完として、現実世界のPDEシミュレーション性能を評価するのに役立つだろうか？
• RQ4 固定問題サイズにおいて、最新のIntelプロセッサは、旧式または代替アーキテクチャ（例：PowerPC、AMD）と比較して、FFTベースのPDEソルバーの性能でどのように差を示すか？
• RQ5 異種HPCシステムにおけるFFTベースのスペクトル法の移植性およびパフォーマンスの移植性について、どのような知見が得られるか？

主な発見

• 同じピーク浮動小数点性能であっても、相互接続遅延および帯域幅の違いにより、フーリエスペクトルソルバーの性能はシステム間で顕著に異なる。
• 同じ計算チップを搭載しているにもかかわらず、VedurはHectorよりも著しく性能が低いことが判明し、相互接続遅延が主要なパフォーマンスボトルネックであることが示された。
• Beacon、Hornet、Marenostrum、Stampedeなどのシステムに搭載された最新のIntelプロセッサは、浮動小数点スループットの高さと優れたメモリコントローラーのおかげで、旧式アーキテクチャ（PowerPC、AMD、レガシーIntel）を上回っている。
• 帯域幅とMPI_Alltoall性能だけでは、観察された性能差を完全に予測できないことから、FLOP/サイクルやメモリコントローラー効率といったマイクロアーキテクチャ的特性の重要性が浮き彫りになった。
• この問題サイズでは、より大きなスーパーコンピュータを構築してもソリューション時間の改善が得られないことが判明し、アルゴリズム的およびアーキテクチャ的選択が支配的になる性能のプラトーに達していることが示された。
• 本研究は、非線形PDEのパラメトリックシミュレーションの計算的要請を反映する、実用的かつ現実世界のHPCシステムベンチマークを提供しており、特にスペクトル法およびFFTに依存する分野に適している。