A Large, Sparse, and Indefinite Generalized Eigenvalue Problem from Fluid Mechanics

Mittelmann, Hans D.; Law, C. C.; Jankowski, D. F.; Neitzel, G. Paul

doi:10.1137/0913022

Cited by 11 publications

(1 citation statement)

References 14 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…For example, eigenvalues in the middle of the spectrum are needed for studying tidal motion [3], adaptive polynomial preconditioning [4,5], and Richardson's iteration [6], and the stability analysis in computational fluid dynamics needs estimates of interior eigenvalues [7]. Other applications include power system simulations [8], weather forecasting models, molecular chemistry, and so on.…”

Section: Introductionmentioning

confidence: 99%

A global harmonic Arnoldi method for large non-Hermitian eigenproblems with an application to multiple eigenvalue problems

Duan¹,

Jia²

2010

Journal of Computational and Applied Mathematics

View full text Add to dashboard Cite

a b s t r a c tThe global Arnoldi method can be used to compute exterior eigenpairs of a large nonHermitian matrix A, but it does not work well for interior eigenvalue problems. Based on the global Arnoldi process that generates an F -orthonormal basis of a matrix Krylov subspace, we propose a global harmonic Arnoldi method for computing certain harmonic F -Ritz pairs that are used to approximate some interior eigenpairs. We propose computing the F -Rayleigh quotients of the large non-Hermitian matrix with respect to harmonic F -Ritz vectors and taking them as new approximate eigenvalues. They are better and more reliable than the harmonic F -Ritz values. The global harmonic Arnoldi method inherits convergence properties of the harmonic Arnoldi method applied to a larger matrix whose distinct eigenvalues are the same as those of the original given matrix. Some properties of the harmonic F -Ritz vectors are presented. As an application, assuming that A is diagonalizable, we show that the global harmonic Arnoldi method is able to solve multiple eigenvalue problems both in theory and in practice. To be practical, we develop an implicitly restarted global harmonic Arnoldi algorithm with certain harmonic F -shifts suggested. In particular, this algorithm can be adaptively used to solve multiple eigenvalue problems. Numerical experiments show that the algorithm is efficient for the eigenproblem and is reliable for quite ill-conditioned multiple eigenproblems.

show abstract