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

Résumé Technique : Géométrie, Conditionnement et Parité dans la Récupération de Courbure Directionnelle Unilatérale

Énoncé du Problème
L'article traite des limites fondamentales de la récupération d'une matrice hessienne (information de dérivée du second ordre) à partir de valeurs de fonctions échantillonnées le long de rayons émanant d'un point unique. Ce scénario survient dans l'optimisation sans dérivée, les approximations de différences finies à proximité des frontières de domaine, et les problèmes inverses où le complément d'un domaine est inaccessible. Dans ces contextes, les directions admissibles sont contraintes par la géométrie (par exemple, des contraintes d'inégalité actives, des structures fines ou des frontières), empêchant la formation de stencils symétriques. La question centrale n'est pas de savoir comment concevoir un estimateur spécifique, mais plutôt ce qui est théoriquement possible compte tenu de ces contraintes géométriques : quelles quantités sont identifiables, comment le conditionnement se dégrade avec un accès restreint, et quel ordre de précision est atteignable.

Méthodologie et Cadre
L'étude adopte une approche géométrique et analytique plutôt que de proposer un nouvel estimateur numérique. La méthodologie centrale comprend :

1. Modèle de Mesure : Utilisation d'une sonde directionnelle unilatérale définie par trois points sur un rayon : le centre $x$, un nœud intérieur $x + R h u$, et un point terminal $x + h u$. La sonde forme une seconde différence divisée pour estimer la courbure directionnelle $u^T H u$.
2. Analyse Géométrique : Analyse de l'identifiabilité de la Hessienne basée sur la distribution des directions échantillonnées. L'auteur utilise l'espace des matrices symétriques $3 \times 3$ ($\text{Sym}_3(\mathbb{R})$) et les harmoniques sphériques pour décomposer l'opérateur de mesure.
3. Décomposition de Parité : Exploitation de la parité (symétrie paire/impaire) des polynômes homogènes sur la sphère. L'expansion de la sonde est divisée en parties paires (portant la Hessienne) et impaires (portant l'erreur de troncature principale).
4. Comparaison Contrôlée : Un protocole de benchmarking rigoureux comparant huit estimateurs (incluant la régression quadratique, les moindres carrés mobiles et diverses schémas directionnels) sous un « protocole d'équité ». Chaque méthode se voit accorder son propre rayon d'échantillonnage oracle-optimal, et toutes partagent des réalisations de bruit et des contraintes de faisabilité identiques.

Contributions Clés et Résultats

1. Limites d'Identifiabilité (Théorème 4, Proposition 5) :

   • Plancher d'Évaluation : Dans $\mathbb{R}^3$, sans a priori de lissité, au moins 10 évaluations de points en position générale sont nécessaires et suffisantes pour déterminer une Hessienne de manière unique.
   • Unicité Directionnelle : Pour les mesures directionnelles, la condition d'unicité n'est pas simplement que les directions engendrent $\mathbb{R}^3$, mais que les directions ne reposent pas sur un cône quadrique commun.
   • Implication de Coût : Un schéma directionnel utilisant $m$ directions nécessite $1 + 2m$ évaluations. Puisque $m \ge 6$ est requis pour l'unicité, le coût minimum est de 13 évaluations, dépassant le plancher de 10 évaluations pour des ensembles de points généraux.

2. Conditionnement sous Accès Restreint (Théorème 10) :

   • Lorsque les directions admissibles sont confinées à un calotte sphérique de demi-angle $\theta$, les valeurs singulières de l'opérateur de récupération évoluent comme $\{1, \theta, \theta, \theta^2, \theta^2, \theta^2\}$.
   • Par conséquent, le nombre de conditionnement croît comme $\kappa(\theta) \asymp \theta^{-2}$. L'article dérive des constantes exactes pour cette croissance (par exemple, $\kappa \theta^2 \to \sqrt{24} \approx 4.899$ pour un échantillonnage uniforme sur la surface).
   • Interprétation : Les cônes étroits mesurent bien la courbure normale mais peu les courbures mixtes et tangentielles (aux ordres $\theta$ et $\theta^2$ respectivement). Une frontière lisse (hémisphère, $\theta \approx 90^\circ$) est bénigne ; les problèmes de conditionnement n'apparaissent qu'aux coins, aux cuspides ou aux ensembles actifs hautement contraints.

