Fair Six-Point Subdivision on Euclidean and Riemannian Surfaces

Resumen Técnico: Subdivisión de Seis Puntos Justa en Superficies Euclídeas y Riemannianas

Planteamiento del Problema
Los esquemas de subdivisión interpolatorios generan curvas suaves a partir de polígonos de control mediante la inserción iterativa de nuevos vértices mientras se retienen los existentes. Aunque el esquema clásico de cuatro puntos de Deslauriers–Dubuc (específicamente el esquema de Dyn–Gregory–Levin o DGL con tensión $\omega=1/16$) es ampliamente utilizado, este produce curvas límite con solo una regularidad $C^1$. Esta baja regularidad suele resultar en perfiles de curvatura altamente oscilatorios, creando artefactos tangibles en procesos posteriores como el cálculo de normales de superficie y la generación de curvas de desplazamiento. Aunque la "justicia" (suavidad de la curvatura) a menudo se impone mediante post-procesamiento, este artículo investiga si la justicia está intrínsecamente integrada en reglas de subdivisión clásicas específicas. Además, el artículo aborda la extensión de estas propiedades de justicia a variedades riemannianas de curvatura constante (la esfera $S^2$ y el plano hiperbólico $H^2$), donde las reglas de subdivisión lineales estándar no se aplican directamente.

Metodología
El artículo emplea una combinación de cálculo variacional, computación simbólica y análisis numérico para analizar el esquema de seis puntos de Deslauriers–Dubuc y su análogo riemanniano.

