Simple and Nearly Optimal Polynomial Root-finding by Means of Root Radii Approximation

We propose a new simple but nearly optimal algorithm for the approximation of all sufficiently well isolated complex roots and root clusters of a univariate polynomial. Quite typically the known root-finders at first compute some crude but reasonably good approximations to well-conditioned roots (that is, those isolated from the other roots) and then refine the approximations very fast, by using Boolean time which is nearly optimal, up to a polylogarithmic factor. By combining and extending some old root-finding techniques, the geometry of the complex plane, and randomized parametrization, we accelerate the initial stage of obtaining crude to all well-conditioned simple and multiple roots as well as isolated root clusters. Our algorithm performs this stage at a Boolean cost dominated by the nearly optimal cost of subsequent refinement of these approximations, which we can perform concurrently, with minimum processor communication and synchronization. Our techniques are quite simple and elementary; their power and application range may increase in their combination with the known efficient root-finding methods.

[1]  Alexandre Ostrowski Recherches sur la méthode de graeffe et les zéros des polynomes et des séries de laurent , 1940 .

[2]  A. Householder Dandelin, Lobacevskii, or Graeffe , 1959 .

[3]  A. Sluis Upperbounds for roots of polynomials , 1970 .

[4]  Ronald L. Graham,et al.  An Efficient Algorithm for Determining the Convex Hull of a Finite Planar Set , 1972, Inf. Process. Lett..

[5]  Henry C. Thacher,et al.  Applied and Computational Complex Analysis. , 1988 .

[6]  A. D. Briuno,et al.  Local methods in nonlinear differential equations , 1989 .

[7]  Victor Y. Pan,et al.  Optimal (up to polylog factors) sequential and parallel algorithms for approximating complex polynomial zeros , 1995, STOC '95.

[8]  Joachim von zur Gathen,et al.  Fast algorithms for Taylor shifts and certain difference equations , 1997, ISSAC.

[9]  Peter Kirrinnis,et al.  Partial Fraction Decomposition in (z) and Simultaneous Newton Iteration for Factorization in C[z] , 1998, J. Complex..

[10]  Victor Y. Pan,et al.  Approximating Complex Polynomial Zeros: Modified Weyl's Quadtree Construction and Improved Newton's Iteration , 2000, J. Complex..

[11]  V. Pan Structured Matrices and Polynomials , 2001 .

[12]  Gregorio Malajovich,et al.  Tangent Graeffe iteration , 2001, Numerische Mathematik.

[13]  V. Pan Structured Matrices and Polynomials: Unified Superfast Algorithms , 2001 .

[14]  Victor Y. Pan,et al.  Numerical methods for roots of polynomials , 2007 .

[15]  Victor Y. Pan,et al.  On the boolean complexity of real root refinement , 2013, ISSAC '13.

[16]  Victor Y. Pan,et al.  Accelerated approximation of the complex roots of a univariate polynomial , 2014, SNC.

[17]  Victor Y. Pan,et al.  Polynomial Real Root Isolation by Means of Root Radii Approximation , 2015, CASC.

[18]  Victor Y. Pan,et al.  Nearly optimal refinement of real roots of a univariate polynomial , 2016, J. Symb. Comput..

[19]  Victor Y. Pan,et al.  Accelerated approximation of the complex roots and factors of a univariate polynomial , 2015, Theor. Comput. Sci..