A certified iterative method for isolated singular roots. (March 2023)

• Record Type: Journal Article
• Authors: Mantzaflaris, Angelos; Mourrain, Bernard; Szanto, Agnes
• Is Part Of: Journal of symbolic computation. Volume 115(2023)
• Page Start: 223
• Page End: 247
• Publication Date: 2023-03
• Subjects: Polynomial equations -- Isolated solution -- Multiple root -- Certification -- Newton iteration -- Deflation; Mathematics -- Data processing -- Periodicals; Numerical analysis -- Data processing -- Periodicals; Automatic programming (Computer science) -- Periodicals; Mathématiques -- Informatique -- Périodiques; Analyse numérique -- Informatique -- Périodiques; Programmation automatique -- Périodiques; Automatic programming (Computer science); Mathematics -- Data processing; Numerical analysis -- Data processing; Periodicals; Electronic journals; 510.285
• DOI: 10.1016/j.jsc.2022.08.006
• Languages: English
• ISSNs: 0747-7171

Abstract

Abstract: In this paper we provide a new method to certify that a nearby polynomial system has a singular isolated root and we compute its multiplicity structure. More precisely, given a polynomial system f = ( f 1, …, f N ) ∈ C [ x 1, …, x n ] N, we present a Newton iteration on an extended deflated system that locally converges, under regularity conditions, to a small deformation of f such that this deformed system has an exact singular root. The iteration simultaneously converges to the coordinates of the singular root and the coefficients of the so-called inverse system that describes the multiplicity structure at the root. We use α -theory test to certify the quadratic convergence, and to give bounds on the size of the deformation and on the approximation error. The approach relies on an analysis of the punctual Hilbert scheme, for which we provide a new description. We show in particular that some of its strata can be rationally parametrized and exploit these parametrizations in the certification. We show in numerical experimentation how the approximate inverse system can be computed as a starting point of the Newton iterations and the fast numerical convergence to the singular root with its multiplicity structure, certified by our criteria.