Geometry, Conditioning, and Parity in One-Sided Directional Curvature Recovery, with a Controlled Comparison of Estimators

Resumo Técnico: Geometria, Condicionamento e Paridade na Recuperação de Curvatura Direcional Unilateral

Enunciado do Problema
O artigo aborda os limites fundamentais da recuperação de uma matriz Hessiana (informação de derivada de segunda ordem) a partir de valores de função amostrados ao longo de raios que emanam de um único ponto. Este cenário surge na otimização sem derivadas (derivative-free optimization), aproximações de diferenças finitas próximas a fronteiras de domínio e problemas inversos onde o complemento de um domínio é inacessível. Nestes contextos, as direções admissíveis são restringidas pela geometria (ex: restrições de desigualdade ativas, estruturas finas ou fronteiras), impedindo a formação de estênciles simétricos. A questão central não é como projetar um estimador específico, mas sim o que é teoricamente possível dadas estas restrições geométricas: quais quantidades são identificáveis, como o condicionamento degrada com o acesso restrito e qual ordem de precisão é atingível.

Metodologia e Estrutura
O estudo adota uma abordagem geométrica e analítica em vez de propor um novo estimador numérico. A metodologia central envolve:

1. Modelo de Medição: Utilização de uma sonda direcional unilateral definida por três pontos em um raio: o centro $x$, um nó interior $x + R h u$ e um ponto final $x + h u$. A sonda forma uma segunda diferença dividida para estimar a curvatura direcional $u^T H u$.
2. Análise Geométrica: Análise da identificabilidade da Hessiana baseada na distribuição das direções amostradas. O autor utiliza o espaço de matrizes simétricas $3 \times 3$ ($\text{Sym}_3(\mathbb{R})$) e harmônicos esféricos para decompor o operador de medição.
3. Decomposição de Paridade: Exploração da paridade (simetria par/ímpar) de polinômios homogêneos na esfera. A expansão da sonda é dividida em partes pares (que carregam a Hessiana) e ímpares (que carregam o erro de truncamento principal).
4. Comparação Controlada: Um protocolo de benchmarking rigoroso comparando oito estimadores (incluindo regressão quadrática, mínimos quadrados móveis e vários esquemas direcionais) sob um "protocolo de justiça" (fairness protocol). Cada método recebe seu próprio raio de amostragem ótimo de oráculo, e todos compartilham as mesmas realizações de ruído e restrições de viabilidade.

Principais Contribuições e Resultados

1. Limites de Identificabilidade (Teorema 4, Proposição 5):

   • Piso de Avaliação: Em $\mathbb{R}^3$, sem priors de suavidade, pelo menos 10 avaliações de pontos em posição geral são necessárias e suficientes para determinar univocamente uma Hessiana.
   • Unicidade Direcional: Para medições direcionais, a condição para unicidade não é meramente que as direções gerem $\mathbb{R}^3$, mas que as direções não residam em um cone quádrico comum.
   • Implicação de Custo: Um esquema direcional usando $m$ direções requer $1 + 2m$ avaliações. Como $m \ge 6$ é necessário para a unicidade, o custo mínimo é de 13 avaliações, excedendo o piso de 10 avaliações para conjuntos de pontos gerais.

2. Condicionamento sob Acesso Restrito (Teorema 10):

   • Quando as direções admissíveis estão confinadas a um cap esférico de semi-ângulo $\theta$, os valores singulares do operador de recuperação escalam como $\{1, \theta, \theta, \theta^2, \theta^2, \theta^2\}$.
   • Consequentemente, o número de condição cresce como $\kappa(\theta) \asymp \theta^{-2}$. O artigo deriva constantes exatas para este crescimento (ex: $\kappa \theta^2 \to \sqrt{24} \approx 4.899$ para amostragem uniforme na superfície).
   • Interpretação: Cones estreitos medem bem a curvatura normal, mas mal as curvaturas mista e tangencial (ordens de $\theta$ e $\theta^2$ respectivamente). Uma fronteira suave (hemisfério, $\theta \approx 90^\circ$) é benigna; problemas de condicionamento surgem apenas em cantos, cúspides ou conjuntos ativos altamente restritos.

