High rayleigh number mantle convection on GPU

David A. Sanchez, Christopher Gonzalez, David A. Yuen, Grady B. Wright, Gregory A. Barnett

Research output: Chapter in Book/Report/Conference proceedingChapter


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 languageEnglish (US)
Title of host publicationLecture Notes in Earth System Sciences
PublisherSpringer International Publishing
Number of pages18
StatePublished - 2013

Publication series

NameLecture Notes in Earth System Sciences
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.


Dive into the research topics of 'High rayleigh number mantle convection on GPU'. Together they form a unique fingerprint.

Cite this