You are here: MIMS > EPrints
MIMS EPrints

2007.224: Stability of Structured Hamiltonian Eigensolvers

2007.224: Françoise Tisseur (2001) Stability of Structured Hamiltonian Eigensolvers. SIAM Journal on Matrix Analysis and Applications, 23 (1). pp. 103-125. ISSN 1095-7162

Full text available as:

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

DOI: 10.1137/S0895479800368007


Various applications give rise to eigenvalue problems for which the matrices are Hamiltonian or skew-Hamiltonian and also symmetric or skew-symmetric. We define structured backward errors that are useful for testing the stability of numerical methods for the solution of these four classes of structured eigenproblems. We introduce the symplectic quasi-QR factorization and show that for three of the classes it enables the structured backward error to be efficiently computed. We also give a detailed rounding error analysis of some recently developed Jacobi-like algorithms of Fassbender, Mackey, and Mackey [Linear Algebra Appl., to appear] for these eigenproblems. Based on the direct solution of 4 × 4, and in one case 8 × 8, structured subproblems these algorithms produce a complete basis of symplectic orthogonal eigenvectors for the two symmetric cases and a symplectic orthogonal basis for all the real invariant subspaces for the two skew-symmetric cases. We prove that, when the rotations are implemented using suitable formulae, the algorithms are strongly backward stable and we show that the QR algorithm does not have this desirable property.

Item Type:Article
Uncontrolled Keywords:Hamiltonian; quaternion rotation; skew-Hamiltonian; symmetric; skew-symmetric; symplectic; backward error; structure-preserving; rounding error; Jacobi algorithm
Subjects:MSC 2000 > 15 Linear and multilinear algebra; matrix theory
MSC 2000 > 65 Numerical analysis
MIMS number:2007.224
Deposited By:Ms Helen Kirkbright
Deposited On:05 December 2007

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