Algorithm 800 swMATH ID: 4405 Software Authors: Benner, Peter; Byers, Ralph; Barth, Eric Description: Algorithm 800. Fortran 77 subroutines for computing the eigenvalues of Hamiltonian matrices. I: The square-reduced method. The authors describe a set of LAPACK-based Fortran 77 subroutines for reducing a Hamiltonian matrix to square-reduced form using orthogonal symplectic transformations and for approximating all of its eigenvalues using an implicit version of Van Loan’s method. This approach requires less than half as many FLOPS as are needed by the general QR algorithm. Routines are provided for computing the square-reduced form, computing the eigenvalues using the output from the first routine, computing the eigenvalues (this routine is easy to use), and symplectic and norm scaling. There is a comprehensive discussion of the numerical tests used, which include several applications in control theory. Tests on randomly generated matrices confirm that the reduction in FLOPS is reflected in a similar reduction in execution times compared to the general QR routine, DGEEVX, from the LAPACK library. The authors’ software also preserves the plus-minus eigenvalue pairs of the Hamiltonian, which are destroyed by rounding errors using DGEEVX. The eigenvalues of larger magnitude are computed to an accuracy in line with their condition numbers; very small eigenvalues may be perturbed by up to the square root of the unit rounding error. Such a loss of accuracy may sometimes (but not always) be avoided by scaling. Of the three scaling strategies mentioned (norm, symplectic, and Hessenberg), norm scaling was found to be the least effective in practice. (review ACM) Homepage: http://dl.acm.org/citation.cfm?id=347837.347852 Programming Languages: Fortran 77 Keywords: algorithms; documentation; performance Related Software: LAPACK; SLICOT; Algorithm 854; CAREX; mctoolbox; ARPACK; EISPACK; Matlab; Yambo; mftoolbox; ELPA; BerkeleyGW; DMSRIC; DGRSVX; DSUBSP; HQR3; JDQZ; CTDSX; JDQR; BLAS Cited in: 11 Publications Standard Articles 1 Publication describing the Software, including 1 Publication in zbMATH Year Algorithm 800. Fortran 77 subroutines for computing the eigenvalues of Hamiltonian matrices. I: The square-reduced method. Zbl 1137.65338Benner, Peter; Byers, Ralph; Barth, Eric 2000 all top 5 Cited by 12 Authors 6 Benner, Peter 4 Kressner, Daniel 2 Mehrmann, Volker 2 Tisseur, Françoise 1 Barth, Eric 1 Burke, James V. 1 Byers, Ralph 1 Lewis, Adrian S. 1 Mackey, D. Steven 1 Mackey, Niloufer 1 Overton, Michael L. 1 Penke, Carolin all top 5 Cited in 6 Serials 2 ACM Transactions on Mathematical Software 2 Linear Algebra and its Applications 1 IMA Journal of Numerical Analysis 1 Journal of Computational and Applied Mathematics 1 SIAM Journal on Matrix Analysis and Applications 1 ELA. The Electronic Journal of Linear Algebra Cited in 5 Fields 11 Numerical analysis (65-XX) 2 History and biography (01-XX) 2 Linear and multilinear algebra; matrix theory (15-XX) 2 Computer science (68-XX) 1 Systems theory; control (93-XX) Citations by Year