[論文レビュー] An algebraic distances measure of AMG strength of connection

Achi Brandt, James Brannick+2 more|arXiv (Cornell University)|Jun 29, 2011

Advanced Numerical Methods in Computational Mathematics参考文献 25被引用数 3

ひとこと要約

本稿では、最小二乗最適化に基づく代数的距離測度を提案し、M行列の仮定に依存せずに、非局所的拡散問題に対しても頑健な粗化と補間を可能にする代数的多重グリッド（AMG）法における接続強度の定義を行う。この手法は、異方性の角度や係数が変化しても一様に有界な収束速度と低複雑性を達成し、グリッドに沿わない問題に対しても頑健であることを示している。

ABSTRACT

Algebraic multigrid is an iterative method that is often optimal for solving the matrix equations that arise in a wide variety of applications, including discretized partial differential equations. It automatically constructs a sequence of increasingly smaller matrix problems that enable efficient resolution of all scales present in the solution. One of the main components of the method is an adequate choice of coarse grids. The current coarsening methodology is based on measuring how a so-called algebraically smooth error value at one point depends on the error values at its neighbors. Such a concept of strength of connection is well understood for operators whose principal part is an M-matrix; however, the strength concept for more general matrices is not yet clearly understood, and this lack of knowledge limits the scope of AMG applicability. The purpose of this paper is to motivate a general definition of strength of connection, based on the notion of algebraic distances, discuss its implementation, and present the results of initial numerical experiments. The algebraic distance measure, we propose, uses as its main tool a least squares functional, which is also applied to define interpolation.

研究の動機と目的

• 古典的AMG法がM行列でない系に適用可能でないという限界を解消すること。特に、非局所的拡散問題に焦点を当てる。
• M行列の性質やヒューリスティックな粗化ルールに依存しない、一般化されたパラメータフリーの接続強度定義を開発すること。
• 任意の異方性方向、特にグリッドに沿わない場合を含む、非局所的問題に対しても効果的な粗グリッド選択と補間を可能にすること。
• さまざまな問題パラメータにわたって、低グリッドおよびオペレータ複雑性を維持しながら、二レベル収束の頑健性を確保すること。

提案手法

• 点間の代数的距離を定義するため、最小二乗関数を用いる。これが接続強度の測度として機能する。
• 同じ最小二乗関数を用いて補間演算子を構築し、接続強度測度と整合性を保つ。
• 代数的距離と適合する緩和法を組み合わせ、代数的に滑らかになる誤差を捉える粗グリッド点を選択する。
• 固定された強度閾値（θ_ad = 0.5）と一様なセットアップパラメータをすべてのテストケースに適用し、頑健性を確保する。
• ガリルキン粗化（A_c = RAP）を用い、代数的距離から導かれた補間Pを用いて粗グリッドレベルの問題を構築する。
• パラメータチューニングを回避し、低複雑性を維持する体系的なセットアップアルゴリズムを採用する。

実験結果

リサーチクエスチョン

• RQ1 最小二乗法で定義された代数的距離は、M行列に依存しない一般化された接続強度測度として有効であるか？
• RQ2 非グリッドに沿った異方性を持つ非局所的拡散問題において、提案手法はどのように性能を示すか？
• RQ3 異方性の角度や係数が変化しても、収束速度とグリッド／オペレータ複雑性が一様に有界に保たれるか？
• RQ4 細粒度グリッドの離散化が異方性方向とずれている場合でも、手法は頑健に機能するか？
• RQ5 この手法は多段階AMGに拡張可能であり、効率性と精度を維持できるか？

主な発見

• 等方的問題（ε = 0, α = 0）では、漸近的二段階収束率がρ = 0.05まで低下し、N = 128でγ_g = 1.3、γ_o = 1.3を達成する。
• ε = 0.1、α = π/4の異方的ケースでは、すべてのメッシュサイズで収束率がρ = 0.35未満を維持し、γ_g ≈ 25、γ_o ≈ 1.5となる。
• グリッドに沿わない異方性（α = −π/4）に対しても、ρ ≤ 0.38、γ_o ≤ 1.9を維持し、方向のずれに対しても頑健であることが示された。
• 固定されたhに対してεおよびαの変動にかかわらず収束率と複雑性が一様に有界であるため、パラメータに敏感でないことが判明した。
• 細粒度グリッドの離散化が異方性方向と逆向きであっても、効果的な粗グリッド演算子が得られ、異方性をよりよく捉えている。
• 初期の細粒度グリッドの向きにかかわらず、適切な離散化と整合する最初の粗グリッド系が得られるため、多段階への拡張が有望である。