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

Resumen Técnico: Geometría, Acondicionamiento y Paridad en la Recuperación de la Curvatura Direccional Unilateral

Planteamiento del Problema
El artículo aborda los límites fundamentales de la recuperación de una matriz Hessiana (información de la derivada de segundo orden) a partir de valores de función muestreados a lo largo de rayos que emanan de un único punto. Este escenario surge en la optimización sin derivadas, en aproximaciones de diferencias finitas cerca de los límites del dominio y en problemas inversos donde el complemento de un dominio es inaccesible. En estos entornos, las direcciones admisibles están restringidas por la geometría (por ejemplo, restricciones de desigualdad activas, estructuras delgadas o límites), lo que impide la formación de estenciles simétricos. La cuestión central no es cómo diseñar un estimador específico, sino qué es teóricamente posible dadas estas restricciones geométricas: qué cantidades son identificables, cómo se degrada el acondicionamiento con el acceso restringido y qué orden de precisión es alcanzable.

Metodología y Marco de Trabajo
El estudio adopta un enfoque geométrico y analítico en lugar de proponer un nuevo estimador numérico. La metodología central consiste en:

1. Modelo de Medición: El uso de una sonda direccional unilateral definida por tres puntos en un rayo: el centro $x$, un nodo interior $x + R h u$ y un punto final $x + h u$. La sonda forma una segunda diferencia dividida para estimar la curvatura direccional $u^T H u$.
2. Análisis Geométrico: Análisis de la identificabilidad de la Hessiana basada en la distribución de las direcciones muestreadas. El autor utiliza el espacio de matrices simétricas $3 \times 3$ ($\text{Sym}_3(\mathbb{R})$) y armónicos esféricos para descomponer el operador de medición.
3. Descomposición de Paridad: Explotación de la paridad (simetría par/impar) de los polinomios homogéneos en la esfera. La expansión de la sonda se divide en partes pares (que portan la Hessiana) e impares (que portan el error de truncamiento principal).
4. Comparación Controlada: Un protocolo de benchmarking riguroso que compara ocho estimadores (incluyendo regresión cuadrática, mínimos cuadrados móviles y varios esquemas direccionales) bajo un "protocolo de equidad". A cada método se le otorga su propio radio de muestreo óptimo de tipo oráculo, y todos comparten las mismas realizaciones de ruido y restricciones de factibilidad.

Contribuciones Clave y Resultados

1. Límites de Identificabilidad (Teorema 4, Proposición 5):

   • Piso de Evaluación: En $\mathbb{R}^3$, sin prior de suavidad, son necesarios y suficientes al menos 10 puntos de evaluación en posición general para determinar unánimemente una Hessiana.
   • Unicidad Direccional: Para mediciones direccionales, la condición para la unicidad no es simplemente que las direcciones abarquen $\mathbb{R}^3$, sino que las direcciones no yacen sobre un cono cuadrático común.
   • Implicación de Costo: Un esquema direccional que utiliza $m$ direcciones requiere $1 + 2m$ evaluaciones. Dado que se requiere $m \ge 6$ para la unicidad, el costo mínimo es de 13 evaluaciones, superando el piso de 10 evaluaciones para conjuntos de puntos generales.

2. Acondicionamiento bajo Acceso Restringido (Teorema 10):

   • Cuando las direcciones admisibles se confinan a un casquete esférico de semiángulo $\theta$, los valores singulares del operador de recuperación escalan como $\{1, \theta, \theta, \theta^2, \theta^2, \theta^2\}$.
   • Consecuentemente, el número de condición crece como $\kappa(\theta) \asymp \theta^{-2}$. El artículo deriva constantes exactas para este crecimiento (por ejemplo, $\kappa \theta^2 \to \sqrt{24} \approx 4.899$ para muestreo uniforme en la superficie).
   • Interpretación: Los conos estrechos miden bien la curvatura normal, pero pobremente las curvaturas mixtas y tangenciales (a órdenes de $\theta$ y $\theta^2$ respectivamente). Un límite suave (hemisferio, $\theta \approx 90^\circ$) es benigno; los problemas de acondicionamiento surgen solo en esquinas, cúspides o conjuntos activos altamente restringidos.

