Preprint des Projektes: SFB 701: Spektrale Strukturen und Topologische Methoden in der Mathematik - Projekt B3
Numerische Analyse äquivarianter Evolutionsgleichungen
11-024 Wolf-Jürgen Beyn, Cedric Effenberger, Daniel Kressner.
Invariant pairs have been proposed as a numerically robust means to represent and compute several eigenvalues along with the corresponding (generalized) eigenvectors for matrix eigenvalue problems that are nonlinear in the eigenvalue parameter. In this work, we consider nonlinear eigenvalue problems that depend on an additional parameter and our interest is to track several eigenvalues as this parameter varies. Based on the concept of invariant pairs, a theoretically sound and reliable numerical continuation procedure is developed. Particular attention is paid to the situation when the procedure approaches a singularity, that is, when eigenvalues included in the invariant pair collide with other eigenvalues. For the real generic case, it is proven that such a singularity only occurs when two eigenvalues collide on the real axis. It is shown how this situation can be handled numerically by an appropriate expansion of the invariant pair. Several numerical examples demonstrate the viability of our continuation procedure.