Numer. Math. (2005) Numerische Digital Object Identifier (DOI) 10.1007/s00211-005-0595-4 Mathematik
Fast and stable QR eigenvalue algorithms for generalized companion matrices and secular equations
Dario A. Bini1, Luca Gemignani1, Victor Y. Pan2 1 Dipartimento di Matematica, Universit`a di Pisa, Via F. Buonarroti 2, 56127 Pisa, Italy; e-mail: {gemignan, [email protected]}dm.unipi.it. 2 Department of Mathematics and Computer Science, Lehman College of CUNY, Bronx, NY 10468, USA; e-mail: [email protected].
Received May 31, 2005 / Revised version received December 29, 2004 Published online xxx 2005, c Springer-Verlag 2005
Summary. We introduce a class Cn of n × n structured matrices which in- cludes three well-known classes of generalized companion matrices: tridi- agonal plus rank-one matrices (comrade matrices), diagonal plus rank-one matrices and arrowhead matrices. Relying on the structure properties of Cn, we show that if A ∈ Cn then A = RQ ∈ Cn, where A = QR is the QR decomposition of A. This allows one to implement the QR iteration for computing the eigenvalues and the eigenvectors of any A ∈ Cn with O(n) arithmetic operations per iteration and with O(n) memory storage. This iter- ation, applied to generalized companion matrices, provides new O(n2) flops algorithms for computing polynomial zeros and for solving the associated (ra- tional) secular equations. Numerical experiments confirm the effectiveness and the robustness of our approach. Mathematics Subject Classification (2000): 65F15, 65H17