Local Consistency and Higher-Order Structure of Spherical Interpolation

Resumo Técnico: Consistência Local e Estrutura de Ordem Superior da Interpolação Esférica

Enunciado do Problema
O artigo aborda o desafio de interpolar dados restritos a uma variedade não linear, especificamente a esfera unitária $S^2$. A interpolação polinomial euclidiana padrão falha neste contexto porque a curva resultante geralmente deixa a variedade, e a projeção de volta altera tanto a geometria quanto as propriedades de aproximação. Embora a Interpolação Linear Esférica (SLERP) forneça um método de primeira ordem geometricamente natural que preserva a restrição de comprimento unitário, a SLERP por partes carece da suavidade necessária para uma reconstrução sistemática de ordem superior. Trabalhos anteriores introduziram o SIDER (Interpolação Esférica de ordem $n$), uma família recursiva de interpolantes inspirada na interpolação de Neville, mas que substitui combinações afins por SLERP. No entanto, uma análise rigorosa de consistência local provando a esperada precisão de ordem superior do SIDER para curvas suaves ainda não havia sido estabelecida.

Metodologia
A análise é conduzida utilizando coordenadas normais geodésicas centradas no ponto de avaliação exato $\gamma(s)$ em uma curva suave. Este sistema de coordenadas permite que a recursão SIDER seja comparada diretamente com a interpolação de Neville euclidiana clássica, enquanto retém explicitamente as correções dependentes da curvatura introduzidas pela SLERP.

A abordagem técnica central envolve:

1. Expansões de Taylor em Coordenadas Normais: Derivação de expansões locais para os pontos de dados amostrados $p_i = \gamma(ih)$ e para a curva exata no espaço tangente $T_{\gamma(s)}S^2$.
2. Expansão Cúbica da SLERP: Estabelecimento de uma expansão cúbica precisa para a operação SLERP em coordenadas normais. Uma descoberta crítica é que, embora a SLERP seja afim em primeira ordem, sua primeira correção não linear aparece na ordem cúbica. Isso é significativo porque o erro líder de uma reconstrução quadrática (SIDER2) também é de ordem $h^3$. Portanto, provas rigorosas para ordens superiores não podem tratar a parte não linear da SLERP como um resto negligenciável; a estrutura específica do termo cúbico deve ser analisada.
3. Expansão Formal Filtrada por Grau: Introdução de um arcabouço onde o erro local é tratado como uma série formal em $h$, com coeficientes que são polinômios na variável de stencil normalizada $\xi = \theta - i$. O artigo prova que a operação SLERP preserva uma propriedade específica de "filtração de grau": se os inputs possuem coeficientes de erro de grau no máximo $r$ em $\xi$ na ordem $h^r$, o output mantém essa estrutura.

Contribuições Principais e Resultados

• SIDER2 (Precisão de Terceira Ordem): O artigo deriva uma fórmula de erro refinada para o SIDER2. Prova que, apesar da geometria não linear, o termo de erro líder é exatamente $O(h^3)$ e possui a mesma estrutura nodal deslocada $\theta(\theta-1)(\theta-2)$ que a interpolação quadrática euclidiana ordinária. As correções de curvatura cúbica das etapas intermediárias da SLERP não alteram o coeficiente do erro líder.
• SIDER3 (Precisão de Quarta Ordem): A análise demonstra que duas reconstruções SIDER2 adjacentes que entram na recursão SIDER3 compartilham um coeficiente cúbico líder comum (o vetor $B$ que representa a terceira derivada da curva em coordenadas normais). Como os pesos na recursão SIDER3 coincidem com os do algoritmo de Neville, esses termos cúbicos líderes se cancelam exatamente. A mistura final por SLERP, aplicada a pontos que já estão a $O(h^3)$ de distância da solução, introduz erros não lineares de ordem $O(h^9)$, que são negligenciáveis. Isso resulta em uma precisão local de $O(h^4)$.
• SIDER4 (Precisão de Quinta Ordem): Ao estender a expansão para a ordem seguinte, o artigo verifica que os coeficientes de erro líderes para reconstruções SIDER3 adjacentes também são compatíveis (invariantes por deslocamento). Isso permite o cancelamento dos termos de ordem $O(h^4)$, provando a precisão de quinta ordem ($O(h^5)$) para o SIDER4.
• Generalização para Consistência de Todas as Ordens (Teorema 4): O artigo introduz um arcabouço de expansão filtrada por grau para generalizar esses resultados. Prova que, para qualquer $n$ fixo, a recursão SIDER-$n$ preserva a estrutura de grau polinomial dos coeficientes de erro. Combinado com as condições de interpolação em $n+1$ nós, isso força os primeiros $n$ coeficientes de erro a desaparecerem identicamente. Consequentemente, a estimativa de consistência local é estabelecida como:
  $$d_{S^2}(\gamma(\theta h), P^{[n]}_i(\theta; h)) = O(h^{n+1})$$
  sob suposições de suavidade ($C^{n+1}$) e uma condição de stencil pequeno (pontos permanecendo em uma vizinhança geodesicamente convexa comum).

Significância e Alegações
O artigo afirma fornecer uma explicação local rigorosa para o comportamento de alta ordem observado na construção SIDER. Sua principal significância reside em resolver a tensão entre o mecanismo de cancelamento algébrico da interpolação de Neville e a não linearidade geométrica da SLERP.

Os autores enfatizam que o sucesso da recursão SIDER não depende meramente do cancelamento algébrico de fatores nodais (como no espaço euclidiano), mas de uma compatibilidade de coeficientes mais profunda: os vetores de erro líderes de reconstruções de ordem inferior adjacentes devem ser idênticos quando expressos no mesmo sistema de coordenadas normais. O artigo prova essa compatibilidade explicitamente para SIDER2 através de SIDER4 e a estabelece genericamente para todas as ordens via o arcabouço de filtração por grau.

Os resultados são estritamente locais e assintóticos, assumindo que o stencil de interpolação e os valores intermediários permaneçam dentro de uma vizinhança geodesicamente convexa para garantir a unicidade das geodésicas e a validade do mapa logaritmo. O artigo não reivindica convergência global ou estabilidade para o método SENO (Esfericamente Essencialmente Não Oscilatório) associado, observando que uma teoria completa de não oscilação exigiria análise adicional do mecanismo de seleção de stencil.

Nota de Divulgação
O artigo divulga explicitamente que a modelagem computacional primária e as derivações matemáticas foram realizadas utilizando GPT-5.5 (Thinking). O autor atua como o investigador principal definindo o escopo, mas reconhece que a profundidade matemática dos resultados excede as capacidades de verificação manual atuais. Consequentemente, os achados são apresentados como hipóteses assistidas por IA destinadas à disseminação como um preprint, sujeitas a futura verificação formal por pares.