[論文レビュー] An inverse theorem for compact Lipschitz regions in $R^d$ using localized kernel bases

Thomas Hangelbroek, Francis J. Narcowich+2 more|arXiv (Cornell University)|Aug 12, 2015

Numerical methods in engineering参考文献 26被引用数 4

ひとこと要約

本稿は、$\mathbb{R}^d$ 内のコンパクトなリプシッツ領域上で、マテRNおよびサーフェススプラインRBFから導かれる局所化されたカーネル基底を用いて、径数関数（RBF）近似の逆推定を確立する。Sobolevノルムが、充填距離 $h$ に依存する $h^{-\sigma}$ スケーリングを施した $L_p$ ノルムによって有界であることを証明し、メッシュレス法における長年のギャップを解消し、PDEソルバーにおけるノルム安定性および条件数の上限を可能にする。

ABSTRACT

While inverse estimates in the context of radial basis function approximation on boundary-free domains have been known for at least ten years, such theorems for the more important and difficult setting of bounded domains have been notably absent. This article develops inverse estimates for finite dimensional spaces arising in radial basis function approximation and meshless methods. The inverse estimates we consider control Sobolev norms of linear combinations of a localized basis by the $L_p$ norm over a bounded domain. The localized basis is generated by forming local Lagrange functions for certain types of RBFs (namely Matérn and surface spline RBFs). In this way it extends the boundary-free construction of Fuselier, Hangelbroek, Narcowich, Ward and Wright.

研究の動機と目的

• メッシュレス法における重要なギャップを埋めるために、有界でコンパクトなリプシッツ領域上での径数関数近似に対する逆推定を開発すること。
• 境界のない設定でのみ知られていた局所化RBF基底の理論を、実用的なPDE応用にとってより関連性の高い境界を有する領域へと拡張すること。
• RBF空間に属する関数のSobolevノルムが、充填距離 $h$ に依存する $h^{-\sigma}$ スケーリングを施した $L_p$ ノルムによって制御されることを証明し、安定性と収束を保証すること。
• ガレルキン法に不可欠な $L_p$ 安定性およびノルム同値性の結果を提供すること、特に混合有限要素法および鞍点問題において重要である。
• 境界付き領域設定におけるトレース推定およびベルンシュタイン＝ニコルスキー不等式を確立することであり、これらはRBFベースの手法においてこれまで得られていなかった。

提案手法

• マテルンおよびサーフェススプラインカーネルを用いて、コンパクトなリプシッツ領域上での局所化されたRBF基底を、局所ラグランジュ関数を用いて構築する。
• ユニット分割と、半径 $\gamma q$ の球 $B_\xi$ 上に台を持つ局所化カーネル $\chi_\xi$ を用いる。ここで $\gamma$ は $\rho$、$d$、$m$、$p$ に依存する。
• 非局所的基底関数間の相互作用を制御するために、準三角不等式と指数的減衰推定を適用する。
• K-関数法と $\ell_p$ ノルム推定を用いて、$\|s\|_{W_p^\sigma(\Omega)}$ を $\|s\|_{L_p(\Omega)}$ の項で上限付ける。
• 基底の $L_p$ 安定性を証明し、$\|s\|_{L_p(\Omega)} \geq c_1 q^{d/p} \|\mathbf{a}\|_{\ell_p(\Xi)}$ を得る。ここで $c_1 \geq C \rho^{\frac{d(d-2m)}{2\epsilon p}} (\log \rho)^{-\frac{d^2}{p\epsilon}}$ である。
• 基底関数間の相互作用を二重指数的分解し、RBFカーネルからの指数的減衰を用いて非対角項を上限付ける。

実験結果

リサーチクエスチョン

• RQ1 有界なリプシッツ領域上でのRBF空間に対して、$\|s\|_{W_p^\sigma(\Omega)} \leq C h^{-\sigma} \|s\|_{L_p(\Omega)}$ の形の逆推定を確立できるか。
• RQ2 ガレルキン法における剛性行列の $L_p$ 安定性と非対角項の減衰を保証するために、どのように局所化RBF基底を構築できるか。
• RQ3 逆推定定数が充填距離 $h$、次元 $d$、滑らかさ $\sigma$ にどのように依存するか。
• RQ4 境界付き領域設定におけるRBFベースの手法において、トレース推定およびベルンシュタイン＝ニコルスキー不等式を導出できるか。これは境界のない場合にのみ既知であった。
• RQ5 境界のない多様体から得られる局所化基底フレームワークを、境界を有する領域にどの程度適応できるか。

主な発見

• 本稿は、コンパクトなリプシッツ領域上でのRBF空間に属する関数に対して、$\|s\|_{W_p^\sigma(\Omega)} \leq C h^{-\sigma} \|s\|_{L_p(\Omega)}$ の形の逆推定を確立する。ここで $C$ は点集合に依存しないが、$d$、$p$、$\sigma$ に依存する。
• $L_p$ 安定性は、$\|s\|_{L_p(\Omega)} \geq c_1 q^{d/p} \|\mathbf{a}\|_{\ell_p(\Xi)}$ を満たし、$c_1 \geq C \rho^{\frac{d(d-2m)}{2\epsilon p}} (\log \rho)^{-\frac{d^2}{p\epsilon}}$（$0 < \epsilon < m - d/2$）であることを証明した。
• 逆推定は、ガレルキン法における剛性行列の条件数が $h^{-2\sigma}$ のオーダーでしか増大しないことを示し、数値的安定性を保証する。
• トレース推定およびベルンシュタイン＝ニコルスキー不等式は、逆推定の直接的結果として得られ、有界領域への適用が可能になる。
• マテルンおよびサーフェススプラインRBFに対して、局所ラグランジュ関数を用いた局所化基底の構築が、境界の存在下でも有効かつ安定であることが示された。
• 本結果は、有界領域におけるメッシュレス法の文脈で、長年にわたり欠落していた逆推定の欠落を解消し、特にRBFガレルキン法およびコロケーション法において重要である。