Help Privacy Policy Disclaimer
  Advanced SearchBrowse




Journal Article

Optimization of the ADER-DG method in GPU applied to linear hyperbolic PDEs

There are no MPG-Authors in the publication available
External Resource
No external resources are shared
Fulltext (restricted access)
There are currently no full texts shared for your IP range.
Fulltext (public)
There are no public fulltexts stored in PuRe
Supplementary Material (public)
There is no public supplementary material available

Castro, C. E., Behrens, J., & Pelties, C. (2016). Optimization of the ADER-DG method in GPU applied to linear hyperbolic PDEs. International Journal for Numerical Methods in Fluids, 81, 195-219. doi:10.1002/fld.4179.

Cite as: https://hdl.handle.net/11858/00-001M-0000-002D-A2B6-3
We optimized the Arbitrary accuracy DErivatives Riemann problem (ADER) - Discontinuous Galerkin (DG) numerical method using the CUDA-C language to run the code in a graphic processing unit (GPU). We focus on solving linear hyperbolic partial-differential equations where the method can be expressed as a combination of precomputed matrix multiplications becoming a good candidate to be used on the GPU hardware. Moreover, the method is arbitrarily high order involving intensive work on local data, a property that is also beneficial for the target hardware. We compare our GPU implementation against CPU versions of the same method observing similar convergence properties up to a threshold where the error remains fixed. This behavior is in agreement with the CPU version, but the threshold is slightly larger than in the CPU case. We also observe a big difference when considering single and double precisions where in the first case, the threshold error is significantly larger. Finally, we did observe a speed-up factor in computational time that depends on the order of the method and the size of the problem. In the best case, our novel GPU implementation runs 23 times faster than the CPU version. We used three partial-differential equation to test the code considering the linear advection equation, the seismic wave equation, and the linear shallow water equation, all of them considering variable coefficients. Copyright (c) 2015 John Wiley & Sons, Ltd.