[論文レビュー] Generalization of the complex shifted Laplacian: on the class of expansion preconditioners for Helmholtz problems

Siegfried Cools, Wim Vanroose|arXiv (Cornell University)|Jan 19, 2015

Matrix Theory and Algorithms参考文献 25被引用数 3

ひとこと要約

本稿では、逆ヘルムホルツ作用素の級数展開に基づき、ヘルムホルツ問題における複素シフトラプラシアン（CSL）プレコンディショナの一般化として、拡張プレコンディショナクラス EX(m) を導入する。EX(m) における高次項の導入により Krylov 反復法の収束が向上するが、応用コストの増加を鑑み、実用的最適解として EX(1)（CSL に等価）が得られる。

ABSTRACT

In this note we propose the class of expansion preconditioners which forms a direct generalization of the CSL preconditioner. The construction of the EX(m) preconditioner is based upon a truncated Taylor series expansion of the original Helmholtz operator inverse. The EX(m) preconditioner is shown to significantly improve Krylov solver convergence rates for the Helmholtz problem for growing values of m. However, the addition of multiple terms in the expansion also increases the computational cost of applying the preconditioner. A thorough cost-benefit analysis of the addition of extra terms in the EX(m) preconditioner proves the CSL or EX(1) preconditioner to be the practically most efficient member of the EX(m) class. We suggest possible extensions to the expansion preconditioner class which are particularly interesting forma theoretical viewpoint.

研究の動機と目的

• ヘルムホルツ問題における複素シフトラプラシアン（CSL）プレコンディショナを、より広範なプレコンディショナクラスに一般化すること。
• プレコンディショナ展開における項の数を増加させた際の収束性の向上と計算コストのトレードオフを分析すること。
• 実用的効率を最大化する EX(m) クラスの最適なメンバーを特定すること。
• 将来的な研究のための拡張プレコンディショナフレームワークの理論的拡張を検討すること。

提案手法

• 逆ヘルムホルツ作用素のテイラー級数展開に基づき、EX(m) プレコンディショナを構築する。
• 逆作用素の展開から得られる有限項の級数として、プレコンディショナを形式的に定義する。
• Krylov部分空間ソルバー内に EX(m) プレコンディショナを適用し、ヘルムホルツ問題の収束を加速する。
• 異なる m 値におけるソルバーの収束速度とプレコンディショナ適用コストを比較するコスト・ベネフィット分析を実施する。
• CSL プレコンディショナをベースラインとし、EX(m) クラスにおける m=1 のケースに相当する。
• 将来的な発展に向けた拡張プレコンディショナクラスの理論的拡張を検討する。

実験結果

リサーチクエスチョン

• RQ1 切り捨てられたテイラー級数展開における項の数を増加させることで、ヘルムホルツ問題における Krylov 反復法の収束速度にどのような影響を与えるか？
• RQ2 高次 EX(m) プレコンディショナを用いる際、収束性の向上と計算コストの増加のトレードオフはどのように変化するか？
• RQ3 実用的に EX(m) プレコンディショナの効率を最大化する最適な m 値が存在するか？
• RQ4 EX(m) プレコンディショナは、既存の CSL プレコンディショナをどのように一般化するか？
• RQ5 拡張プレコンディショナフレームワークにどのような理論的拡張が可能であり、その応用性を高めうるか？

主な発見

• EX(m) プレコンディショナ族は、逆ヘルムホルツ作用素のテイラー級数展開からの高次項を組み込むことで、CSL プレコンディショナを一般化する。
• m の値を増加させることで、ヘルムホルツ問題における Krylov 反復法の収束が、より良好なプレコンディショニングにより向上する。
• しかし、プレコンディショナの適用コストは m の増加に伴い増加し、特に複数の項を評価する必要があるため顕著になる。
• 収束性の向上は見られるが、高次 EX(m) プレコンディショナの高いコストは、実用的効率を低下させる。
• EX(1) プレコンディショナ（CSL 法に等価）は、実用的に EX(m) クラスの中で最も効率的なメンバーであると判明した。
• 本研究では、将来的な研究に価値ある理論的拡張が可能であると示唆している。