QR factoring to compute the GCD of univariate approximate polynomials

  • Authors:
  • R.M. Corless;S.M. Watt;Lihong Zhi

  • Affiliations:
  • Ontario Res. Centre for Comput. Algebra, Univ. of Western Ontario, London, Ont., Canada;-;-

  • Venue:
  • IEEE Transactions on Signal Processing
  • Year:
  • 2004

Quantified Score

Hi-index 35.69

Visualization

Abstract

We present a stable and practical algorithm that uses QR factors of the Sylvester matrix to compute the greatest common divisor (GCD) of univariate approximate polynomials over R[x] or C[x]. An approximate polynomial is a polynomial with coefficients that are not known with certainty. The algorithm of this paper improves over previously published algorithms by handling the case when common roots are near to or outside the unit circle, by splitting and reversal if necessary. The algorithm has been tested on thousands of examples, including pairs of polynomials of up to degree 1000, and is now distributed as the program QRGCD in the SNAP package of Maple 9.