Efficient isolation of polynomial's real roots

  • Authors:
  • Fabrice Rouillier;Paul Zimmermann

  • Affiliations:
  • LORIA/INRIA/LIP6, LIP6, BP 168, 8 rue du capitaine Scott, F-75015 Paris, France;LORIA/INRIA, 615 rue du Jardin Botanique, F-54602 Villers-lès-Nancy Cedex, France

  • Venue:
  • Journal of Computational and Applied Mathematics - Special issue: Proceedings of the international conference on linear algebra and arithmetic, Rabat, Morocco, 28-31 May 2001
  • Year:
  • 2004

Quantified Score

Hi-index 0.01

Visualization

Abstract

This paper revisits an algorithm isolating the real roots of a univariate polynomial using Descartes' rule of signs. It follows work of Vincent, Uspensky, Collins and Akritas, Johnson, Krandick.Our first contribution is a generic algorithm which enables one to describe all the known algorithms based on Descartes' rule of sign and the bisection strategy in a unified framework.Using that framework, a new algorithm is presented, which is optimal in terms of memory usage and as fast as both Collins and Akritas' algorithm and Krandick's variant, independently of the input polynomial. From this new algorithm, we derive an adaptive semi-numerical version, using multi-precision interval arithmetic.We finally show that these critical optimizations have important consequences since our new algorithm still works with huge polynomials, including orthogonal polynomials of degree 1000 and more, which were out of reach of previous methods.