Abstract:An algorithm for computing polynomial zeros, based on Aberth's method, is presented. The starting approximations are chosen by means of a suitable application of Rouché's theorem. More precisely, an integerq ≥ 1 and a set of annuliAi,i=1,...,q, in the complex plane, are determined together with the numberki of zeros of the polynomial contained in each annulusAi. As starting approximations we chooseki complex numbers lying on a suitable circle contained in the annulusAi, fori=1,...,q. The computation of Newton's correction is performed in such a way that overflow situations are removed. A suitable stop condition, based on a rigorous backward rounding error analysis, guarantees that the computed approximations are the exact zeros of a “nearby” polynomial. This implies the backward stability of our algorithm. We provide a Fortran 77 implementation of the algorithm which is robust against overflow and allows us to deal with polynomials of any degree, not necessarily monic, whose zeros and coefficients are representable as floating point numbers. In all the tests performed with more than 1000 polynomials having degrees from 10 up to 25,600 and randomly generated coefficients, the Fortran 77 implementation of our algorithm computed approximations to all the zeros within the relative precision allowed by the classical conditioning theorems with 11.1 average iterations. In the worst case the number of iterations needed has been at most 17. Comparisons with available public domain software and with the algorithm PA16AD of Harwell are performed and show the effectiveness of our approach. A multiprecision implementation in MATHEMATICA™ is presented together with the results of the numerical tests performed.

On computing zeros of a bivariate Bernstein polynomial

Approximate Degree Reduction of Interval Triangular Bézier Surfaces

Bernstein Type Theorems for Minimal Surfaces in $(\Alpha,\beta)$-Space

Bernstein Polynomial Composition Through Interpolation and Its Applications in Curves and Surfaces

Zeros of univariate interval polynomials

On the zeros of certain composite polynomials and an operator preserving inequalities

How to solve a binary cubic equation in integers

Parallel computation of real solving bivariate polynomial systems by zero-matching method

A Bernstein Polynomial Collocation Method for the Solution of Elliptic Boundary Value Problems

Fast Algorithm For Composition Of The Bernstein Polynomials

Bernstein-Sato polynomials for projective hypersurfaces with weighted homogeneous isolated singularities

Integral zeros of quadratic polynomials avoiding sublattices

Iterated integrals and Borwein-Chen-Dilcher polynomials

Distribution of the zeros of polynomials near the unit circle

Zero loci of Bernstein-Sato ideals

Bernstein--Sato polynomials of locally quasi-homogeneous divisors in $\mathbb{C}^{3}$

Counting solutions of a polynomial system locally and exactly

Zeros of the Selberg zeta function for symmetric infinite area hyperbolic surfaces

Special cubic zeros and the dual variety

A Fast Algorithm for Computing Macaulay Null Spaces of Bivariate Polynomial Systems

Numerical computation of polynomial zeros by means of Aberth's method