[論文レビュー] A two-grid accelerated sweeping preconditioner for the Helmholtz equation

Christiaan C. Stolk|arXiv (Cornell University)|Dec 1, 2014

Electromagnetic Scattering and Analysis参考文献 13被引用数 2

ひとこと要約

本稿では、ヘルムホルツ方程式に対する2グリッド高速化スイーピング予約法を提案する。この手法は、スイーピング予約法を8分の1のスケールに粗化された問題にのみ適用することで、計算コストを低減する。新規の粗スケール作用素と並列化戦略を組み合わせることで、1波長あたりの格子点数が少ない場合でも高速収束を実現し、10^8自由度にまで達する大規模問題において、最先端の性能と顕著なコスト削減を達成する。

ABSTRACT

Helmholtz solvers based on sweeping preconditioners have been applied successfully in large 3-D examples. However, such solvers are relatively expensive, due to the high cost per iteration. We propose to reduce this cost by using the sweeping preconditioner at the coarse level of a two-grid method, so that it is applied to an eight times smaller problem. For this to be attractive, the two-grid method must converge rapidly, even when the discretization is done using few points per wavelength. new coarse scale operator was developed to accomplish this, using the ideas of [C.C.Stolk, M. Ahmed and S.K. Bhowmik, A multigrid method for the Helmholtz equation with optimized coarse grid corrections, DOI:10.1137/13092349X]. Several strategies to parallellize the algorithm are proposed. Experiments with a 3-D, MPI implementation show state-of-the-art performance, and a large cost reduction compared to pure sweeping methods for problems with up to 10^8 degrees of freedom.

研究の動機と目的

• 大規模3次元問題におけるヘルムホルツソルバーにおけるスイーピング予約法の高い1反復あたりの計算コストを低減すること。
• 1波長あたりの格子点数が少ない場合でも、2グリッド法の収束を迅速に達成できることを実現すること。
• 粗化された問題においても精度と効率を維持できる新しい粗スケール作用素を開発すること。
• 2グリッドアルゴリズムのスケーラビリティを向上させるために並列化戦略を設計すること。
• 標準的なスイーピング予約法と比較して顕著な計算コスト削減を達成すること。

提案手法

• 1反復あたりの計算負荷を軽減するため、スイーピング予約法を自由度が1/8の粗化された問題にのみ適用する。
• マルチグリッド法にインspiredされた最適化された粗グリッド補正を備えた新規の粗スケール作用素を導入し、粗グリッド上での収束を維持する。
• 2グリッドフレームワークを採用し、粗グリッド補正を縮小された問題にスイーピング予約法を適用することで計算する。
• マルチグリッドにインspiredされた粗作用素を活用することで、1波長あたりの格子点数が少ない場合でも、強固で高速な収束を確保する。
• MPIを用いた並列実行により、大規模3次元問題におけるスケーラビリティを実現する。
• 粗グリッドの解法と予約法の演算をプロセッサ間で効率的に分散するための並列化戦略を統合する。

実験結果

リサーチクエスチョン

• RQ1 2グリッド法の粗いグリッドレベルでのみスイーピング予約法を適用することで、全体の計算コストを著しく低減できるか？
• RQ2 1波長あたりの格子点数が少ないヘルムホルツ問題においても、2グリッド法の収束を迅速に達成できるように、新しい粗スケール作用素をどのように設計できるか？
• RQ3 大規模3次元ヘルムホルツ問題において、標準的なスイーピング予約法と比較して、どの程度の性能向上とコスト削減が達成できるか？
• RQ4 大規模問題に対して、MPIを用いた2グリッドスイーピング予約法をどの程度効率的に並列化できるか？
• RQ5 離散化が粗い（1波長あたりの格子点数が少ない）場合でも、提案手法が精度と強健性を維持できるか？

主な発見

• 2グリッド高速化スイーピング予約法は、自由度が10^8に達する3次元ヘルムホルツ問題において、最先端の性能を達成する。
• 高価な予約法を粗化された問題にのみ適用することで、標準的なスイーピング予約法と比較して計算コストを顕著に低減する。
• 新規の粗スケール作用素により、1波長あたりの格子点数が少ない場合でも迅速な収束が実現され、強健性が維持される。
• MPIによる実装は効果的な並列化を示し、大規模シミュレーションをサポートする。
• 大規模3次元問題において、精度とスケーラビリティを保ちながら顕著なコスト削減を達成する。
• 特に高周波数および大規模問題において、純粋なスイーピング予約法と比較して、計算効率が優れている。