A continuation approach to the computation of essential and absolute spectra of differential operators on the real line is presented. The advantages of this approach, compared with direct eigenvalue computations for the discretized operator, are the efficient and accurate computation of selected parts of the spectrum (typically those near the imaginary axis) and the option to compute nonlinear travelling waves and selected eigenvalues or other stability indicators simultaneously in order to locate accurately the onset to instability. We also discuss the implementation and usage of this approach with the software package auto and provide example computations for the FitzHugh-Nagumo and the complex Ginzburg-Landau equation.
Bibliographical noteFunding Information:
J. Rademacher acknowledges hospitality at the Free University of Berlin and support from NSF grant DMS-0203301, a PIMS fellowship and DFG priority program SPP 1095. B. Sandstede acknowledges support from the NSF through grant DMS-0203854 and from a Royal Society Wolfson Research Merit Award. A. Scheel was partially supported by the NSF through grants DMS-0203301 and DMS-0504271.
- Absolute spectrum
- Instability thresholds
- Reaction-diffusion systems
- Spectral stability