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
