Accelerating a 3D finite-difference wave propagation code using GPU graphics cards

David Michéa 1 Dimitri Komatitsch 2
2 Magique 3D - Advanced 3D Numerical Modeling in Geophysics
LMAP - Laboratoire de Mathématiques et de leurs Applications [Pau], Inria Bordeaux - Sud-Ouest
Abstract : We accelerate a three-dimensional finite-difference in the time domain (FDTD) wave propagation code by a factor between about 20 and 60 compared to a serial implementation using Graphics Processing Unit (GPU) computing on NVIDIA graphics cards with the CUDA programming language. We describe the implementation of the code in CUDA to simulate the propagation of seismic waves in a heterogeneous elastic medium. We also implement Convolution Perfectly Matched Layers (CPMLs) on the graphics cards to efficiently absorb outgoing waves on the fictitious edges of the grid. We show that the code that runs on a graphics card gives the expected results by comparing our results to those obtained by running the same simulation on a classical processor core. The methodology that we present can be used for Maxwell's equations as well because their form is similar to that of the seismic wave equation written in velocity vector and stress tensor.
Complete list of metadatas

Cited literature [72 references]  Display  Hide  Download
Contributor : Dimitri Komatitsch <>
Submitted on : Tuesday, November 30, 2010 - 1:04:47 PM
Last modification on : Thursday, May 9, 2019 - 9:40:04 AM
Long-term archiving on : Tuesday, March 1, 2011 - 2:27:31 AM


Files produced by the author(s)




David Michéa, Dimitri Komatitsch. Accelerating a 3D finite-difference wave propagation code using GPU graphics cards. Geophysical Journal International, Oxford University Press (OUP), 2010, 182 (1), pp.389-402. ⟨10.1111/j.1365-246X.2010.04616.x⟩. ⟨inria-00528487⟩



Record views


Files downloads