English

Avoiding discretization issues for nonlinear eigenvalue problems

Numerical Analysis 2023-05-04 v1 Numerical Analysis Systems and Control Systems and Control Mathematical Physics math.MP Spectral Theory

Abstract

The first step when solving an infinite-dimensional eigenvalue problem is often to discretize it. We show that one must be extremely careful when discretizing nonlinear eigenvalue problems. Using examples, we show that discretization can: (1) introduce spurious eigenvalues, (2) entirely miss spectra, and (3) bring in severe ill-conditioning. While there are many eigensolvers for solving matrix nonlinear eigenvalue problems, we propose a solver for general holomorphic infinite-dimensional nonlinear eigenvalue problems that avoids discretization issues, which we prove is stable and converges. Moreover, we provide an algorithm that computes the problem's pseudospectra with explicit error control, allowing verification of computed spectra. The algorithm and numerical examples are publicly available in infNEP\texttt{infNEP}, which is a software package written in MATLAB.

Keywords

Cite

@article{arxiv.2305.01691,
  title  = {Avoiding discretization issues for nonlinear eigenvalue problems},
  author = {Matthew J. Colbrook and Alex Townsend},
  journal= {arXiv preprint arXiv:2305.01691},
  year   = {2023}
}