2013.49: NLEIGS: A class of robust fully rational Krylov methods for nonlinear eigenvalue problems
2013.49: Stefan Güttel, Roel Van Beeumen, Karl Meerbergen and Wim Michiels (2013) NLEIGS: A class of robust fully rational Krylov methods for nonlinear eigenvalue problems.
Full text available as:
|PDF - Requires a PDF viewer such as GSview, Xpdf or Adobe Acrobat Reader|
A new rational Krylov method for the efficient solution of nonlinear eigenvalue problems is proposed. This iterative method, called fully rational Krylov method for nonlinear eigenvalue problems (abbreviated as NLEIGS), is based on linear rational interpolation and generalizes the Newton rational Krylov method proposed in [R. Van Beeumen, K. Meerbergen, and W. Michiels, SIAM J. Sci. Comput., 35 (2013), pp. A327-A350]. NLEIGS utilizes a dynamically constructed rational interpolant of the nonlinear operator and a new companion-type linearization for obtaining a generalized eigenvalue problem with special structure. This structure is particularly suited for the rational Krylov method. A new approach for the computation of rational divided differences using matrix functions is presented. It is shown that NLEIGS has a computational cost comparable to the Newton rational Krylov method but converges more reliably, in particular, if the nonlinear operator has singularities nearby the target set. Moreover, NLEIGS implements an automatic scaling procedure which makes it work robustly independent of the location and shape of the target set, and it also features low-rank approximation techniques for increased computational efficiency. Small- and large-scale numerical examples are included.
|Item Type:||MIMS Preprint|
|Uncontrolled Keywords:||nonlinear eigensolver, rational Krylov, linear rational interpolation|
|Subjects:||MSC 2000 > 41 Approximations and expansions|
MSC 2000 > 47 Operator theory
MSC 2000 > 65 Numerical analysis
|Deposited By:||Stefan Güttel|
|Deposited On:||30 August 2013|