Recursive integral method with Cayley transformation

Abstract

Summary The recently developed RIM (recursive integral method) finds eigenvalues in a region of the complex plane. It computes an indicator to test if the region contains eigenvalues using an approximate spectral projection. If the indicator shows that the region contains eigenvalues, it is subdivided into smaller regions, which are then recursively tested and subdivided until any eigenvalues are isolated to a specified precision. We propose an enhancement to RIM that uses Cayley transformations and Arnoldi's method to greatly decrease the number of factorizations required to solve the linear systems that arise from a discretized contour integral to compute the indicator, which substantially reduces the computational cost. We demonstrate the efficacy of our method with numerical examples and compare with the implicitly restarted Arnoldi method, as implemented in eigs in MATLAB.

A new finite element approach for the Dirichlet eigenvalue problem

“…We use the linear Lagrange element on a series of uniformly refined meshes for discretization. The spectral indicator method [5,6,4] is employed to compute the smallest eigenvalue of (3.9). The results are shown in Table 1, which confirms the second order convergence.…”

“…The new approach has the following characteristics: 1) it provides a new finite element methodology which is different than the classic Babuška-Osborn theory; 2) it can be applied to a large class of nonlinear eigenvalue problems [4]; and 3) combined with the spectral indicator method [5,6,4], it can be parallelized to compute many eigenvalues effectively.…”

“…The linear Lagrange finite element is used for discretization and the convergence is proved using the abstract approximation results of Karma [7,8]. In Section 4, the spectral indicator method [5,6,4] is applied to compute the eigenvalues of the unit square.…”

An inverse medium problem using Stekloff eigenvalues and a Bayesian approach

“…Thus the numerical methods need to compute Stekloff eigenvalues of (2.2) in a region on C close to the origin. A recently developed spectral indicator method (SIM) is a good fit for this case [17,18]. Given reconstructed eigenvalues λ, a rectangular region containing these eigenvalues is chosen.…”

A Multilevel Spectral Indicator Method for Eigenvalues of Large Non-Hermitian Matrices

“…Recently, a family of eigensolvers, called the spectral indicator methods (SIMs), was proposed [6,7,14]. The idea of SIMs is different from the classical eigensolvers.…”