3. Décomposition de Parité et Troncature (Propositions 13, 14, Corollaire 15) :

   • La sonde unilatérale se développe en composantes d'harmoniques sphériques paires (liées à la Hessienne) et impaires (liées à la troncature).
   • Mesure Sans Hessienne : La projection impaire isole la dérivée directionnelle du troisième ordre ($D^3 f$) indépendamment de la Hessienne, avec une précision de $O(h^2)$. Cela permet de mesurer le terme de troncature sans connaissance préalable de la Hessienne.
   • Ordre Atteignable (Proposition 17) : Une précision de second ordre dans l'estimation agrégée de la Hessienne est atteignable si et seulement si l'ensemble des directions faisables est antipodalement symétrique ($u \in U \implies -u \in U$). Si l'ensemble est véritablement unilatéral (pas de paires antipodales), l'erreur reste de $O(h)$ quel que soit le nombre de directions ou les schémas de pondération.

4. Bruit et Conception de Pilote (Propositions 19, 21, 22) :

   • L'article quantifie la variance de la sonde et la limite inférieure résultante sur la taille de pas utilisable (étendue). L'étendue optimale évolue en $\sigma^{1/3}$, cohérent avec la théorie standard des différences finies.
   • L'extraction de la constante de troncature via la projection de parité impaire nécessite une phase « pilote ». L'auteur note que ce pilote est coûteux, nécessitant au moins 18 directions pour être stable, ce qui peut excéder le coût de l'estimation de la Hessienne elle-même.

5. Comparaison Contrôlée des Estimateurs :

   • Sous le protocole d'équité, les schémas directionnels n'ont pas surpassé la régression quadratique standard ou l'interpolation $\Lambda$-poisée sur les mêmes points faisables.
   • Là où des stencils symétriques sont réalisables, les différences centrales sont les plus précises. Là où ils ne le sont pas, la régression quadratique et les moindres carrés mobiles égalent ou dépassent les schémas directionnels.
   • L'échec des schémas directionnels à obtenir un avantage est attribué au fait qu'ils ne peuvent surmonter l'erreur de troncature de $O(h)$ inhérente à l'accès unilatéral (Proposition 17) et souffrent des mêmes limites de conditionnement géométrique que les autres méthodes.

Signification et Revendications
L'article cadre explicitement ses contributions comme des limites structurelles sur la géométrie de la mesure plutôt que comme des améliorations de algorithmes spécifiques.

• Contraintes sur l'Optimalité : Les résultats définissent la région de l'« impossible » : aucun estimateur ne peut atteindre une précision de second ordre sans symétrie antipodale, et aucun estimateur ne peut contourner la pénalité de conditionnement $\theta^{-2}$ des cônes étroits.
• Structure Inexploitée : L'identification de la composante d'harmonie sphérique impaire comme une mesure de la troncature sans la Hessienne est soulignée comme une intuition structurelle, bien que l'auteur mette en garde contre le fait que l'extraction de celle-ci est coûteuse en calcul.
• Rigueur Méthodologique : Une partie importante de l'article est consacrée à la documentation de neuf artefacts expérimentaux spécifiques (ex: erreurs de propagation de paramètres, bugs de génération de géométrie et hypothèses d'équité erronées) qui ont dû être éliminés pour produire des résultats dignes de confiance. L'auteur soutient que le résultat négatif — que les schémas directionnels n'offrent aucun avantage sur la régression sous des conditions équitables — est plus précieux qu'un résultat positif potentiellement spécieux.

Modestie et Limites
L'auteur prend soin de distinguer les théorèmes prouvés des observations numériques. Il ne prétend pas que l'échantillonnage directionnel unilatéral est supérieur près des frontières ; en fait, ses résultats et la comparaison suggèrent que ce n'est pas le cas. Il reconnaît que ses résultats de décomposition de parité peuvent être des instances de la théorie de Slepian existante ou de la littérature sur l'encastrement de Veronese restreint, les présentant comme des questions de positionnement ouvert plutôt que comme des découvertes nouvelles. L'article conclut que bien que la géométrie de la récupération unilatérale soit bien caractérisée, l'utilité pratique des schéments directionnels est limitée par les compromis fondamentaux entre géométrie, troncature et bruit.