Abstract:Isolating real roots of a square-free polynomial in a given interval is a fundamental problem in computational algebra. Subdivision based algorithms are a standard approach to solve this problem. For instance, Sturm's method, or various algorithms based on the Descartes's rule of signs. For the benchmark problem of isolating all the real roots of a polynomial of degree n and root separation σ, the size of the subdivision tree of most of these algorithms is bounded by O(log⁡1/σ) (assume σ<1). Moreover, it is known that this is optimal for subdivision algorithms that perform uniform subdivision, i.e., the width of the interval decreases by some constant. Recently Sagraloff (2012) and Sagraloff–Mehlhorn (2016) have developed algorithms for real root isolation that combine subdivision with Newton iteration to reduce the size of the subdivision tree to O(n(log⁡(nlog⁡1/σ))).We describe a subroutine that reduces the size of the subdivision tree of any subdivision algorithm for real root isolation. The subdivision tree size of our algorithm using predicates based on either the Descartes's rule of signs or Sturm sequences is bounded by O(nlog⁡n), which is close to the optimal value of O(n). The corresponding bound for the algorithm EVAL, which uses certain interval-arithmetic based predicates, is O(n2log⁡n). Our analysis differs in two key aspects from earlier approaches. First, we use the general technique of continuous amortization from Burr–Krahmer–Yap (2009), and extend it to handle non-uniform subdivisions; second, we use the geometry of clusters of roots instead of root bounds. The latter aspect enables us to derive a bound on the subdivision tree that is independent of the root separation σ. The number of Newton iterations is bounded by O(nlog⁡log⁡(1/σ)).

Numerical Instability of Algebraic Rootfinding Methods

Extraneous Fixed Points of Euler Iteration and Corresponding Sullivan’s Basin

A Parallel Algorithm For Finding All The Roots Of A Polynomial

Near optimal subdivision algorithms for real root isolation

On the stability analysis of numerical schemes for solving non-linear polynomials arises in engineering problems

Root Radii and Subdivision for Polynomial Root-Finding

A novel Newton-Raphson style root finding algorithm

Condition Number Based Complexity Estimate for Solving Polynomial Systems

Detecting Potential Instabilities of Numerical Algorithms

Survey of Real Root Finding of Univariate Polynomial Equation in CAGD/CG

New practical advances in polynomial root clustering

Continued Fraction Expansion of Real Roots of Polynomial Systems

Algebraic approach to stability results for Erdős-Ko-Rado theorem

Polynomial Root Isolation by Means of Root Radii Approximation

Improved Algebraic Degeneracy Testing

Dynamical and numerical analysis of classical multiple roots finding methods applied for different multiplicities

Robust Numerical Algebraic Geometry

Algebraic Conditions for Stability in Runge-Kutta Methods and Their Certification via Semidefinite Programming

A Further Analysis of Backward Error in Polynomial Deflation

The Bit Complexity of Dynamic Algebraic Formulas and their Determinants

Rank-stability of polynomial equations