[論文レビュー] Hierarchical interpolative factorization for elliptic operators: integral equations

Kenneth L. Ho, Lexing Ying|arXiv (Cornell University)|Jul 10, 2013

Electromagnetic Scattering and Analysis被引用数 4

ひとこと要約

本稿では、密な離散化された積分作用素の再帰的で木構造的な因子分解により、楕円型偏微分方程式を効率的に解くための、積分方程式における階層的補間因子分解（HIF-IE）を提案する。補間分解（ID）を用いた行列のスパース化と再帰的次元削減を組み合わせることで、推定された準線形計算量を持つ近似一般化LU分解を実現し、高速多重極法、直接解法、または preconditioner として作用素とその逆作用素の高速応用を可能にする。

ABSTRACT

This paper introduces the hierarchical interpolative factorization for integral equations (HIF-IE) associated with elliptic problems in two and three dimensions. This factorization takes the form of an approximate generalized LU decomposition that permits the efficient application of the discretized operator and its inverse. HIF-IE is based on the recursive skeletonization algorithm but incorporates a novel combination of two key features: (1) a matrix factorization framework for sparsifying structured dense matrices and (2) a recursive dimensional reduction strategy to decrease the cost. Thus, higher-dimensional problems are effectively mapped to one dimension, and we conjecture that constructing, applying, and inverting the factorization all have linear or quasilinear complexity. Numerical experiments support this claim and further demonstrate the performance of our algorithm as a generalized fast multipole method, direct solver, and preconditioner. HIF-IE is compatible with geometric adaptivity and can handle both boundary and volume problems. MATLAB codes are freely available.

研究の動機と目的

• 2次元および3次元における、離散化された楕円型積分方程式から生じる密行列線形系の高い計算コストに対処すること。
• 高対比または複雑な幾何形状において、従来の直接解法（O(N³)）や反復法（病的な条件数の問題では収束が遅い）の限界を克服すること。
• 幾何的適応性を備え、境界積分方程式および体積積分方程式の両方に対してスケーラブルで、堅牢かつ高精度な直接解法を開発すること。
• 因子分解、作用素の応用、逆作用素の計算において、線形または準線形の計算量を達成し、時間積み上げ法、逆問題、複数の右辺に対する問題への応用を可能にすること。
• 単一の因子分解構造を通じて、高速多重極法、直接解法、preconditioner と互換性を持つ統一的なフレームワークを提供すること。

提案手法

• 楕円型積分方程式から得られる密行列系に対する、階層的かつ木構造的な補間因子分解であるHIF-IEを提案する。
• 補間分解（ID）に基づく新しい行列スパース化フレームワークを組み込んだ再帰的スケルトン化アルゴリズム（RS）を用い、非対角ブロックの圧縮を実現する。
• 高次元問題を低次元の部分問題に写像する再帰的次元削減を導入し、計算コストを低減する。
• 因子分解を階層的木構造に格納することで、効率的な行列-ベクトル乗算を可能にする、近似一般化LU分解を構築する。
• スケルトン化作用素を用いて低ランク構造を適応的に特定し、IDによる基底ベクトルの選択により圧縮効率と精度を最適化する。
• 幾何的適応性をサポートし、ラプラス、ヘルムホルツ、および他の楕円型偏微分方程式に由来する特異核を含む、第1種および第2種のフレドホルム積分方程式の両方を処理できる。

実験結果

リサーチクエスチョン

• RQ1 階層的補間因子分解は、2次元および3次元の楕円型偏微分方程式から生じる密な積分方程式に対して、準線形計算量を達成できるか？
• RQ2 補間分解による行列スパース化と再帰的次元削減の組み合わせが、直接解法の効率性とスケーラビリティをどのように向上させるか？
• RQ3 HIF-IEは、異なる問題タイプにおいて一貫した性能を示しながら、高速多重極法、直接解法、preconditioner の3つとしての役割を果たせるか、その範囲はどの程度か？
• RQ4 HIF-IEの実験的および理論的計算量は何か？ また、H行列やFMMなど既存手法と比較してどうなるか？
• RQ5 特に高精度において、3次元でスケルトンサイズを顕著に削減しつつも、HIF-IEは高い精度を維持できるか？

主な発見

• HIF-IEは、離散化作用素の因子分解、応用、逆作用素の計算において、推定された準線形計算量を達成しており、数値実験により近似線形スケーリングが確認されている。
• 標準的な再帰的スケルトン化（RS）やrskelf3と比較して、特に高精度な状況下で速度と精度の両面で優れている。hifie3xは、より大きなスケルトンサイズを犠牲にすることで、完全な精度を回復する。
• 数値実験では、HIF-IEが問題サイズに対して良好にスケーリングしており、実験的計算量は理論的推定値とよく一致しているが、漸近的収束はまだ達成されていない。
• HIF-IEは境界積分方程式および体積積分方程式の両方を効果的に処理でき、病的な条件数の第1種問題や良好な条件数の第2種問題を含む。
• アルゴリズムは幾何的適応性に対応でき、複数の右辺や低ランク更新を効率的に処理でき、時間積み上げ法や逆問題に適している。
• MATLAB実装は公開されており、科学計算分野への再現性と広範な採用を促進する。