Lyapunov inverse iteration for computing a few rightmost eigenvalues of large generalized eigenvalue problems

Howard C. Elman, Minghao Wu

Research output: Contribution to journalArticle

8 Scopus citations

Abstract

In linear stability analysis of a large-scale dynamical system, we need to compute the rightmost eigenvalue(s) for a series of large generalized eigenvalue problems. Existing iterative eigenvalue solvers are not robust when no estimate of the rightmost eigenvalue(s) is available. In this study, we show that such an estimate can be obtained from Lyapunov inverse iteration applied to a special eigenvalue problem of Lyapunov structure. An analysis that explains the fast convergence of this algorithm observed in numerical experiments is provided, based on which we propose a more efficient and robust algorithm. Furthermore, we generalize the same idea to a deflated version of this Lyapunov eigenvalue problem and propose an algorithm that computes a few rightmost eigenvalues for the eigenvalue problems arising from linear stability analysis.

Original languageEnglish (US)
Pages (from-to)1685-1707
Number of pages23
JournalSIAM Journal on Matrix Analysis and Applications
Volume34
Issue number4
DOIs
StatePublished - Dec 1 2013

    Fingerprint

Keywords

  • Inverse iteration
  • Linear stability analysis
  • Lyapunov solvers

ASJC Scopus subject areas

  • Analysis

Cite this