[論文レビュー] Solving Elliptic Finite Element Systems in Near-Linear Time with Support Preconditioners

Erik G. Boman, Bruce Hendrickson+1 more|ArXiv.org|Jul 9, 2004

Matrix Theory and Algorithms参考文献 14被引用数 5

ひとこと要約

本稿では、サポート理論を用いて剛性行列を対角優勢行列で近似することにより、スカラー楕円型PDEの有限要素離散化に対する近線形時間解法を提案する。主な貢献は、プリコンディショニング系の条件数が問題のサイズや領域形状に依存せず、メッシュの品質と問題パラメータにのみ依存することを証明したことである。これにより、サポートプリコンディショナーを用いた効率的な反復解法が可能になる。

ABSTRACT

We consider linear systems arising from the use of the finite element method for solving scalar linear elliptic problems. Our main result is that these linear systems, which are symmetric and positive semidefinite, are well approximated by symmetric diagonally dominant matrices. Our framework for defining matrix approximation is support theory. Significant graph theoretic work has already been developed in the support framework for preconditioners in the diagonally dominant case, and in particular it is known that such systems can be solved with iterative methods in nearly linear time. Thus, our approximation result implies that these graph theoretic techniques can also solve a class of finite element problems in nearly linear time. We show that the support number bounds, which control the number of iterations in the preconditioned iterative solver, depend on mesh quality measures but not on the problem size or shape of the domain.

研究の動機と目的

• スカラー楕円型PDEの有限要素離散化から生じる大規模なスパース線形系に対する高速でスケーラブルな反復解法の開発を目的とする。
• 従来、対角優勢行列に対して有効であったサポート理論に基づくプリコンディショニングを、有限要素剛性行列へ拡張することを目的とする。
• プリコンディショニング系の条件数に対する厳密な境界を、問題のサイズや領域形状に依存しない形で確立することを目的とする。
• 近似の品質が、導電率の変動や領域形状に依存せず、メッシュ品質の指標、要素種別、多項式次数にのみ依存することを示すこと。

提案手法

• 剛性行列 $ K $ を $ K = A^T \bar{D}^{1/2} H \bar{D}^{1/2} A $ と因数分解する。ここで $ A $ は要素からグローバルへの接続行列、$ \bar{D} $ は要素のガウス求積重みの対角行列である。
• 行列 $ H $ を単位行列に置き換えることで、対角優勢行列 $ \bar{K} = A^T \bar{D} A $ を構築する。これにより、対角優勢系に既存のサポート理論を適用可能にする。
• 近似の品質はサポート数 $ \sigma(\bar{K}, K) $ を用いて分析され、これは $ \kappa(H) $ で上限が与えられる。$ \kappa(H) $ はメッシュ品質と要素種別にのみ依存する。
• 解析により、$ \kappa(H) $ が要素数、メッシュサイズ、領域形状に依存しないことが示され、要素行列の条件数およびアスペクト比などのメッシュ品質指標にのみ依存することが判明した。
• サポートツリーを用いた既存の対角優勢系に対するほぼ線形時間解法が、近似行列 $ \bar{K} $ に適用可能である。
• 計算実験により、条件数の理論的境界が妥当であることが検証され、さまざまなメッシュタイプと要素次数において、本手法の実用的スケーラビリティが示された。

実験結果

リサーチクエスチョン

• RQ1 スカラー楕円型PDEの有限要素剛性行列は、近線形時間解法を保ちつつ、対角優勢行列でよく近似可能か？
• RQ2 プリコンディショニング系の条件数は、問題のサイズや領域形状に依存せず、メッシュ品質にのみ依存するか？
• RQ3 サポート理論は、従来の対角優勢行列に限らず、標準的なFEM離散化から生じる有限要素剛性行列へ拡張可能か？
• RQ4 要素のアスペクト比や要素行列の条件数といったメッシュ品質指標は、近似の品質と解法収束性にどのように影響するか？
• RQ5 提案されたプリコンディショナーは、導電率分布や要素種別が変化する大規模問題に対しても、効率的に実装・スケーリング可能か？

主な発見

• スカラー楕円型PDEの有限要素離散化における剛性行列 $ K $ は、対角優勢行列 $ \bar{K} $ で近似可能であり、サポート数 $ \sigma(\bar{K}, K) $ は $ \kappa(H) $ で上限が与えられる。$ \kappa(H) $ はメッシュ品質と要素種別にのみ依存する。
• プリコンディショニング系の条件数は、要素数、メッシュサイズ、領域形状に依存せず、空間次元 $ d $、多項式次数 $ p $、求積則、メッシュ品質指標にのみ依存する。
• テスト事例において、最悪メッシュでもプリコンディショニング系の条件数は20.9未満に抑えられ、反復解法の収束性が強く保証された。
• 導電率場 $ \theta $ が0および無限大から離れている限り、サポート数の境界は $ \theta $ に依存せず、メッシュが形状正則である限り成立する。
• 計算結果から、三角形メッシュおよび四面体メッシュを含むさまざまなメッシュタイプと多項式次数において、近似品質が安定しており、メッシュ歪みに伴い条件数が僅かに増加するにとどまった。
• 既存の対角優勢行列用サポートツリー基盤解法を用いることで、有限要素系のほぼ線形時間解法が可能となり、その適用範囲が広範な楕円型PDEにまで拡張された。