Abstract:This work is a continuation of "Fast and backward stable computation of roots of polynomials" by J.L. Aurentz, T. Mach, R. Vandebril, and D.S. Watkins, SIAM Journal on Matrix Analysis and Applications, 36(3): 942--973, 2015. In that paper we introduced a companion QR algorithm that finds the roots of a polynomial by computing the eigenvalues of the companion matrix in $O(n^{2})$ time using $O(n)$ memory. We proved that the method is backward stable. Here we introduce, as an alternative, a companion QZ algorithm that solves a generalized eigenvalue problem for a companion pencil. More importantly, we provide an improved backward error analysis that takes advantage of the special structure of the problem. The improvement is also due, in part, to an improvement in the accuracy (in both theory and practice) of the turnover operation, which is the key component of our algorithms. We prove that for the companion QR algorithm, the backward error on the polynomial coefficients varies linearly with the norm of the polynomial's vector of coefficients. Thus the companion QR algorithm has a smaller backward error than the unstructured QR algorithm (used by MATLAB's \texttt{roots} command, for example), for which the backward error on the polynomial coefficients grows quadratically with the norm of the coefficient vector. The companion QZ algorithm has the same favorable backward error as companion QR, provided that the polynomial coefficients are properly scaled.

A rectangular additive convolution for polynomials

Finite free convolutions via Weingarten calculus

The Newman algorithm for constructing polynomials with restricted coefficients and many real roots

Grothendieck inequalities characterize converses to the polynomial method

Expanding polynomials for sets with additive structure

Cumulants in rectangular finite free probability and beta-deformed singular values

A Unified Approach to Higher Order Convolutions Within a Certain Subset of Appell Polynomials

On converses to the polynomial method

Constant-Depth Arithmetic Circuits for Linear Algebra Problems

On the degree of polynomials computing square roots mod p

The fractional free convolution of R -diagonal elements and random polynomials under repeated differentiation

A better than $3/2$ exponent for iterated sums and products over $\mathbb R$

A polynomial with a root mod $p$ for every $p$ has a real root

On the Hadamard and Vandermonde Determinants and the Bernoulli–Euler–Lagrange–Aitken Method for Calculating the Roots of Polynomials

Generalized Convolution Properties Based on the Modified Mittag-Leffler Function

Convolution of periodic multiplicative functions and the divisor problem

Hankel determinants of convolution powers of Catalan numbers revisited

An asymptotic refinement of the Gauss-Lucas Theorem for random polynomials with i.i.d. roots

An analogue of the Rademacher function for generalized Dedekind sums in higher dimension

On the Log-Concavity of Hilbert Series of Veronese Subrings and Ehrhart Series

Fast and backward stable computation of roots of polynomials, Part II: backward error analysis; companion matrix and companion pencil