[論文レビュー] Efficient Numerical Method for Models Driven by Lévy Process via Hierarchical Matrices

Kailai Xu, Eric Darve|arXiv (Cornell University)|Dec 18, 2018

Fractional Differential Equations Solutions参考文献 8被引用数 6

ひとこと要約

本稿では、Lévy過程によって駆動される対流拡散方程式のための効率的な数値解法を、階層行列（H行列）を用いて提案する。H行列による圧縮と、無条件安定なCrank-Nicolsonスキームを活用することで、直接法のO(N²)およびO(N³)のコストに対し、O(N)のストレージおよび計算複雑度を達成し、中〜大規模な問題において顕著な高速化を実現した。

ABSTRACT

Modeling via fractional partial differential equations or a Lévy process has been an active area of research and has many applications. However, the lack of efficient numerical computation methods for general nonlocal operators impedes people from adopting such modeling tools. We proposed an efficient solver for the convection-diffusion equation whose operator is the infinitesimal generator of a Lévy process based on $\mathcal{H}$-matrix technique. The proposed Crank Nicolson scheme is unconditionally stable and has a theoretical $\mathcal{O}(h^2+Δt^2)$ convergence rate. The $\mathcal{H}$-matrix technique has theoretical $\mathcal{O}(N)$ space and computational complexity compared to $\mathcal{O}(N^2)$ and $\mathcal{O}(N^3)$ respectively for the direct method. Numerical experiments demonstrate the efficiency of the new algorithm.

研究の動機と目的

• 非局所的Lévy作用素を伴う対流拡散方程式を解く際の計算ボトルネック、特に密で悪条件な剛性行列に起因する問題を解決すること。
• Lévy過程に由来する分数マージンPDE（特に分数ラプラシアンを含むもの）に対して、高速でスケーラブルな数値法を開発すること。
• 階層行列技術を用いて、一般のLévy駆動モデルの効率的数値解法を可能にし、ストレージおよび計算コストをO(N²)およびO(N³)からO(N)に削減すること。
• 核構造や低ランク展開に関する最小限のユーザー入力で済む、ブラックボックス的で自動的なH行列構築およびLU因子分解を提供すること。

提案手法

• Lévy過程によって駆動される対流拡散方程式に対して、Crank-Nicolson時間離散化を採用し、無条件安定性と理論的O(Δt² + h²)収束を保証する。
• 非局所作用素の密な剛性行列を、非対角ブロックに低ランク構造を有する階層行列（H行列）で表現する。
• 四分木などの階層的分割戦略を用い、再帰的にブロックを分割し、小さなブロックには完全行列を、許容可能なブロックには低ランク近似を格納する。
• 自動的H行列構築を実施：核の展開に関する明示的なユーザー指定を必要とせず、最適なブロック分割と低ランク圧縮をアルゴリズムが自動で決定する。
• H行列算術を用いて、再帰的に行列-ベクトル積およびLU因子分解を実行し、完全行列、低ランク、階層的サブブロックそれぞれに特化したソルバーを適用する。
• 系列展開ベースおよびブラックボックスFMMベースのH行列構築を両方サポートし、ユーザーは核関数または密行列を直接入力可能である。

実験結果

リサーチクエスチョン

• RQ1 一般のLévy測度を有するLévy過程によって駆動される対流拡散方程式に対して、階層行列が効果的に適用可能か？
• RQ2 直接法と比較して、H行列ベースのソルバーの計算およびストレージ複雑度はどの程度達成可能か？
• RQ3 近似的線形複雑度を達成しつつ、2次時間的および空間的収束性を維持できるか？
• RQ4 手動による低ランク展開を必要とせずに、さまざまなLévy測度および核タイプに対して自動H行列構築がどの程度効果的か？
• RQ5 実際の計算において、H行列法が直接LU因子分解を上回る性能を示すのはどのスケールか（クロスオーバー点）？

主な発見

• 提案されたCrank-Nicolsonスキームは無条件安定であり、理論的収束速度O(h² + Δt²)を達成する。
• H行列ベースのソルバーにより、ストレージおよび計算複雑度をO(N²)およびO(N³)からO(N)に削減し、大規模問題のスケーラブルな解法を可能にする。
• 数値実験により、直接法と比較して顕著な高速化が確認され、1次元においてLU因子分解のクロスオーバー点は約N = 1100で発生する。
• 特異的または緩やかに減衰するLévy測度を含む広範なLévy測度に対して、弱い仮定のもとで有効である。
• 自動H行列構築は、系列展開モードおよびブラックボックスFMMモードの両方で機能し、ユーザー入力が最小限で済む。
• 逆剛性行列の近似において、ブロックランクに関して指数的収束を達成しており、低ランク構造の理論的期待と整合的である。