This paper describes the implementation of a solver for systems of semi-linear parabolic partial differential equations in two and three space dimensions. The solver is based on a parallel implementation of a non-linear Alternating Direction Implicit (ADI) scheme which uses a Cartesian grid in space and an implicit time-stepping algorithm. Various reordering strategies for the linearized equations are used to reduce the stride and improve the overall effectiveness of the parallel implementation. We have successfully used this solver for large-scale reaction-diffusion problems in computational biology and medicine in which the desired solution is a traveling wave that may contain rapid transitions. A number of examples that illustrate the efficiency and accuracy of the method are given here; the theoretical analysis will be presented.
[1]
P. Monk.
Two Dimensional Wave Propagation in a Model of Dictyostelium Discoideum
,
1991
.
[2]
Hans Othmer,et al.
A computational study of wave propagation in a model for anisotropic cardiac ventricular tissue
,
1995,
HPCN Europe.
[3]
S. Lennart Johnsson,et al.
Alternating direction methods on multiprocessors
,
1987
.
[4]
Robert Dillon,et al.
Pattern formation in generalized Turing systems
,
1994
.
[5]
Yousef Saad,et al.
Hybrid Krylov Methods for Nonlinear Systems of Equations
,
1990,
SIAM J. Sci. Comput..
[6]
Mark J. Ablowitz,et al.
Explicit solutions of Fisher's equation for a special wave speed
,
1979
.