Resumo Técnico: Computação de Soluções Singulares de Sistemas Polinomiais

Enunciado do Problema

O artigo aborda a computação numérica de soluções singulares isoladas para sistemas de equações polinomiais $f(z) = 0$ no âmbito da Geometria Algébrica Numérica (GAN). As abordagens padrão utilizam a continuação de homotopia, rastreando uma curva de solução $z(t)$ definida por $h(z(t), t) = 0$ de um sistema inicial conhecido para o sistema alvo conforme $t \to 0$.

O desafio principal surge no regime de "endgame" (final de jogo), onde $t$ se aproxima de zero e a solução $z^* = z(0)$ é singular (ou seja, o Jacobiano $f_z(z^*)$ é de posto deficiente). Nesta região, os corretores padrão baseados em Newton falham devido à singularidade do Jacobiano, e os métodos existentes frequentemente lutam com taxas de convergência ou exigem derivadas de ordem superior computacionalmente caras. Os autores focam em sistemas onde a solução é isolada e a homotopia satisfaz condições específicas de transversalidade e regularidade (Assunção 1), garantindo a existência de uma expansão de série de Puiseux para a curva de solução.

Metodologia

A abordagem proposta, denominada Arclength Endgame (Endgame de Comprimento de Arco), combina técnicas de continuação de pseudo-comprimento de arco com a estimativa de parâmetros de séries de Puiseux. O método é dividido em duas fases: uma fase de preditor e uma fase de corretor.

1. Fase de Preditor

O objetivo é estimar os termos principais da série de Puiseux $z(t) = z^* + \sum a_j t^{k_j/c}$ para prever a solução em $t=0$ (ou um $t$ pequeno).

• Estimativa de Expoentes: Os autores introduzem um novo método, cRATIO, para estimar a razão de expoente fracionário líder $k_1/c$. Este método utiliza a razão de normas $\|t \dot{z}(t) - t_1 \dot{z}(t_1)\| / \|z(t) - z(t_1)\|$, que converge para $k_1/c$ conforme $t \to 0$. Eles comparam este método contra métodos existentes como cSORT (tentativa e erro) e cLOG (amostragem de sequência geométrica), observando a estabilidade empírica do cRATIO.
• Estimativa de Ordem Superior: Para sistemas que exigem maior precisão, o artigo propõe um método recursivo para estimar as razões subsequentes $k_i/c$ (para $i \ge 2$) aplicando a lógica de estimativa de $k_1/c$ a funções transformadas derivadas da série original.
• Estimativa de Coeficientes: Os autores apresentam um novo esquema para estimar coeficientes $a_{k_i}$ usando a inversão de uma matriz de Vandermonde escalonada. Diferente da interpolação de Hermite clássica, que assume expoentes densos $(1, 2, 3)$, este método acomoda expoentes esparsos e utiliza derivadas de ordem superior (aproximadas via Diferenças Finitas de Hermite para evitar o cálculo da Hessiana).

2. Fase de Corretor

Uma vez estimado um ponto preditor $\hat{x}$, o sistema é aumentado para restaurar a não-singularidade para o método de Newton.

• Corank-1 ($\kappa=1$): O sistema $h(x)=0$ é aumentado com uma única equação linear $q \cdot (x - \hat{x}) = 0$, onde $q$ é um vetor ortogonal à tangente da curva de solução (derivado do núcleo do Jacobiano aumentado $[h_z, h_t]$). Isso cria um sistema quadrado bem condicionado, solúvel pelo método de Newton.
• Corank $\kappa \ge 2$ (Lifted Arclength Endgame): Para corank superior, o método introduz $\kappa-1$ variáveis auxiliares e $\kappa$ hiperplanos para regularizar o sistema. Como o ponto resultante pode não estar na homotopia original, os autores propõem uma estratégia de "elevação" (lifted): a homotopia é dinamicamente resetada para uma homotopia de Newton $h_v(z, t) = f(z) - t f(v)$, onde $v$ é a solução refinada. Isso permite que o processo continue refinando a solução próximo à singularidade.

Contribuições Principais

1. Arclength Endgame para Corank-1: Os autores introduzem um novo algoritmo para sistemas de corank-1 que alcança convergência superlinear em uma vizinhança da raiz singular. Crucialmente, este método depende apenas de avaliações do sistema e de seu Jacobiano, evitando o cálculo de derivadas de ordem superior exigidas por métodos de deflação com convergência superlinear comprovada.
2. Lifted Arclength Endgame: Uma extensão heurística para sistemas com corank $\kappa \ge 2$. Embora provas formais de convergência não sejam fornecidas para este caso, os resultados experimentais demonstram um comportamento promissor de convergência superlinear.
3. Melhoria na Estimativa de Séries de Puiseux:
   • Um novo método estável (cRATIO) para estimar o expoente fracionário líder $k_1/c$.
   • Um novo framework recursivo para estimar expoentes fracionários de ordem superior ($k_i/c$ para $i \ge 2$) e seus coeficientes associados, permitindo preditores de ordem superior para séries de Puiseux esparsas.
4. Garantias Teóricas: O artigo fornece provas formais (Teorema 1 e Corolário 2) estabelecendo a taxa de convergência superlinear para o caso de corank-1, dependente da precisão da estimativa do expoente e da vizinhança da raiz.

Resultados

Os autores testam seus métodos (AL para corank-1, LAL para corank $\ge 2$) contra o endgame clássico de série de potência (CLASSIC) utilizando uma suíte de sistemas polinomiais, incluindo os sistemas Griewank-Osborne, Lecerf e Caprasse, bem como sistemas gerados customizados com coranks e números de enrolamento (winding numbers) específicos.

• Convergência: Os gráficos (Figuras 1–12) ilustram o resíduo de rastreamento $\log_{10} \|f(z)\|$ versus o número cumulativo de inversões de matriz. Os métodos propostos geralmente demonstram convergência mais rápida (descida mais íngreme no resíduo) e capacidade de atingir resíduos mais profundos com menos inversões de matriz em comparação com a abordagem clássica, particularmente no regime profundo de endgame.
• Estabilidade: A comparação dos métodos de estimativa de expoente (Figuras 13–18) sugere que o método cRATIO proposto oferece estabilidade competitiva ou superior ao cSORT e cLOG em várias instâncias de teste.
• Robustez: A abordagem de elevação (LAL) lida com sucesso com sistemas de corank mais alto onde os endgames padrão tipicamente falham ou estagnam, mantendo a convergência através do reset dinâmico da homotopia.

Significância e Alegações

O artigo afirma que o Arclength Endgame oferece um avanço significativo no cálculo de soluções singulares ao alcançar convergência superlinear sem o custo computacional das técnicas de deflação (que exigem derivadas de ordem superior e podem aumentar significativamente a dimensão do sistema).

Os autores enfatizam que seu método é "consciente do endgame" (endgame-aware), utilizando explicitamente as propriedades geométricas da curva de solução (série de Puiseux) para regularizar o sistema. Eles posicionam o trabalho como uma melhoria prática para a Geometria Algébrica Numérica, fornecendo um método que é tanto teoricamente fundamentado (para corank-1) quanto empiricamente eficaz (para corank superior), enquanto depende apenas de avaliações padrão de função e Jacobiano. O trabalho não pretende resolver todas as singularidades universalmente, mas visa o regime específico de raízes singulares isoladas sob suposições de homotopia genéricas.