High-order discontinuous galerkin methods by GPU metaprogramming

Andreas Klöckner, Timothy Warburton, Jan S. Hesthaven

Research output: Chapter in Book/Report/Conference proceedingChapter


Discontinuous Galerkin (DG) methods for the numerical solution of partial differential equations have enjoyed considerable success because they are both flexible and robust: They allow arbitrary unstructured geometries and easy control of accuracy without compromising simulation stability. In a recent publication, we have shown that DG methods also adapt readily to execution on modern, massively parallel graphics processors (GPUs). A number of qualities of the method contribute to this suitability, reaching from locality of reference, through regularity of access patterns, to high arithmetic intensity. In this article, we illuminate a few of the more practical aspects of bringing DG onto a GPU, including the use of a Python-based metaprogramming infrastructure that was created specifically to support DG, but has found many uses across all disciplines of computational science.

Original languageEnglish (US)
Title of host publicationLecture Notes in Earth System Sciences
Number of pages22
StatePublished - 2013
Externally publishedYes

Publication series

NameLecture Notes in Earth System Sciences
ISSN (Print)2193-8571
ISSN (Electronic)2193-858X

ASJC Scopus subject areas

  • Computers in Earth Sciences
  • General Earth and Planetary Sciences


Dive into the research topics of 'High-order discontinuous galerkin methods by GPU metaprogramming'. Together they form a unique fingerprint.

Cite this