A Parallel Exponential Integrator for Large-Scale Discretizations of Advection-Diffusion Models

We propose a parallel implementation of the ReLPM (Real Leja Points Method) for the exponential integration of large sparse systems of ODEs, generated by Finite Element discretizations of 3D advection-diffusion models. The performance of our parallel exponential integrator is compared with that of a parallelized Crank-Nicolson (CN) integrator, where the local linear solver is a parallel BiCGstab accelerated with the approximate inverse preconditioner FSAI. We developed message passing codes written in Fortran 90 and using the MPI standard. Results on SP5 and CLX machines show that the parallel efficiency raised by the two algorithms is comparable. ReLPM turns out to be from 3 to 5 times faster than CN in solving realistic advection-diffusion problems, depending on the number of processors employed.

[1]  Message P Forum,et al.  MPI: A Message-Passing Interface Standard , 1994 .

[2]  Marco Vianello,et al.  The ReLPM Exponential Integrator for FE Discretizations of Advection-Diffusion Equations , 2004, International Conference on Computational Science.

[3]  L. Bergamaschi,et al.  Interpolating discrete advection-diffusion propagators at Leja sequences , 2004 .

[4]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..

[5]  Marlis Hochbruck,et al.  Exponential Integrators for Large Systems of Differential Equations , 1998, SIAM J. Sci. Comput..

[6]  L. Kolotilina,et al.  Factorized Sparse Approximate Inverse Preconditionings I. Theory , 1993, SIAM J. Matrix Anal. Appl..

[7]  Luca Bergamaschi,et al.  Parallel Acceleration of Krylov Solvers by Factorized Approximate Inverse Preconditioners , 2004, VECPAR.

[8]  Luca Bergamaschi,et al.  Efficient Parallelization of Preconditioned Conjugate Gradient Scheme for Matrices Arising from Discretizations of Diffusion Equations , 1999, PPSC.

[9]  Message Passing Interface Forum.  MPI: A message - passing interface standard , 1994 .

[10]  A. A. Nikishin,et al.  Factorized sparse approximate inverse preconditionings. III. Iterative construction of preconditioners , 2000 .

[11]  Jack Dongarra,et al.  High Performance Computing for Computational Science - VECPAR 2004, 6th International Conference, Valencia, Spain, June 28-30, 2004, Revised Selected and Invited Papers , 2005, VECPAR.

[12]  Stephen Gilmore,et al.  Evaluating the Performance of Skeleton-Based High Level Parallel Programs , 2004, International Conference on Computational Science.