Resumen Técnico: Computación de Soluciones Singulares de Sistemas Polinómicos

Planteamiento del Problema

El artículo aborda la computación numérica de soluciones singulares aisladas a sistemas de ecuaciones polinómicas $f(z) = 0$ dentro del marco de la Geometría Algebraica Numérica (NAG). Los enfoques estándar utilizan la continuación de homotopía, rastreando una curva de solución $z(t)$ definida por $h(z(t), t) = 0$ desde un sistema inicial conocido hacia el sistema objetivo a medida que $t \to 0$.

El desafío principal surge en el régimen del "endgame" (final de trayectoria), donde $t$ se aproxima a cero y la solución $z^* = z(0)$ es singular (es decir, el Jacobiano $f_z(z^*)$ tiene rango deficiente). En esta región, los correctores estándar basados en Newton fallan debido a la singularidad del Jacobiano, y los métodos existentes a menudo tienen dificultades con las tasas de convergencia o requieren derivadas de orden superior computacionalmente costosas. Los autores se centran en sistemas donde la solución es aislada y la homotopía satisface condiciones específicas de transversalidad y regularidad (Suposición 1), lo que garantiza la existencia de una expansión de serie de Puiseux para la curva de solución.

Metodología

El enfoque propuesto, denominado Arclength Endgame (Endgame de Longitud de Arco), combina técnicas de continuación de pseudo-longitud de arco con la estimación de los parámetros de la serie de Puiseux. El método se divide en dos fases: una fase de predictor y una fase de corrector.

1. Fase de Predictor

El objetivo es estimar los términos principales de la serie de Puiseux $z(t) = z^* + \sum a_j t^{k_j/c}$ para predecir la solución en $t=0$ (o un $t$ pequeño).

• Estimación de Exponentes: Los autores introducen un nuevo método, cRATIO, para estimar la razón del exponente fraccionario principal $k_1/c$. Este método utiliza la razón de las normas $\|t \dot{z}(t) - t_1 \dot{z}(t_1)\| / \|z(t) - z(t_1)\|$, la cual converge a $k_1/c$ cuando $t \to 0$. Comparan este método contra otros existentes como cSORT (ensayo y error) y cLOG (muestreo de secuencia geométrica), señalando la estabilidad empírica de cRATIO.
• Estimación de Orden Superior: Para sistemas que requieren precisión de orden superior, el artículo propone un método recursivo para estimar las razones subsecuentes $k_i/c$ (para $i \ge 2$) aplicando la lógica de estimación de $k_1/c$ a funciones transformadas derivadas de la serie original.
• Estimación de Coeficientes: Los autores presentan un nuevo esquema para estimar los coeficientes $a_{k_i}$ mediante la inversión de una matriz de Vandermonde escalada. A diferencia de la interpolación de Hermite clásica que asume exponentes densos $(1, 2, 3)$, este método admite exponentes dispersos y utiliza derivadas de orden superior (aproximadas mediante Diferencias Finitas de Hermite para evitar el cálculo de la Hessiana).

2. Fase de Corrector

Una vez estimado un punto predictor $\hat{x}$, el sistema se aumenta para restaurar la no-singularidad para el método de Newton.

• Corrank-1 ($\kappa=1$): El sistema $h(x)=0$ se aumenta con una única ecuación lineal $q \cdot (x - \hat{x}) = 0$, donde $q$ es un vector ortogonal a la tangente de la curva de solución (derivado del espacio nulo del Jacobiano aumentado $[h_z, h_t]$). Esto crea un sistema cuadrado bien condicionado, resoluble mediante el método de Newton.
• Corank $\kappa \ge 2$ (Lifted Arclength Endgame): Para corankes más altos, el método introduce $\kappa-1$ variables auxiliares y $\kappa$ hiperplanos para regularizar el sistema. Dado que el punto resultante puede no pertenecer a la homotopía original, los autores proponen una estrategia "elevada" (lifted): la homotopía se reinicia dinámicamente a una homotopía de Newton $h_v(z, t) = f(z) - t f(v)$, donde $v$ es la solución refinada. Esto permite que el proceso continúe refinando la solución cerca de la singularidad.

Contribuciones Clave

1. Arclength Endgame para Corank-1: Los autores introducen un nuevo algoritmo para sistemas de corank-1 que logra una convergencia superlineal en un entorno de la raíz singular. Crucialmente, este método depende únicamente de las evaluaciones del sistema y su Jacobiano, evitando el cálculo de derivadas de orden superior requeridas por los métodos de deflación con convergencia superlineal probada.
2. Lifted Arclength Endgame: Una extensión heurística para sistemas con corank $\kappa \ge 2$. Aunque no se proporcionan pruebas formales de convergencia para este caso, los resultados experimentales demuestran un comportamiento de convergencia superlineal prometedor.
3. Mejora en la Estimación de la Serie de Puiseux:
   • Un nuevo método estable (cRATIO) para estimar el exponente fraccionario principal $k_1/c$.
   • Un nuevo marco recursivo para estimar exponentes fraccionarios de orden superior ($k_i/c$ para $i \ge 2$) y sus coeficientes asociados, permitiendo predictores de orden superior para series de Puiseux dispersas.
4. Garantías Teóricas: El artículo proporciona pruebas formales (Teorema 1 y Corolario 2) que establecen la tasa de convergencia superlineal para el caso de corank-1, dependiendo de la precisión de la estimación del exponente y el entorno de la raíz.

Resultados

Los autores comparan sus métodos (AL para corank-1, LAL para corank $\ge 2$) contra el final de trayectoria de potencia clásica (CLASSIC) utilizando una suite de sistemas polinómicos, incluyendo los sistemas Griewank-Osborne, Lecerf y Caprasse, así como sistemas generados a medida con corankes y números de giro específicos.

• Convergencia: Los gráficos (Figuras 1–12) ilustran el residuo de seguimiento $\log_{10} \|f(z)\|$ frente al número acumulado de inversiones de matrices. Los métodos propuestos generalmente demuestran una convergencia más rápida (descenso más pronunciado en el residuo) y la capacidad de alcanzar residuos más profundos con menos inversiones de matrices en comparación con el enfoque clásico, particularmente en el régimen profundo del endgame.
• Estabilidad: La comparación de los métodos de estimación de exponentes (Figuras 13–18) sugiere que el método cRATIO propuesto ofrece una estabilidad competitiva o superior respecto a cSORT y cLOG en diversas instancias de prueba.
• Robustez: El enfoque elevado (LAL) maneja con éxito sistemas con corankes más altos donde los endgames estándar típicamente fallan o se estancan, manteniendo la convergencia mediante el reinicio dinámico de la homotopía.

Significación y Reivindicaciones

El artículo afirma que el Arclength Endgame ofrece un avance significativo en la computación de soluciones singulares al lograr una convergencia superlineal sin la sobrecarga computacional de las técnicas de deflación (que requieren derivadas de orden superior y pueden aumentar significativamente la dimensión del sistema).

Los autores enfatizan que su método es "consciente del endgame" (endgame-aware), utilizando explícitamente las propiedades geométricas de la curva de solución (serie de Puiseux) para regularizar el sistema. Posicionan este trabajo como una mejora práctica para la Geometría Algebraica Numérica, proporcionando un método que es tanto teóricamente fundamentado (para corank-1) como empíricamente efectivo (para corank superior), apoyándose en evaluaciones estándar de la función y su Jacobiano. El trabajo no pretende resolver todas las singularidades universalmente, sino que apunta al régimen específico de raíces singulares aisladas bajo supuestos de homotopía genéricos.