Analytic First Derivatives of SIDER Interpolation

Resumo Técnico: Derivadas Analíticas de Primeira Ordem de Interpolação SIDER

Enunciado do Problema
Na análise numérica e modelagem geométrica, há uma necessidade crescente de interpolar dados restritos a variedades não lineares, especificamente a esfera unitária $S^2 \subset \mathbb{R}^3$. A interpolação polinomial euclidiana padrão falha em preservar a restrição esférica entre os pontos de dados, muitas vezes exigindo normalização pontual, o que complica a precisão e a análise de derivadas. Embora a Interpolação Esférica de ordem $n$ (SIDER-$n$) tenha sido introduzida anteriormente como uma ferramenta de reconstrução de alta ordem para métodos de Esfericamente Essencialmente Não Oscilatórios (SENO), ela carecia de uma teoria analítica para derivadas de primeira ordem. A ausência de informações de derivada limita a aplicação do SIDER em esquemas numéricos que requerem inclinações locais, mudanças direcionais ou quantidades diferenciais (por exemplo, resolvedores de Riemann generalizados, procedimentos de Lax–Wendroff inverso e quadratura de alta ordem).

Metodologia
O artigo desenvolve fórmulas analíticas de derivadas de primeira ordem para curvas SIDER de ordem arbitrária, explorando a estrutura recursiva da interpolação. A metodologia central baseia-se nos seguintes princípios:

1. Composição Recursiva: O SIDER-$n$ é construído inteiramente a partir de operações de Interpolação Linear Esférica (SLERP), formando uma árvore binária de mapas SLERP.
2. Propagação da Regra da Cadeia: Em vez de diferenciar uma expressão de alta ordem totalmente expandida, os autores aplicam a regra da cadeia diretamente através da definição recursiva. A derivada de uma curva SIDER-$n$ é expressa como a derivada total de uma operação SLERP aplicada a dois ramos SIDER-$(n-1)$ de ordem inferior.
3. Derivada Total do SLERP: Um bloco fundamental é derivado para a derivada total do SLERP com extremidades móveis e um parâmetro de interpolação móvel. Esta fórmula contabiliza a variação do ângulo geodésico e das extremidades.
4. Especialização: A fórmula recursiva geral é especializada para:
   • Nós de Interpolação Interiores: Onde os dois ramos de ordem inferior se encontram no mesmo ponto de dado, permitindo uma combinação ponderada simplificada de derivadas.
   • Extremidades: Onde um ramo é avaliado fora do intervalo de interpolação padrão (extrapolação), introduzindo termos de mapa logarítmico.
   • Pontos Médios: Localizações situadas no meio do caminho entre a amostragem de dados, que são críticas para algoritmos de reconstrução baseados em ordem superior.

Principais Contribuições

• Fórmulas de Derivada Recursivas: O artigo fornece um algoritmo recursivo compacto e pronto para implementação para computar $D_n(t) = \frac{d}{dt}S_n(t)$ para qualquer ordem $n$. Isso evita expressões excessivamente complexas para curvas de alta ordem.
• Casos Base Explícitos: O caso base SIDER2 é tratado explicitamente, e o SIDER3 e SIDER4 são usados para ilustrar o mecanismo recursivo e o tratamento específico de extremidades e nós interiores.
• Preservação de Tangência: Uma prova rigorosa é fornecida demonstrando que o vetor derivada derivado analiticamente é inerentemente tangente à esfera em cada ponto reconstruído (tanto nos nós de amostragem quanto nos pontos médios). Esta propriedade ($S_n(t) \cdot D_n(t) = 0$) surge da estrutura algébrica da derivada do SLERP e não requer pós-processamento de projeção.
• Derivadas de Pontos Médios: O artigo deriva fórmulas específicas para derivadas em coordenadas de índice de meio-inteiro ($\xi = j + 1/2$), que são essenciais para avaliar reconstruções entre amostras de dados em esquemas de volume finito e do tipo WENO.

Resultos

• Recursão Geral: A derivada $D_n$ é computada avaliando a derivada total do SLERP nos dois ramos de ordem inferior $A(t)$ e $B(t)$, escalonados pelo fator de reparametrização $n/(n-1)$.
• Simplificação de Nós: Nos nós interiores $t_j = j/n$, a derivada simplifica-se para uma soma ponderada das derivadas dos dois ramos de ordem inferior que se encontram em $p_j$. Os pesos somam $n/(n-1)$, refletindo o escalonamento da parametrização.
• Comportamento nas Extremidades: Nas extremidades, a derivada inclui termos envolvendo o mapa logarítmico ($\text{Log}_p$) para contabilizar a direção dos valores do ramo extrapolado.
• Consistência Geométrica: Visualizações numéricas confirmam que os vetores derivada computados residem exatamente no espaço tangente da esfera nos pontos interpolados, validando a prova teórica de tangência.

Significância e Escopo
O artigo afirma que estas fórmulas estendem o framework original SIDER/SENO ao fornecer as informações diferenciais necessárias para reconstruções com valores na esfera. A significância reside em permitir o uso do SIDER em métodos de volume finito de alta ordem, ENO/WENO e do tipo SENO para leis de conservação e problemas de evolução onde o campo reconstruído deve ser diferenciado enquanto preserva estritamente a restrição da variedade. O trabalho é apresentado como uma teoria analítica que interpreta o cálculo da derivada como uma diferenciação de modo direto (forward-mode differentiation) do algoritmo de avaliação SIDER.

Os autores observam que a implementação permanece compacta e segue diretamente a definição recursiva. O artigo conclui reconhecendo que, embora a profundidade matemática dos resultados seja significativa, o trabalho é destinado como um preprint para disseminação e futura verificação formal por pares, com potencial para trabalhos futuros envolvendo espaçamento de parâmetros não uniformes e extensões para outras variedades, como grupos de rotação.