3. Descomposición de Paridad y Truncamiento (Proposiciones 13, 14, Corolario 15):

   • La sonda unilateral se expande en componentes de armónicos esféricos pares (relacionados con la Hessiana) e impares (relacionados con el truncamiento).
   • Medición Libre de la Hessiana: La proyección impar aísla la derivada direccional de tercer orden ($D^3 f$) independientemente de la Hessiana, con una precisión de $O(h^2)$. Esto permite la medición del término de truncamiento sin conocimiento previo de la Hessiana.
   • Orden Alcanzable (Proposición 17): La precisión de segundo orden en la estimación agregada de la Hessiana es alcanzable si y solo si el conjunto de direcciones factibles es antipodialmente simétrico ($u \in U \implies -u \in U$). Si el conjunto es genuinamente unilateral (sin pares antipodales), el error permanece en $O(h)$ independientemente del número de direcciones o de los esquemas de ponderación.

4. Ruido y Diseño de Piloto (Proposiciones 19, 21, 22):

   • El artículo cuantifica la varianza de la sonda y el resultante límite inferior del paso utilizable (extensión). La extensión óptima escala como $\sigma^{1/3}$, consistente con la teoría estándar de diferencias finitas.
   • La extracción de la constante de truncamiento mediante la proyección de paridad impar requiere una fase "piloto". El autor señala que este piloto es costoso, requiriendo al menos 18 direcciones para ser estable, lo que puede exceder el costo de la estimación de la Hessiana en sí misma.

5. Comparación Controlada de Estimadores:

   • Bajo el protocolo de equidad, los esquemas direccionales no superaron a la regresión cuadrática estándar o a la interpolación $\Lambda$-poised sobre los mismos puntos factibles.
   • Donde los estenciles simétricos son factibles, las diferencias centrales son más precisas. Donde no lo son, la regresión cuadrática y los mínimos cuadrados móviles igualan o superan a los esquemas direccionales.
   • El hecho de que los esquemas direccionales no logren una ventaja se atribuye a que no pueden superar el error de truncamiento de $O(h)$ inherente al acceso unilateral (Proposición 17) y sufren los mismos límites de acondicionamiento geométrico que otros métodos.

Significado y Reivindicaciones
El artículo enmarca explícitamente sus contribuciones como límites estructurales de la geometría de la medición, más que como mejoras a algoritmos específicos.

• Restricciones sobre la Optimidad: Los resultados definen la región "imposible": ningún estimador puede lograr una precisión de segundo orden sin simetría antipodal, y ningún estimador puede eludir la penalización de acondicionamiento $\theta^{-2}$ de los conos estrechos.
• Estructura No Explotada: Se destaca la identificación del componente de armónicos esféricos impares como una medida de la terminación libre de la Hessiana como un hallazgo estructural, aunque el autor advierte que extraerlo es computacionalmente costoso.
• Rigor Metodológico: Una parte significativa del documento se dedica a documentar nueve artefactos experimentales específicos (por ejemplo, errores de propagación de parámetros, errores en la generación de geometría y supuestos de equidad defectuosos) que tuvieron que ser eliminados para producir resultados fiables. El autor argumenta que el resultado negativo —que los esquemas direccionales no ofrecen ventaja sobre la regresión bajo condiciones de equidad— es más valioso que un posible resultado positivo espurio.

Modestia y Limitaciones
El autor es cuidadoso al distinguir entre teoremas probados y observaciones numéricas. No afirma que el muestreo direccional unilateral sea superior cerca de los límites; de hecho, sus resultados y la comparación sugieren que no lo es. Reconoce que sus resultados de descomposición de paridad pueden ser instancias de la teoría de Slepian existente o de la literatura de incrustación de Veronese restringida, enmarcándolos como "cuestiones de posicionamiento abierto" en lugar de descubrimientos novedosos. El artículo concluye que, si bien la geometría de la recuperación unilateral está bien caracterizada, la utilidad práctica de los esquemas direccionales está limitada por los compromisos fundamentales entre geometría, truncamiento y ruido.