## 3D finite difference computation on GPUs using CUDA

NVIDIA, 2701 San Tomas Expressway, Santa Clara, CA 95050

In GPGPU-2: Proceedings of 2nd Workshop on General Purpose Processing on Graphics Processing Units (2009), pp. 79-84.

@conference{micikevicius20093d,

title={3D finite difference computation on GPUs using CUDA},

author={Micikevicius, P.},

booktitle={Proceedings of 2nd Workshop on General Purpose Processing on Graphics Processing Units},

pages={79–84},

year={2009},

organization={ACM}

}

In this paper we describe a GPU parallelization of the 3D finite difference computation using CUDA. Data access redundancy is used as the metric to determine the optimal implementation for both the stencil-only computation, as well as the discretization of the wave equation, which is currently of great interest in seismic computing. For the larger stencils, the described approach achieves the throughput of between 2,400 to over 3,000 million of output points per second on a single Tesla 10-series GPU. This is roughly an order of magnitude higher than a 4-core Harpertown CPU running a similar code from seismic industry. Multi-GPU parallelization is also described, achieving linear scaling with GPUs by overlapping inter-GPU communication with computation.

November 2, 2010 by hgpu

Similar papers:

- Modeling the propagation of elastic waves using spectral elements on a cluster of 192 GPUs
- High-order finite-element seismic wave propagation modeling with MPI on a large GPU cluster
- Scalable Multi-GPU 3-D FFT for TSUBAME 2.0 Supercomputer
- Polymer Field-Theory Simulations on Graphics Processing Units
- Solving 3D Anisotropic Elastic Wave Equations on Parallel GPU Devices