[論文レビュー] Compact Finite Differences and Cubic Splines

Robert M. Corless|arXiv (Cornell University)|May 19, 2018

Advanced Numerical Analysis Techniques参考文献 25被引用数 3

ひとこと要約

本稿は、一様メッシュ上における4次精度コンパクト差分と3次スプラインの間の新しい数学的関係を確立し、Padéに基づくコンパクト差分スキームが3次スプラインにおける2回連続微分可能性に必要な正確な導関数を生成することを示している。境界で高次の精度を達成する新しいコンパクト3次補間関数を導入し、一様またはやや非一様なメッシュ上でも、関連する行列が完全非負かつ正定値であることを証明している。

ABSTRACT

In this paper I uncover and explain---using contour integrals and residues---a connection between cubic splines and a popular compact finite difference formula. The connection is that on a uniform mesh the simplest Padé scheme for generating fourth-order accurate compact finite differences gives extsl{exactly} the derivatives at the interior nodes needed to guarantee twice-continuous differentiability for cubic splines. %I found this connection surprising, because the two problems being solved are different. I also introduce an apparently new spline-like interpolant that I call a compact cubic interpolant; this is similar to one introduced in 1972 by Swartz and Varga, but has higher order accuracy at the edges. I argue that for mildly nonuniform meshes the compact cubic approach offers some potential advantages, and even for uniform meshes offers a simple way to treat the edge conditions, relieving the user of the burden of deciding to use one of the three standard options: free (natural), complete (clamped), or "not-a-knot" conditions. Finally, I establish that the matrices defining the compact cubic splines (equivalently, the fourth-order compact finite difference formulæ) are positive definite, and in fact totally nonnegative, if all mesh widths are the same sign.

研究の動機と目的

• 一様メッシュ上における4次精度コンパクト差分スキームと3次スプラインの導関数の間で、これまで注目されなかった数学的関係を明らかにすること。
• 従来の3次スプラインよりも境界で高次の精度を達成する、新しいスプラインに類似した補間関数（コンパクト3次補間関数）を構築すること。
• 自然境界条件、固定境界条件、またはノット・ア・ノット条件といった選択肢をユーザーが選ぶ必要なく、境界条件を統一的に扱うフレームワークを提供すること。
• メッシュ幅が同じ符号を持つ場合に、コンパクト3次スプラインの係数行列が正定値かつ完全非負であることを証明すること。

提案手法

• コンtour積分と留数計算を用いて、コンパクト差分式と3次スプラインに必要な導関数値との間の等価性を解析的に導出し、解釈している。
• Padé近似を用いてコンパクト差分スキームを導出し、内部ノードにおいて4次精度を保証している。
• 境界で4次精度 $O(h^4)$ を達成できるように、4ノードの修正されたコンパクト差分式を導入し、わずかに非一様なメッシュ上でも有効である。
• 得られた区分的3次補間関数の安定的かつ効率的な評価のために、重心形式を適用している。
• メッシュ幅がすべて正またはすべて負の場合に、係数行列が完全非負かつ正定値であることを証明している。
• 新規のコンパクト3次補間関数を古典的スプラインおよびpchipと比較し、ノードにおける導関数の精度向上と境界処理の簡素化に焦点を当てている。

実験結果

リサーチクエスチョン

• RQ1 一様メッシュ上における4次精度コンパクト差分スキームが、なぜ3次スプラインに必要な正確な導関数値を生成するのか？
• RQ2 3次スプラインの滑らかさを保ちつつ、境界で高次の精度を達成する新たな補間関数を構築可能か？
• RQ3 非一様メッシュ下で、コンパクト3次補間関数の安定性と数値的条件数は、古典的スプライン法と比べてどうか？
• RQ4 コンパクト3次スプラインの係数行列が持つ代数的性質（例：定値性）は何か？
• RQ5 コンパクト差分と3次スプラインの間の観察された等価性は偶然に過ぎないのか、それともより深い解析的説明が存在するのか？

主な発見

• 一様メッシュ上では、Padé近似に基づく4次精度コンパクト差分スキームが、3次スプラインに必要な導関数と同一であることが確認され、正確な数学的等価性が裏付けられた。
• 提案されたコンパクト3次補間関数は、端点における導関数について $O(h^4)$ の精度を達成し、古典的スプラインにおける端点条件の $O(h^3)$ の精度を上回っている。
• すべてのメッシュ幅が同じ符号を持つ場合、コンパクト3次スプラインの係数行列が正定値かつ完全非負であることが証明され、数値的安定性が保証された。
• コンパクト差分と3次スプラインの関係は、コンtour積分と留数を用いて解釈され、そうでない驚きの数値的偶然に、複素解析的根拠を与えた。
• スプラインを純粋な3次多項式と $|x - x_0|^3$ 項に分解する別の説明も等価性を説明できるが、留数法がより広範な理論的意義を有するとされる。
• MATLABで実装された `compactcubic.m` を通じて、実用的で安定的かつユーザーフレンドリーな、境界処理を自動で行う従来のスプライン構築法の代替手段が提供された。