Newton's method with deflation for isolated singularities of polynomial systems

Anton Leykin, Jan Verschelde, Ailing Zhao

Published 2004-08-30, updated 2004-10-13Version 2

We present a modification of Newton's method to restore quadratic convergence for isolated singular solutions of polynomial systems. Our method is symbolic-numeric: we produce a new polynomial system which has the original multiple solution as a regular root. Using standard bases, a tool for the symbolic computation of multiplicities, we show that the number of deflation stages is bounded by the multiplicity of the isolated root. Our implementation performs well on a large class of applications.