TY - JOUR

T1 - Numerical evaluation of the evans function by Magnus integration

AU - Aparicio, Nairo D.

AU - Malham, S. J A

AU - Oliver, Marcel

PY - 2005/6

Y1 - 2005/6

N2 - We use Magnus methods to compute the Evans function for spectral problems as arise when determining the linear stability of travelling wave solutions to reaction-diffusion and related partial differential equations. In a typical application scenario, we need to repeatedly sample the solution to a system of linear non-autonomous ordinary differential equations for different values of one or more parameters as we detect and locate the zeros of the Evans function in the right half of the complex plane. In this situation, a substantial portion of the computational effort-the numerical evaluation of the iterated integrals which appear in the Magnus series-can be performed independent of the parameters and hence needs to be done only once. More importantly, for any given tolerance Magnus integrators possess lower bounds on the step size which are uniform across large regions of parameter space and which can be estimated a priori. We demonstrate, analytically as well as through numerical experiment, that these features render Magnus integrators extremely robust and, depending on the regime of interest, efficient in comparison with standard ODE solvers. © Springer 2005.

AB - We use Magnus methods to compute the Evans function for spectral problems as arise when determining the linear stability of travelling wave solutions to reaction-diffusion and related partial differential equations. In a typical application scenario, we need to repeatedly sample the solution to a system of linear non-autonomous ordinary differential equations for different values of one or more parameters as we detect and locate the zeros of the Evans function in the right half of the complex plane. In this situation, a substantial portion of the computational effort-the numerical evaluation of the iterated integrals which appear in the Magnus series-can be performed independent of the parameters and hence needs to be done only once. More importantly, for any given tolerance Magnus integrators possess lower bounds on the step size which are uniform across large regions of parameter space and which can be estimated a priori. We demonstrate, analytically as well as through numerical experiment, that these features render Magnus integrators extremely robust and, depending on the regime of interest, efficient in comparison with standard ODE solvers. © Springer 2005.

KW - Evans function

KW - Magnus expansion

KW - Neumann expansion

UR - http://www.scopus.com/inward/record.url?scp=27144551687&partnerID=8YFLogxK

U2 - 10.1007/s10543-005-0001-8

DO - 10.1007/s10543-005-0001-8

M3 - Article

VL - 45

SP - 219

EP - 258

JO - BIT Numerical Mathematics

JF - BIT Numerical Mathematics

SN - 0006-3835

IS - 2

ER -