Computing roots of polynomials by quadratic clipping

  • Authors:
  • Michael Bartoň;Bert Jüttler

  • Affiliations:
  • Johann Radon Institute for Computational and Applied Mathematics, Austria;Johannes Kepler University Linz, Institute of Applied Geometry, Austria

  • Venue:
  • Computer Aided Geometric Design
  • Year:
  • 2007

Quantified Score

Hi-index 0.00

Visualization

Abstract

We present an algorithm which is able to compute all roots of a given univariate polynomial within a given interval. In each step, we use degree reduction to generate a strip bounded by two quadratic polynomials which encloses the graph of the polynomial within the interval of interest. The new interval(s) containing the root(s) is (are) obtained by intersecting this strip with the abscissa axis. In the case of single roots, the sequence of the lengths of the intervals converging towards the root has the convergence rate 3. For double roots, the convergence rate is still superlinear (32). We show that the new technique compares favorably with the classical technique of Bezier clipping.