Analytic First Derivatives of SIDER Interpolation

Technische Samenvatting: Analytische Eerste Afgeleiden van SIDER-interpolatie

Probleemstelling
In de numerieke analyse en geometrische modellering is er een groeiende behoefte aan interpolatie van gegevens die beperkt zijn tot niet-lineaire manifolds, specifiek de eenheidssfeer $S^2 \subset \mathbb{R}^3$. Standaard Euclidische polynoominterpolatie slaagt er niet in om de sferische beperking tussen datapunten te behouden, wat vaak puntgewijze normalisatie vereist die de nauwkeurigheid en de analyse van afgeleiden bemoeilijkt. Hoewel Spherical Interpolation of orDER $n$ (SIDER-$n$) eerder werd geïntroduceerd als een hoogwaardig reconstructie-instrument voor Spherical Essentially Non-Oscillatory (SENO) methoden, ontbrak het aan een analytische theorie voor eerste afgeleiden. Het ontbreken van afgeleide-informatie beperkt de toepassing van SIDER in numerieke schema's die lokale hellingen, directionele veranderingen of differentiële grootheden vereisen (bijv. gegeneraliseerde Riemann-solvers, inverse Lax–Wendroff-procedures of hoogwaardige kwadratuur).

Methodologie
Het artikel ontwikkelt analytische formules voor de eerste afgeleiden van SIDER-curven van willekeurige orde door gebruik te maken van de recursieve structuur van de interpolatie. De kernmethodologie berust op de volgende principes:

1. Recursieve Compositie: SIDER-$n$ wordt volledig geconstrueerd uit Spherical Linear Interpolation (SLERP)-operaties, wat een binaire boom van SLERP-afbeeldingen vormt.
2. Kettingregel-propagatie: In plaats van een volledig uitgebreide hoogwaardige expressie te differentiëren, passen de auteurs de kettingregel direct toe via de recursieve definitie. De afgeleide van een SIDER-$n$ curve wordt uitgedrukt als de totale afgeleide van een SLERP-operatie toegepast op twee lagere orde SIDER-$(n-1)$ takken.
3. Totale SLERP-afgeleide: Een fundamenteel bouwblok is afgeleid voor de totale afgeleide van SLERP met bewegende eindpunten en een bewegende interpolatieparameter. Deze formule houdt rekening met de variatie van de geodetische hoek en de eindpunten.
4. Specialisatie: De algemene recursieve formule wordt gespecialiseerd voor:
   • Interne Interpolatienodes: Waar de twee lagere orde takken samenkomen bij hetzelfde datapunt, wat een vereenvoudigde gewogen combinatie van afgeleiden mogelijk maakt.
   • Eindpunten: Waar een tak buiten het standaard interpolatie-interval wordt geëvalueerd (extrapolatie), wat termen met de logaritmische kaart introduceert.
   • Middelpunten: Locaties halverwege de bemonstering, die cruciaal zijn voor hoogwaardige reconstructie-gebaseerde algoritmen.

Belangrijkste Bijdragen

• Recursieve Afleidingsformules: Het artikel biedt een compact, implementatieklaar recursief algoritme voor het berekenen van $D_n(t) = \frac{d}{dt}S_n(t)$ voor elke orde $n$. Dit voorkomt onhandelbare, uitgebreide expressies voor hoogwaardige curven.
• Expliciete Basiskansen: De basisgeval SIDER2 wordt expliciet behandeld, en SIDER3 en SIDER4 worden gebruikt om het recursieve mechanisme en de specifieke afhandeling van eindpunten en interne nodes te illustreren.
• Tangenstiebehoud: Er wordt een rigoureus bewijs geleverd dat de analytisch afgeleide afgeleidevector inherent tangent aan de sfeer is op elk gereconstrueerd punt (zowel bij bemonsteringsnodes als middelpunten). Deze eigenschap ($S_n(t) \cdot D_n(t) = 0$) vloeit voort uit de algebraïsche structuur van de SLERP-afgeleide en vereist geen post-processing projectie.
• Middelpunt-afgeleiden: Het artikel leidt specifieke formules af voor afgeleiden bij half-integer indexcoördinaten ($\xi = j + 1/2$), die essentieel zijn voor het evalueren van reconstructies tussen datapunten in finite-volume en WENO-type schema's.

Resultaten

• Algemene Recursie: De afgeleide $D_n$ wordt berekend door de totale SLERP-afgeleide te evalueren op de twee lagere orde takken $A(t)$ en $B(t)$, geschaald met de herparameteriseringsfactor $n/(n-1)$.
• Node Vereenvoudiging: Bij interne nodes $t_j = j/n$ vereenvoudigt de afgeleide tot een gewogen som van de afgeleiden van de twee lagere orde takken die samenkomen bij $p_j$. De gewichten tellen op tot $n/(n-1)$, wat de schaling van de parametrisatie weerspiegelt.
• Eindpuntgedrag: Bij eindpunten bevat de afgeleide termen die de logaritmische kaart ($\text{Log}_p$) involleren om de richting van de geëxtrapoleerde takwaarden te accounten.
• Geometrische Consistentie: Numerieke visualisaties bevestigen dat de berekende afgeleidevectoren exact in de raum van de sfeer liggen op de geïntripeerde punten, wat de theoretische bewijsvoering van de tangentie valideert.

Betekenis en Reikwijdte
Het artikel stelt dat deze formules de oorspronkelijke SIDER/SENO-framework uitbreiden door de noodzakelijke differentiële informatie te leveren voor sfeer-waardige reconstructies. De betekenis ligt in het mogelijk maken van het gebruik van SIDER in hoogwaardige finite-volume, ENO/WENO- en SENO-type methoden voor behoudswetten en evolutieproblemen waarbij het gereconstrueerde veld gedifferentieerd moet worden terwijl de manifold-beperking strikt behouden blijft. Het werk wordt gepresenteerd als een analytische theorie die de berekening van de afgeleide interpreteert als een forward-mode differentiatie van het SIDER-evaluatiealgoritme.

De auteurs merken op dat de implementatie compact blijft en de recursieve definitie direct volgt. Het artikel concludeert door te erkennen dat hoewel de wiskundige diepgang van de resultaten aanzienlijk is, het werk bedoeld is als een preprint voor verspreiding en toekomstige formele peer-verificatie, met potentieel toekomstig werk betreffende niet-uniforme parametrisatie-afstand en uitbreidingen naar andere manifolds zoals rotatiegroepen.