Stability of Structured Hamiltonian Eigensolvers

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

PDF tisseur4.pdf Download (227kB)

Abstract

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 2010, the AMS's Mathematics Subject Classification > 15 Linear and multilinear algebra; matrix theory MSC 2010, the AMS's Mathematics Subject Classification > 65 Numerical analysis
Depositing User:	Ms Lucy van Russelt
Date Deposited:	05 Dec 2007
Last Modified:	20 Oct 2017 14:12
URI:	https://eprints.maths.manchester.ac.uk/id/eprint/980

Actions (login required)

View Item