Deflation and certified isolation of singular zeros of polynomial systems

We develop a new symbolic-numeric algorithm for the certification of singular isolated points, using their associated local ring structure and certified numerical computations. An improvement of an existing method to compute inverse systems is presented, which avoids redundant computation and reduces the size of the intermediate linear systems to solve. We derive a one-step deflation technique, from the description of the multiplicity structure in terms of differentials. The deflated system can be used in Newton-based iterative schemes with quadratic convergence. Starting from a polynomial system and a sufficiently small neighborhood, we obtain a criterion for the existence and uniqueness of a singular root of a given multiplicity structure, applying a well-chosen symbolic perturbation. Standard verification methods, based e.g. on interval arithmetic and a fixed point theorem, are employed to certify that there exists a unique perturbed system with a singular root in the domain. Applications to topological degree computation and to the analysis of real branches of an implicit curve illustrate the method.

[1]  S. Smale,et al.  Complexity of Bézout’s theorem. I. Geometric aspects , 1993 .

[2]  Maria Grazia Marinari,et al.  Gröbner duality and multiplicities in polynomial system solving , 1995, ISSAC '95.

[3]  Bernard Mourrain,et al.  Topology and arrangement computation of semi-algebraic planar curves , 2008, Comput. Aided Geom. Des..

[4]  Lihong Zhi,et al.  Determining singular solutions of polynomial systems via symbolic-numeric reduction to geometric involutive forms , 2012, J. Symb. Comput..

[5]  Siegfried M. Rump,et al.  Verified error bounds for multiple roots of systems of nonlinear equations , 2010, Numerical Algorithms.

[6]  Michael Francis Atiyah,et al.  Introduction to commutative algebra , 1969 .

[7]  Marc Giusti,et al.  On Location and Approximation of Clusters of Zeros: Case of Embedding Dimension One , 2007, Found. Comput. Math..

[8]  Anton Leykin,et al.  Higher-Order Deflation for Polynomial Systems With Isolated Singular Solutions , 2006, math/0602031.

[9]  Zhonggang Zeng,et al.  Computing the multiplicity structure in solving polynomial systems , 2005, ISSAC.

[10]  Daniel Lazard,et al.  Resolution des Systemes d'Equations Algebriques , 1981, Theor. Comput. Sci..

[11]  ZHONGGANG ZENG The Closedness Subspace Method for Computing the Multiplicity Structure of a Polynomial System , 2008 .

[12]  Anton Leykin,et al.  Newton's method with deflation for isolated singularities of polynomial systems , 2006, Theor. Comput. Sci..

[13]  B. Mourrain Isolated points, duality and residues , 1997 .

[14]  Ágnes Szántó,et al.  Nearest multivariate system with given root multiplicities , 2009, J. Symb. Comput..

[15]  Bernard Mourrain,et al.  Continued fraction expansion of real roots of polynomial systems , 2009, SNC '09.

[16]  Lihong Zhi,et al.  Computing the multiplicity structure from geometric involutive form , 2008, ISSAC '08.

[17]  Grégoire Lecerf Quadratic Newton Iteration for Systems with Multiplicity , 2002, Found. Comput. Math..

[18]  Bernard Mourrain,et al.  Subdivision methods for solving polynomial equations , 2009, J. Symb. Comput..

[19]  Zbigniew Szafraniec Topological degree and quadratic forms , 1999 .

[20]  T. Ojika,et al.  Deflation algorithm for the multiple roots of a system of nonlinear equations , 1983 .

[21]  David Eisenbud,et al.  An Algebraic Formula for the Degree of a C ∞ Map Germ , 1977 .

[22]  Hans J. Stetter Analysis of zero clusters in multivariate polynomial systems , 1996, ISSAC '96.

[23]  F. S. Macaulay,et al.  The Algebraic Theory of Modular Systems , 1972 .