You are here: MIMS > EPrints
MIMS EPrints

2014.49: On the stability of computing polynomial roots via confederate linearizations

2014.49: Yuji Nakatsukasa and Vanni Noferini (2014) On the stability of computing polynomial roots via confederate linearizations.

Full text available as:

PDF - Requires a PDF viewer such as GSview, Xpdf or Adobe Acrobat Reader
607 Kb


A common way of computing the roots of a polynomial is to nd the eigenvalues of a linearization, such as the companion (when the polynomial is expressed in the monomial basis), colleague (Chebyshev basis) or comrade matrix (general orthogonal polynomial basis). For the monomial case, many studies exist on the stability of linearization-based rootnding algorithms. By contrast, little seems to be known for other polynomial bases. This paper studies the stability of algorithms that compute the roots via linearization in nonmonomial bases, and has three goals. First we prove its normwise stability when the polynomial is properly scaled and the QZ algorithm (as opposed to the more commonly used QR algorithm) is applied to a comrade pencil associated with a Jacobi orthogonal polynomial. Second, we extend a result by Arnold that leads to a rst-order expansion of the backward error when the eigenvalues are computed via QR, which shows that the method can be unstable. Finally, we focus on the special case of Chebyshev basis, in particular the Chebfun rootnder: we discuss its stability and describe an optional functionality, made available for improved stability, for computing the roots of a general continuous function f(x), implemented in the recently updated version 5. The main message is that to guarantee backward stability QZ applied to a properly scaled pencil is necessary

Item Type:MIMS Preprint
Uncontrolled Keywords:polynomial, roots, Chebyshev polynomial, companion matrix, colleague matrix, comrade matrix, QZ, QR, recurrence relation, perturbation, orthogonal polynomial
Subjects:MSC 2000 > 65 Numerical analysis
MIMS number:2014.49
Deposited By:Dr V Noferini
Deposited On:29 September 2014

Download Statistics: last 4 weeks
Repository Staff Only: edit this item