Analytic First Derivatives of SIDER Interpolation

Resumen Técnico: Primeras Derivadas Analíticas de la Interpolación SIDER

Planteamiento del Problema
En el análisis numérico y el modelado geométrico, existe una necesidad creciente de interpolar datos restringidos a variedades no lineales, específicamente la esfera unidad $S^2 \subset \mathbb{R}^3$. La interpolación polinómica euclídea estándar no logra preservar la restricción esférica entre los puntos de datos, lo que a menudo requiere una normalización puntual que complica la precisión y el análisis de las derivadas. Si bien la Interpolación Esférica de orden $n$ (SIDER-$n$) fue introducida previamente como una herramienta de reconstrucción de alto orden para métodos de Esencialmente No Oscilatorios Esféricos (SENO), carecía de una teoría analítica para las primeras derivadas. La ausencia de información sobre las derivadas limita la aplicación de SIDER en esquemas numéricos que requieren pendientes locales, cambios direccionales o cantidades diferenciales (por ejemplo, resolvedores de Riemann generalizados, procedimientos de Lax–Wendroff inversos o cuadratura de alto orden).

Metodología
El artículo desarrolla fórmulas analíticas para las primeras derivadas de curvas SIDER de orden arbitrario explotando la estructura recursiva de la interpolación. La metodología central se basa en los siguientes principios:

1. Composición Recursiva: SIDER-$n$ se construye íntegramente a partir de operaciones de Interpolación Lineal Esférica (SLERP), formando un árbol binario de mapas SLERP.
2. Propagación de la Regla de la Cadena: En lugar de derivar una expresión de alto orden totalmente expandida, los autores aplican la regla de la cadena directamente a través de la definición recursiva. La derivada de una curva SIDER-$n$ se expresa como la derivada total de una operación SLERP aplicada a dos ramas SIDER-$(n-1)$ de orden inferior.
3. Derivada Total de SLERP: Se deriva un bloque de construcción fundamental para la derivada total de un SLERP con extremos móviles y un parámetro de interpolación móvil. Esta fórmula tiene en cuenta la variación del ángulo geodésico y de los extremos.
4. Especialización: La fórmula recursiva general se especializa para:
   • Nodos de Interpolación Interiores: Donde las dos ramas de orden inferior se encuentran en el mismo punto de datos, permitiendo una combinación simplificada de derivadas.
   • Extremos: Donde una rama se evalúa fuera del intervalo de interpolación estándar (extrapolación), introduciendo términos de mapa logarítmico.
   • Puntos Medios: Ubicaciones situadas a mitad de camino entre el muestreo de datos, que son críticas para algoritmos de reconstrucción de alto orden.

Contribuciones Clave

• Fórmulas de Derivada Recursivas: El artículo proporciona un algoritmo recursivo compacto y listo para su implementación para computar $D_n(t) = \frac{d}{dt}S_n(t)$ para cualquier orden $n$. Esto evita expresiones excesivamente complejas para curvas de alto orden.
• Casos Base Explícitos: El caso base SIDER2 se trata explícitamente, y SIDER3 y SIDER4 se utilizan para ilustrar el mecanismo recursivo y el manejo específico de extremos y nodos interiores.
• Preservación de la Tangencia: Se proporciona una prueba rigurosa que demuestra que el vector derivada derivado analíticamente es inherentemente tangente a la esfera en cada punto reconstruido (tanto en los nodos de muestreo como en los puntos medios). Esta propiedad ($S_n(t) \cdot D_n(t) = 0$) surge de la estructura algebraica de la derivada de SLERP y no requiere un post-procesamiento de proyección.
• Derivadas en Puntos Medios: El artículo deriva fórmulas específicas para derivadas en coordenadas de índice de medio entero ($\xi = j + 1/2$), las cuales son esenciales para evaluar reconstrucciones entre muestras de datos en esquemas de volumen finito y de tipo WENO.

Resultos

• Recursión General: La derivada $D_n$ se computa evaluando la derivada total de SLERP en las dos ramas de orden inferior $A(t)$ y $B(t)$, escalada por el factor de reparametrización $n/(n-1)$.
• Simplificación de Nodos: En los nodos interiores $t_j = j/n$, la derivada se simplifica a una suma ponderada de las derivadas de las dos ramas de orden inferior que convergen en $p_j$. Los pesos suman $n/(n-1)$, reflejando el escalamiento de la parametrización.
• Comportamiento en los Extremos: En los extremos, la derivada incluye términos que involucran el mapa logarítmico ($\text{Log}_p$) para dar cuenta de la dirección de los valores de la rama extrapolada.
• Consistencia Geométrica: Las visualizaciones numéricas confirman que los vectores derivada computados yacen exactamente en el espacio tangente de la esfera en los puntos interpolados, validando la prueba teórica de tangencia.

Significado y Alcance
El artículo sostiene que estas fórmulas extienden el marco original de SIDER/SENO al suministrar la información diferencial necesaria para reconstrucciones con valores en la esfera. La importancia radica en permitir el uso de SIDER en métodos de volumen finito de alto orden, y en métodos de tipo ENO/WENO y SENO para leyes de conservación y problemas de evolución donde el campo reconstruido debe diferenciarse mientras preserva estrictamente la restricción de la variedad. El trabajo se presenta como una teoría analítica que interpreta el cálculo de la derivada como una diferenciación de modo hacia adelante del algoritmo de evaluación de SIDER.

Los autores señalan que la implementación sigue siendo compacta y sigue directamente la definición recursiva. El artículo concluye reconociendo que, si bien la profundidad matemática de los resultados es significativa, el trabajo está destinado a ser un preprint para su difusión y futura verificación formal por pares, con posibles trabajos futuros que involucren espaciado de parámetros no uniformes y extensiones a otras variedades como los grupos de rotación.