Abstract
We implemented two- and three-dimensional Rayleigh–Benard convection on Nvidia GPUs by utilizing a 2nd-order finite difference method. By exploiting the massive parallelism of GPU using both CUDA for C and optimized CUBLAS routines, we have on a single Fermi GPU run simulations of Rayleigh number up to 6 × 1010 (on a mesh of 2000 × 4000 uniform grid points) in two dimensions and up to 107 (on a mesh of 450 × 450 × 225 uniform grid points) for three dimensions. On Nvidia Tesla C2070 GPUs, these implementations enjoy single-precision performance of 535 GFLOP/s and 100 GFLOP/s respectively, and double-precision performance of 230 GFLOP/s and 70 GFLOP/s respectively.
Original language | English (US) |
---|---|
Title of host publication | Lecture Notes in Earth System Sciences |
Publisher | Springer International Publishing |
Pages | 335-352 |
Number of pages | 18 |
Edition | 9783642164040 |
DOIs | |
State | Published - 2013 |
Publication series
Name | Lecture Notes in Earth System Sciences |
---|---|
Number | 9783642164040 |
Volume | 0 |
ISSN (Print) | 2193-8571 |
ISSN (Electronic) | 2193-858X |
Bibliographical note
Funding Information:Fig.22.9 Nuhistogramcorrespondingtothedatain(Fig.22.3).Ra3×1010isinblueand6×1010 is in purple Acknowledgments We thank Matt Knepley for stimulating discussions on GPU. This research has been supported by NSF CMG grant.
Publisher Copyright:
© Springer-Verlag Berlin Heidelberg 2013.