Premium
Matlab code for sorting real Schur forms
Author(s) -
Brandts J. H.
Publication year - 2002
Publication title -
numerical linear algebra with applications
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.02
H-Index - 53
eISSN - 1099-1506
pISSN - 1070-5325
DOI - 10.1002/nla.274
Subject(s) - matlab , mathematics , eigenvalues and eigenvectors , computation , schur product theorem , schur's theorem , schur decomposition , algebra over a field , diagonal , matrix (chemical analysis) , linear algebra , sort , algorithm , computer science , pure mathematics , arithmetic , schur complement , geometry , classical orthogonal polynomials , physics , gegenbauer polynomials , materials science , quantum mechanics , composite material , orthogonal polynomials , operating system
In Matlab 6, there exists a command to generate a real Schur form, wheras another transforms a real Schur form into a complex one. There do not exist commands to prescribe the order in which the eigenvalues appear on the diagonal of the upper (quasi‐) triangular factor T . For the complex case, a routine is sketched in Golub and Van Loan (Matrix Computations (3rd edn). The John Hopkins University Press: Baltimore and London, 1996), that orders the diagonal of T according to their distance to a target value τ . In this technical note, we give a Matlab routine to sort real Schur forms in Matlab. It is based on a block‐swapping procedure by Bai and Demmel (Linear Algebra and Its Applications 1993; 186 : 73) We also describe how to compute a partial real Schur form (see Saad (Numerical methods for large eigenvalue problems. Manchester University Press: Manchester, 1992.)) in case the matrix A is very large. Sorting real Schur forms, both partially and completely, has important applications in the computation of real invariant subspaces. Copyright © 2002 by John Wiley & Sons, Ltd.