Digital Repository at the University of Maryland (DRUM)  >
College of Computer, Mathematical & Natural Sciences  >
Computer Science  >
Technical Reports from UMIACS 

Please use this identifier to cite or link to this item: http://hdl.handle.net/1903/12475

Title: Lyapunov Inverse Iteration for Computing a few Rightmost Eigenvalues of Large Generalized Eigenvalue Problems
Authors: Elman, Howard C.
Wu, Minghao
Type: Technical Report
Issue Date: 20-Apr-2012
Series/Report no.: UM Computer Science Department;CS-TR-5009
UMIACS;UMIACS-TR-2012-07
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. We also show that Lyapunov inverse iteration will always converge in only two steps if the Lyapunov equation in the first step is solved accurately enough. Furthermore, we generalize the analysis 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. Numerical experiments demonstrate the robustness of the algorithm.
URI: http://hdl.handle.net/1903/12475
Appears in Collections:Technical Reports of the Computer Science Department
Technical Reports from UMIACS

Files in This Item:

File Description SizeFormatNo. of Downloads
CS-TR-5009.pdf523.99 kBAdobe PDF132View/Open

All items in DRUM are protected by copyright, with all rights reserved.

 

DRUM is brought to you by the University of Maryland Libraries
University of Maryland, College Park, MD 20742-7011 (301)314-1328.
Please send us your comments