[論文レビュー] Galerkin--Chebyshev approximation of Gaussian random fields on compact Riemannian manifolds

Annika Lang, Mike Pereira|arXiv (Cornell University)|Jul 6, 2021

Soil Geostatistics and Mapping参考文献 57被引用数 4

ひとこと要約

本稿では、ラプラシアン–ベルトラミ作用素のスペクトル分解とチェビシェフ多項式近似を用いて、コンパクトなリーマン多様体上でのガウス型確率場の効率的サンプリングを可能にするガラーキン–チェビシェフ近似法を提案する。この手法はメッシュサイズに対して線形の計算コストと記憶コストを達成し、ガラーキン近似およびガラーキン–チェビシェフ近似の両方で最適な収束速度を示す。球面および双曲面における数値的検証も行われている。

ABSTRACT

A new numerical approximation method for a class of Gaussian random fields on compact connected oriented Riemannian manifolds is introduced. This class of random fields is characterized by the Laplace--Beltrami operator on the manifold. A Galerkin approximation is combined with a polynomial approximation using Chebyshev series. This so-called Galerkin--Chebyshev approximation scheme yields efficient and generic sampling algorithms for Gaussian random fields on manifolds. Strong and weak orders of convergence for the Galerkin approximation and strong convergence orders for the Galerkin--Chebyshev approximation are shown and confirmed through numerical experiments.

研究の動機と目的

• 球以外のコンパクトなリーマン多様体上でのガウス型確率場のサンプリングに一般的で効率的かつスケーラブルな手法を開発すること。
• ガラーキン近似とチェビシェフ多項式展開を組み合わせることで、有限次元表現における相関係数の高速計算を可能にすること。
• ガラーキン近似の強い収束および弱い収束速度、およびガラーキン–チェビシェフスキームの強い収束速度を確立すること。
• モンテカルロ共分散推定を用いて、球面および双曲面の両方で理論的収束速度を数値的に検証すること。
• 2次元多様体上でのマテルン型場のε-精度でのサンプリングにかかる計算複雑性を𝒪(ε⁻²⁄ρ)に抑えること。

提案手法

• 本手法は、コンパクトなリーマン多様体ℳ上のガウス型確率場𝒁を、ラプラシアン–ベルトラミ作用素−Δℳの関数を含む級数展開によって構築する。
• ガラーキン法を用いて𝒁を有限次元空間Vₙ（例えば線形有限要素空間）に射影することで、ガラーキン近似ℤₙを形成する。
• ℤₙの基底表現における係数ベクトルの共分散行列を閉形式で導出でき、相関係数の明示的サンプリングが可能になる。
• 行列-ベクトル積による場の実現値の高速計算を実現するため、スペクトル関数のチェビシェフ多項式近似を適用する。
• マス・ランプイングと多項式の根の特定を活用することで、収束速度を維持したまま効率性を保つ。
• 理論的収束は平均二乗誤差および共分散誤差の観点から分析され、数値実験により予測された収束率が確認されている。

実験結果

リサーチクエスチョン

• RQ1 一般のコンパクトなリーマン多様体（球に限らない）におけるガウス型確率場の一般的で効率的なサンプリングアルゴリズムを開発可能か？
• RQ2 平均二乗誤差および共分散誤差の観点から、このような場におけるガラーキン近似の収束速度は何か？
• RQ3 スペクトル関数のチェビシェフ近似は、場のサンプリングにおける精度と計算コストにどのように影響するか？
• RQ4 提案されたガラーキン–チェビシェフスキームは、自由度の数に対して線形スケーリングを達成しつつ、最適な収束速度を維持できるか？
• RQ5 理論的収束速度は、球面以外の多様体（例えば双曲面）に対しても実際の応用で成立するか？

主な発見

• 標準的な仮定の下で、ガラーキン近似ℤₙは真の場𝒁に対して平均二乗誤差収束速度𝒪(n⁻¹)で収束する。
• ℤₙと𝒁間の共分散誤差は2αβ − 1 = νの速度で収束し、ウィットル–マテルン場における理論的予測を確認する。
• 球面上では、マテルン場の滑らかさパラメータνに対して、ガラーキン–チェビシェフ近似はルート平均二乗誤差収束速度𝒪(n⁻ν)を達成する。
• 2次元多様体では、精度ε > 0の場のサンプリングにかかる計算コストは𝒪(ε⁻²⁄ρ)に比例し、マテルン場ではρ = νである。
• 球面および双曲面上の数値実験により、理論的収束速度が確認された：ν = 1の場合、共分散誤差は𝒪(n⁻¹)で減少する。
• 精度|cₖ/cₘₐₓ| < 10⁻¹²を達成するための多項式近似次数Kは低く保たれ、高Kに対しても1サンプルあたりの計算時間は0.5秒未満にとどまる。