Fortran 77 Subroutines for Computing the Eigenvalues of Hamiltonian Matrices II
Please always quote using this URN:urn:nbn:de:0296-matheon-1679
- This article describes Fortran 77 subroutines for computing eigenvalues and invariant subspaces of Hamiltonian and skew-Hamiltonian matrices. The implemented algorithms are based on orthogonal symplectic decompositions, implying numerical backward stability as well as symmetry preservation for the computed eigenvalues. These algorithms are supplemented with balancing and block algorithms, which can lead to considerable accuracy and performance improvements. As a by-product, an efficient implementation for computing symplectic QR decompositions is provided. We demonstrate the usefulness of the subroutines for several, practically relevant examples.
Author: | Peter Benner, Daniel Kressner |
---|---|
URN: | urn:nbn:de:0296-matheon-1679 |
Referee: | Peter Deuflhard |
Document Type: | Preprint, Research Center Matheon |
Language: | English |
Date of first Publication: | 2005/01/24 |
Release Date: | 2004/05/10 |
Preprint Number: | 204 |