[論文レビュー] $\mathcal{H}$-matrix approximability of inverses of discretizations of the fractional Laplacian

Michael Karkulik, Jens Markus Melenk|arXiv (Cornell University)|Aug 13, 2018

Advanced Mathematical Modeling in Engineering参考文献 1被引用数 3

ひとこと要約

本稿は、有界領域上での分数ラプラシアンのガレルキン離散化から得られる剛性行列の逆行列が、ブロックランクに関して指数的収束率で階層行列（$\mathcal{H}$-行列）によって近似可能であることを証明する。この手法は、Caffarelli-Silvestre拡張と特異値分解によるブロック別低ランク近似に依存し、$\mathcal{H}$-行列算術を用いた分数PDEの効率的解法について理論的・数値的根拠を提示する。

ABSTRACT

The integral version of the fractional Laplacian on a bounded domain is discretized by a Galerkin approximation based on piecewise linear functions on a quasi-uniform mesh. We show that the inverse of the associated stiffness matrix can be approximated by blockwise low-rank matrices at an exponential rate in the block rank.

研究の動機と目的

• 分数ラプラシアンのガレルキン離散化から得られる剛性行列の逆行列が$\mathcal{H}$-行列によって近似可能である理論的基盤を確立すること。
• 従来、古典的境界要素行列に用いられてきた$\mathcal{H}$-行列技術を、非局所的かつ完全に密な分数微分方程式から生じるシステムへと拡張すること。
• 逆行列の$\mathcal{H}$-行列圧縮を通じて、分数PDEの文脈における効率的反復ソルバーやプリコンディショナの基盤を提供すること。
• 異なるドメインと分数量 $s \in \{0.25, 0.5, 0.75\}$ に対して理論的予測を数値的に検証し、実際の指数的収束率を示すこと。

提案手法

• 準均一な三角形分割上での一次連続有限要素法を用いた分数量ラプラシアンの積分形の離散化。
• 分数量ラプラシアンを$\mathbb{R}^{d+1}$ 上の退化楕円型PDEのディリクレ・ツー・ノイマン写像として表現するためのCaffarelli-Silvestre拡張の適用。
• 幾何的クラスタリングを用いたクラスターツリーおよびブロッククラスターツリーの構築、同定パラメータ $\eta=2$。
• 適切なブロックに対して特異値分解を用いて逆行列のブロック別低ランク近似を計算。
• 標準的な$\mathcal{H}$-行列構造と同定基準（2.11）を用い、データスパース表現を保証。
• 正方形およびL字型ドメインにおける数値実験を、分数量 $s \in \{0.25, 0.5, 0.75\}$ およびメッシュサイズの変化を伴って実施し、指数的収束率の妥当性を検証。

実験結果

リサーチクエスチョン

• RQ1 ガレルキン離散化から生じる分数量ラプラシアンの剛性行列の逆行列は、ブロックランクに関して指数的収束率で$\mathcal{H}$-行列によって近似可能か？
• RQ2 Caffarelli-Silvestre拡張は、分数量ラプラシアンの逆行列に対する$\mathcal{H}$-行列近似可能性を理論的に裏付けることができるか？
• RQ3 $\mathcal{H}$-行列近似の収束速度は、解析によって予測された理論的上限と比較してどの程度か？
• RQ4 $\mathcal{H}$-行列形式を$\mathcal{H}^2$-行列や$\mathcal{H}$-LU分解へと分数量ラプラシアンに拡張可能か？
• RQ5 異なるドメインおよび分数量における$\mathcal{H}$-行列近似の実効的性能はいかがなものか？

主な発見

• 分数量ラプラシアンの剛性行列の逆行列は、ブロックランクに関して指数的収束率で$\mathcal{H}$-行列近似可能であり、理論的誤差境界 $\|\mathbf{A}^{-1} - \mathbf{B}_{\mathcal{H}}^{r}\|_2 \lesssim e^{-br^{1/4}}$ を有する。
• 数値実験では理論的予測を上回るより速い経験的収束率が観測され、誤差は $\|\mathbf{A}^{-1} - \mathbf{B}_{\mathcal{H}}^{r}\|_2 \sim e^{-10r^{1/3}}$ のように減少する。
• 正方形およびL字型ドメインを含むさまざまなドメイン、および分数量 $s = 0.25, 0.5, 0.75$ に対して、本手法は安定である。
• $\mathcal{H}$-行列近似による逆行列表現は、対数線形計算量の算術を可能とし、プリコンディショニングまたは直接法による分数量PDEの効率的解法を支援する。
• 解析は、固定次数 $p$ の多項式による高次元有限要素（分離可能な多項式）へと一般化可能であり、$\mathcal{H}^2$-行列や$\mathcal{H}$-LU分解へも拡張可能である。
• 数値結果は、$\mathcal{H}$-行列形式が分数量ラプラシアンの逆行列に対して計算的に実用可能であり、中程度のブロックランクでも誤差が急速に減少することを確認している。