1. Identidad Variacional Local: El autor establece una identidad local sobre una configuración de seis puntos. Fija los cuatro vecinos de semientero mediante interpolación de Lagrange de grado 5 a partir de los seis puntos de referencia. Luego, define un funcional de energía discreta que penaliza la variación de las segundas diferencias (una linealización de la variación de la curvatura). Demuestra que el único minimizador de esta energía es exactamente el punto medio de seis puntos de DD.
2. Extensión a la Variedad: El artículo utiliza el marco de elevación "log-exponencial" (Wallner y Dyn; Rahman et al.). El esquema mapea seis puntos vecinos al plano tangente en el punto medio geodésico del segmento, aplica la máscara lineal de seis puntos y mapea el resultado de vuelta a la variedad mediante el mapa exponencial.
3. Modelado de la Variación de la Curvatura: El autor deriva la ecuación completa de Euler–Lagrange para la energía de variación de la curvatura $\int (\kappa'_g)^2 ds$ en formas de espacio de curvatura constante $K$. Linealiza esta ecuación para obtener un modelo de dos puntos para perfiles de curvatura geodésica, produciendo soluciones trigonométricas en la esfera y soluciones hiperbólicas-trigonométricas en el plano hiperbólico.
4. Diagnósticos Numéricos: El artículo implementa los esquemas en dominios euclídeos, esféricos e hiperbólicos. Utiliza diagnósticos de razón de diferencia para estimar exponentes de Hölder críticos y calcula una energía de variación de curvatura discreta consistente para comparar la "justicia" de los esquemas de cuatro, seis y ocho puntos.

Contribuciones Clave
El artículo realiza cinco contribuciones específicas, distinguiendo entre resultados nuevos y relatos consolidados de la literatura previa:

1. Identidad Variacional Local (Teorema 1): El autor demuestra que, en una configuración de seis puntos con vecinos anclados por Lagrange, el punto medio de DD de seis puntos es el único minimizador de una energía de variación de segunda diferencia discreta. Esto proporciona una "lectura de justicia" estructural de los pesos de DD, vinculándolos con la minimización de la variación de la curvatura en el régimen de pequeñas deflexiones.
2. Análisis de Regularidad: El artículo ofrece un relato consolidado de la regularidad para los esquemas DD de cuatro, seis y ocho puntos, confirmando los exponentes de Hölder críticos de $\approx 2$ ($C^1$), $\approx 2.830$ ($C^2$) y $\approx 3.55$ ($C^3$) respectivamente. Valida esto mediante diagnósticos numéricos de razón de diferencia.
3. Ecuación de Euler–Lagrange en Formas de Espacio: El artículo deriva la ecuación completa de Euler–Lagrange para la energía de variación de la curvatura en superficies de curvatura constante $K$, reteniendo todos los términos de variación. Presenta el modelo linealizado y sus soluciones en forma cerrada, identificando la resonancia geométrica en la distancia antipodal en la esfera.
4. Implementación del Esquema en la Variedad: El artículo identifica explícitamente la elevación log-exponencial del punto medio geodésico como el análogo establecido de Deslauriers–Dubuc en el ámbito riemanniano (citando a Wallner [43]). Proporciona implementaciones autónomas para $S^2$ y el disco de Poincaré ($H^2$), un límite de proximidad elemental en la esfera y una explicación precisa de por qué el argumento de la norma ambiental falla en el hiperboloide.
5. Evaluación Comparativa y Diagnósticos: El artículo presenta una batería de pruebas comparativas entre los esquemas de cuatro, seis y ocho puntos. Valida las implementaciones en variedades contra la teoría conocida de convergencia y regularidad, y demuestra que el esquema de seis puntos equilibra mejor la justicia, la localidad y la robustez que sus contrapartes de cuatro y ocho puntos.

Resultados

• Justicia (Fairness): Los resultados experimentales muestran que el esquema de seis puntos produce perfiles de curvatura significativamente más justos que el esquema de cuatro puntos DGL. En polígonos cerrados, la energía de variación de curvatura discreta del esquema de seis puntos es de 5 a casi dos órdenes de magnitud inferior a la del esquema DGL.
• Regularidad: El esquema de seis puntos alcanza una regularidad $C^2$ (exponente de Hölder crítico $\approx 2.830$), un nivel de clase superior a la regularidad $C^1$ del DGL, sin el soporte más amplio y los potenciales artefactos de oscilación (ringing) del esquema de ocho puntos.
• Comportamiento en la Variedad: El esquema en la variedad preserva la regularidad euclídea completa ($C^2$) tanto en la esfera como en el plano hiperbólico. El análisis de proximidad confirma que el esquema satisface la condición de proximidad de segundo orden requerida para la convergencia.
• Resultado Nulo sobre Ángulos de Giro: El artículo documenta un "resultado nulo de advertencia": la medición de los ángulos de giro de media arista en círculos muestreados no logra detectar el término dependiente de la curvatura predicho por el modelo linealizado al orden $O(h^3)$. Esto se debe a que el defecto de cuerda-arco independiente de $K$ entra en el mismo orden, enmascarando la señal. Los efectos de la curvatura son, en cambio, visibles a través de las desviaciones de proximidad.

Significancia y Reivindicaciones
El artículo es modesto en sus reivindicaciones de novedad, atribuyendo explícitamente la construcción del esquema en la variedad, su convergencia y su regularidad a trabajos previos (Wallner, Dyn, Rahman, et al.). La significancia de este trabajo radica en:

• Proporcionar una interpretación variacional local exacta para el esquema DD de seis puntos, caracterizándolo como un minimizador de la variación de la curvatura discreta.
• Ofrecer implementaciones explícitas y autónomas y un límite de proximidad para el esquema en superficies de curvatura constante.
• Demostrar, mediante pruebas empíricas, que el esquema de seis puntos ofrece un equilibrio óptimo entre justicia (suavidad de la curvatura), localidad (tamaño de soporte) y robustez (manejo de datos no uniformes) en comparación con las alternativas de cuatro y ocho puntos.
• Clarificar la consistencia geométrica de los modelos de variación de la curvatura en formas de espacio, incluyendo la ubicación de las resonancias y las limitaciones de ciertas mediciones diagnósticas.

El artículo concluye que, si bien el esquema de seis puntos no es una construcción nueva, su identificación como una regla que minimiza la justicia y su rendimiento validado en variedades de Riemann proporcionan una base sólida para el diseño geométrico de alta calidad tanto en espacios euclídeos como curvos.