Fast and backward stable computation of the eigenvalues of matrix polynomials
Loading...
Date
2016-11-30
Authors
Aurentz, Jared
Mach, Thomas
Robol, Leonardo
Vandebril, Raf
Watkins, David S.
Journal Title
Journal ISSN
Volume Title
Publisher
arXiv
Abstract
In the last decade matrix polynomials have been investigated with the primary focus on adequate linearizations and good scaling techniques for computing their eigenvalues. In this article we propose a new backward stable method for computing a factored Schur form of the associated companion pencil. The algorithm has a quadratic cost in the degree of the polynomial and a cubic one in the size of the coefficient matrices. The algorithm is a variant of Francis's implicitly shifted QR algorithm applied on the associated companion pencil. A preprocessing unitary equivalence is executed on the matrix polynomial to simultaneously bring the leading matrix coefficient and the constant matrix term to triangular form before forming the companion pencil. The resulting structure allows us to stably factor both matrices of the pencil into $k$ matrices which are of unitary-plus-rank-one form admitting cheap and numerically reliable storage. The problem is then solved as a product core chasing eigenvalue problem. The numerical experiments illustrate stability and efficiency of the proposed methods.
Description
Keywords
math.NA
Citation
Aurentz, J., Mach, T., Robol, L., Vandebril, R., & Watkins, D. S. (2016). Fast and backward stable computation of the eigenvalues of matrix polynomials. arXiv, 1611(10142).