Spiral endgames for computing singular endpoints of homotopy paths

When the solution is singular, prediction-correction path tracking fails at the endpoint. Instead of tracking the path all the way to its end, endgames sample nonsingular points on the incoming path and approximate the endpoint by interpolating the samples with a truncated Puiseux series.

By collecting sample points along a logarithmic spiral in the complex plane, the proposed spiral endgame balances between the superior numerical conditioning of a circular sample (known as the Cauchy endgame) and the update efficiency of a linear sample (known as the power series endgame), these being tied together as two extremes of the new approach. The numerical conditioning of the methods are analyzed and compared, and computational examples illustrate the effectiveness of the new spiral sampling method compared with the circular sampling method of the Cauchy endgame.</p>