Parallel Multi-Step Contour Integral Methods for Nonlinear Eigenvalue Problems

We consider nonlinear eigenvalue problems to compute all eigenvalues in a bounded region on the complex plane. Based on domain decomposition and contour integrals, two robust and scalable parallel multi-step methods are proposed. The first method 1) uses the spectral indicator method to find eigenvalues and 2) calls a linear eigensolver to compute the associated eigenvectors. The second method 1) divides the region into subregions and uses the spectral indicator method to decide candidate regions that contain eigenvalues, 2) computes eigenvalues in each candidate subregion using Beyn's method; and 3) verifies each eigenvalue by substituting it back to the system and computes the smallest eigenvalue. Each step of the two methods is carried out in parallel. Both methods are robust, accurate, and does not require prior knowledge of the number and distribution of the eigenvalues in the region. Examples are presented to show the performance of the two methods.