Spectral algorithms for reaction-diffusion equations
This is a practical resource for practitioners needing efficient numerical solutions for reaction-diffusion equations, but it is an incremental contribution as it packages existing methods.
The paper provides MATLAB and Fortran 77 codes for solving reaction-diffusion equations using spectral methods to handle stiffness, enabling explicit high-order timestepping. It includes examples, timings, and error comparisons with standard methods.
A collection of codes (in MATLAB & Fortran 77), and examples, for solving reaction-diffusion equations in one and two space dimensions is presented. In areas of the mathematical community spectral methods are used to remove the stiffness associated with the diffusive terms in a reaction-diffusion model allowing explicit high order timestepping to be used. This is particularly valuable for two (and higher) space dimension problems. Our aim here is to provide codes, together with examples, to allow practioners to easily utilize, understand and implement these ideas; we incorporate recent theoretical advances such as exponential time differencing methods and provide timings and error comparisons with other more standard approaches. The examples are chosen from the literature to illustrate points and queries that naturally arise.