  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.
