A Quaternion--BCH Framework for the Local Accuracy of SIDER Interpolation

技術要約：SIDER補間における局所精度のためのクォータニオン–BCHフレームワーク

問題提起
本論文は、方向、方位、剛体運動を扱う科学計算において一般的な課題である、単位球面 $S^2$ 上に制約されたデータの補間について取り組んでいる。標準的なユークリッド多項式補間は、球面制約を維持できない。また、結果を再び球面へ投影する手法は、幾何学的性質や近似特性を変化させてしまう。球面線形補間（SLERP）は、単位制約を維持し、幾何学的に自然な手法であるが、区分的なSLERP曲線は滑らかさに欠け、精度が低いという課題がある。先行研究 [1] で導入されたSIDER（Spherical Interpolation of orDER $n$）法は、アフィンステップをSLERPに置き換えることで、ベジェおよびネヴィル型の補間に着想を得た再帰的な高次構成を提供する。しかし、SLERP操作に固有の曲率依存の補正項が存在するため、滑らかなデータに対するSIDERの局所精度次数 $O(h^{n+1})$ の厳密な証明は困難であった。

手法
本論文は、SIDERの局所的な一貫性を分析するために、クォータニオン–リー代数フレームワークを開発している。核心となる手法の転換は、点を純単位クォータニオン $q = (0, p)$ として表現し、固定基底クォータニオン対数座標において補間誤差を測定することにある。

1. クォータニオン表現: データ点 $p \in S^2$ は、純単位クォータニオン $q = (0, p)$ として埋め込まれる。リー代数 $\text{Im } \mathbb{H}$ は $\mathbb{R}^3$ と同一視され、そこでのリー括弧は外積の2倍に対応する。
2. 厳密なBCH表現: 左基本（left-trivialized）対数座標において、SLERP操作はベーカー–キャンベル–ハウスドルフ（BCH）公式を用いた厳密な表現が可能であるという洞察：
   $$\log(q_\theta^{-1} \text{SLERP}(q_\theta e^X, q_\theta e^Y, \lambda)) = \text{BCH}(X, \lambda \text{BCH}(-X, Y))$$
   これにより、幾何学的な補間問題が、フィルタリングされたリー多項式展開を用いた代数問題へと変換される。
3. フィルタリングされた分割誤差の計算: 証明には、古典的な分割誤差（divided-error）議論の「フィルタリングされた」類似手法が用いられる。誤差はステップサイズ $h$ の形式的冪級数として展開され、係数はシフトされたステンシル変数 $\xi = \theta - i$ の多項式となる。このフィルタリングは、べき乗に対する多項式の次数を追跡する。

主要な貢献と構造的結果
本解析により、このフレームワーク内におけるSLERP操作の2つの重要な構造的特性が明らかになった。

• アフィン一次項: SLERP写像は、対数座標において一次項に関してアフィンである。
• 二次補正の不在: 展開には二次項が含まれない。
• 三次交換子補正: 最初の非線形補正は、三次であり、交換子値を持つ（$[Y-X, [X, Y]]$ に比例する）。純虚クォータニオンの交換子は外積に対応するため、この項は終点の誤差が共線的な場合に消失する。

主要な結果
本論文は、「スモール・ステンシル」仮定（すべての点が共通の正規近傍内に留まることを保証する）の下で、等間隔のパラメータ値でサンプリングされた滑らかな球面曲線に対するSIDER-$n$ の局所精度を証明する。

1. ベースケース (SIDER2): 3つのデータ点から構成されるSIDER2補間子は、3次局所誤差（$O(h^3)$）を示す。証明によれば、主要な誤差項は、ユークリッド空間における通常の二次補間と同じである。三次BCH補正は、最低次の対数データが曲線の接ベクトルに対する共線的な倍数であるため、交換子項がゼロとなり、この次数で消失する。
2. 再帰ステップ: 本論文は、次数を上げる再帰的補題を確立する。隣接する2つのSIDER-$(k-1)$ 候補が「フィルタリングされた分割誤差形式」（すなわち、$O(h^k)$ かつ特定のノード因子を持つ）の誤差を持つ場合、次のSIDER-$k$ ステップは次数を $O(h^{k+1})$ へと引き上げる。
   • アフィン部分の再帰は、ネヴィル型の差分消去を実行し、余分な因子 $h$ と新しいノード積 $\Pi_k(\xi)$ を生成する。
   • 非線形BCH剰余は、鋭いフィルタリング次数構造を保持することが示される。対角変数で表現した場合、非線形項は入力間の差の十分なべき乗を含んでおり、補間パラメータのべき乗を補償して、誤差が要求される多項式次数内に留まることを保証する。
3. 一般定理: すべての固定された $n \ge 2$ に対して、SIDER-$n$ 再構成は以下を満たす：
   $$d_{S^2}(\gamma(\theta h), P_i^{[n]}(\theta; h)) = O(h^{n+1})$$
   証明は、ネヴィル型の消去が持続することを可能にする代数的適合条件として、主要な分割誤差係数のシフト不変性を特定している。

意義と範囲
本論文は、BCH展開を活用することで、曲率依存の項の問題を解決し、SIDER法の期待される局所精度に対する厳密かつ代数的な証明を提供する。これは、球面の非可換幾何学が、データが滑らかであり、かつステンシルが小さい限り、再帰的補間の主要な精度を低下させないことを明確にしている。

解析は厳密に局所的かつ漸近的である。以下の前提に基づいている：

• 等間隔のパラメータ値。
• データ点および中間値が共通の正規近傍内に存在すること（対蹠点の曖昧さや大きな角度ジャンプを回避する）。
• 基礎となる曲線の滑らかさ。

本論文は、滑らかなSIDER補間演算子を分析していることを明示しており、SENO（Spherical Essentially Non-Oscillatory）のステンシル選択手順における大域的な安定性、単調性、または非振動的挙動については扱っていない。ただし、使用される滑らかな候補の整合性の基礎となる理論を提供している。本フレームワークは現在、一様ステンシルに限定されており、非一様ステンシルへの拡張には、非一様重みを用いた分割誤差計算の再定式化が必要となる。

注：本論文は、計算モデリングおよび数学的導出がGPT-5.5 (Thinking) の支援を受けて行われたことを開示しており、その知見はAI支援による仮説として提示され、将来的な形式的な査読検証の対象となる。