Computing the field of values and pseudospectra using the Lanczos method with continuation

Braconnier, Thierry and Higham, Nicholas J. (1996) Computing the field of values and pseudospectra using the Lanczos method with continuation. BIT Numerical Mathematics, 36:3. pp. 422-440. ISSN 1572-9125

[thumbnail of l5001475853137wq.pdf] PDF
l5001475853137wq.pdf
Restricted to Repository staff only

Download (989kB)

Abstract

The field of values and pseudospectra are useful tools for understanding the behaviour of various matrix processes. To compute these subsets of the complex plane it is necessary to estimate one or two eigenvalues of a large number of parametrized Hermitian matrices; these computations are prohibitively expensive for large, possibly sparse, matrices, if done by use of the QR algorithm. We describe an approach based on the Lanczos method with selective reorthogonalization and Chebyshev acceleration that, when combined with continuation and a shift and invert technique, enables efficient and reliable computation of the field of values and pseudospectra for large matrices. The idea of using the Lanczos method with continuation to compute pseudospectra is not new, but in experiments reported here our algorithm is faster and more accurate than existing algorithms of this type.

Item Type: Article
Uncontrolled Keywords: Field of values - pseudospectra - Lanczos method - Chebyshev acceleration - continuation
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: 28 Jun 2006
Last Modified: 20 Oct 2017 14:12
URI: https://eprints.maths.manchester.ac.uk/id/eprint/332

Actions (login required)

View Item View Item