3. Decomposição de Paridade e Truncamento (Proposições 13, 14, Corolário 15):

   • A sonda unilateral expande-se em componentes de harmônicos esféricos pares (relacionados à Hessiana) e ímpares (relacionados ao truncamento).
   • Medição Livre de Hessiana: A projeção ímpar isola a derivada direcional de terceira ordem ($D^3 f$) independente da Hessiana, com precisão de $O(h^2)$. Isso permite a medição do termo de truncamento sem conhecimento prévio da Hessiana.
   • Ordem Atingível (Proposição 17): A precisão de segunda ordem na estimativa agregada da Hessiana é atingível se, e somente se, o conjunto de direções viáveis for antipodalmente simétrico ($u \in U \implies -u \in U$). Se o conjunto for genuinamente unilateral (sem pares antipodais), o erro permanece $O(h)$, independentemente do número de direções ou esquemas de ponderação.

4. Ruído e Design de Piloto (Proposições 19, 21, 22):

   • O artigo quantifica a variância da sonda e o consequente limite inferior para o passo utilizável (span). O span ótimo escala como $\sigma^{1/3}$, consistente com a teoria padrão de diferenças finitas.
   • A extração da constante de truncamento via projeção de paridade ímpar requer uma fase "piloto". O artigo observa que este piloto é caro, exigindo pelo menos 18 direções para ser estável, o que pode exceder o custo da própria estimativa da Hessiana.

5. Comparação Controlada de Estimadores:

   • Sob o protocolo de justiça, esquemas direcionais não superaram a regressão quadrática padrão ou a interpolação $\Lambda$-poised nos mesmos pontos viáveis.
   • Onde estênciles simétricos são viáveis, as diferenças centrais são mais precisas. Onde não são, a regressão quadrática e os mínimos quadrados móveis igualam ou superam os esquemas direcionais.
   • A falha dos esquemas direcionais em obter vantagem é atribuída ao fato de que eles não podem superar o erro de truncamento $O(h)$ inerente ao acesso unilateral (Proposição 17) e sofrem dos mesmos limites de condicionamento geomético que outros métodos.

Significância e Alegações
O artigo enquadra explicitamente suas contribuições como limites estruturais da geometria de medição, em vez de melhorias em algoritmos específicos.

• Restrições à Otimalidade: Os resultados definem a região do "impossível": nenhum estimador pode atingir precisão de segunda ordem sem simetria antipodal, e nenhum estimador pode contornar a penalidade de condicionamento $\theta^{-2}$ de cones estreitos.
• Estrutura Não Explorada: A identificação do componente de harmônico esférico ímpar como uma medida livre de Hessiana do termo de truncamento é destacada como um insight estrutural, embora o artigo alerte que extraí-lo é computacionalmente caro.
• Rigor Metodológico: Uma parte significativa do artigo é dedicada à documentação de nove artefatos experimentais específicos (ex: erros de propagação de parâmetros, bugs de geração de geometria e suposições de justiça falhas) que tiveram que ser removidos para produzir resultados confiáveis. O autor argumenta que o resultado negativo — que os esquemas direcionais não oferecem vantagem sobre a regressão sob condições justas — é mais valioso do que um potencial resultado positivo espúrio.

Modéstia e Limitações
O autor é cuidadoso ao distinguir entre teoremas provados e observações numéricas. Eles não alegam que a amostragem direcional unilateral seja superior perto de fronteiras; de fato, seus resultados e a comparação sugerem que não é. Eles reconhecem que seus resultados de decomposição de paridade podem ser instâncias da teoria de Slepian existente ou da literatura de imersão de Veronese restrita, enquadrando-os como "questões de posicionamento aberto" em vez de descobertas inéditas. O artigo conclui que, embora a geometria da recuperação unilateral esteja bem caracterizada, a utilidade prática dos esquemas direcionais é limitada pelos trade-offs fundamentais entre geometria, truncamento e ruído.