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.


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

Cite this