This work considers the Real Leja Points Method (ReLPM) for the exponential integration of large-scale sparse systems of ODEs, generated by Finite Element or Finite Difference discretizations of 3-D advection-diffusion models. We present an efficient parallel implementation of ReLPM for polynomial interpolation of the matrix exponential propagators. A scalability analysis of the most important computational kernel inside the code, the parallel sparse matrix–vector product, has been performed, as well as an experimental study of the communication overhead. As a result of this study an optimized parallel sparse matrix–vector product routine has been implemented. The resulting code shows good scaling behavior even when using more than one thousand processors. The numerical results presented on a number of very large test cases gives experimental evidence that ReLPM is a reliable and efficient tool for the simulation of complex hydrodynamic processes on parallel architectures.

A massively parallel exponential integrator for advection-diffusion models

CALIARI, Marco;
2009-01-01

Abstract

This work considers the Real Leja Points Method (ReLPM) for the exponential integration of large-scale sparse systems of ODEs, generated by Finite Element or Finite Difference discretizations of 3-D advection-diffusion models. We present an efficient parallel implementation of ReLPM for polynomial interpolation of the matrix exponential propagators. A scalability analysis of the most important computational kernel inside the code, the parallel sparse matrix–vector product, has been performed, as well as an experimental study of the communication overhead. As a result of this study an optimized parallel sparse matrix–vector product routine has been implemented. The resulting code shows good scaling behavior even when using more than one thousand processors. The numerical results presented on a number of very large test cases gives experimental evidence that ReLPM is a reliable and efficient tool for the simulation of complex hydrodynamic processes on parallel architectures.
2009
advection diffusion equations; exponential integrators; parallel implementation
File in questo prodotto:
File Dimensione Formato  
preMBCV09.pdf

accesso aperto

Tipologia: Documento in Pre-print
Licenza: Dominio pubblico
Dimensione 960.95 kB
Formato Adobe PDF
960.95 kB Adobe PDF Visualizza/Apri
MBCV09.pdf

non disponibili

Tipologia: Versione dell'editore
Licenza: Accesso ristretto
Dimensione 928.21 kB
Formato Adobe PDF
928.21 kB Adobe PDF   Visualizza/Apri   Richiedi una copia

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11562/470966
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 28
  • ???jsp.display-item.citation.isi??? 26
social impact