Lyapunov Inverse Iteration for Computing a few Rightmost Eigenvalues of Large Generalized Eigenvalue Problems

dc.contributor.authorElman, Howard C.
dc.contributor.authorWu, Minghao
dc.date.accessioned2012-04-27T01:09:23Z
dc.date.available2012-04-27T01:09:23Z
dc.date.issued2012-04-20
dc.description.abstractIn 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.en_US
dc.identifier.urihttp://hdl.handle.net/1903/12475
dc.language.isoen_USen_US
dc.relation.ispartofseriesUM Computer Science Department;CS-TR-5009
dc.relation.ispartofseriesUMIACS;UMIACS-TR-2012-07
dc.titleLyapunov Inverse Iteration for Computing a few Rightmost Eigenvalues of Large Generalized Eigenvalue Problemsen_US
dc.typeTechnical Reporten_US

Files

Original bundle
Now showing 1 - 1 of 1
Loading...
Thumbnail Image
Name:
CS-TR-5009.pdf
Size:
523.99 KB
Format:
Adobe Portable Document Format
License bundle
Now showing 1 - 1 of 1
No Thumbnail Available
Name:
license.txt
Size:
1.57 KB
Format:
Item-specific license agreed upon to submission
Description: