Techniques, Tricks and Algorithms for Efficient GPU-Based Processing of Higher Order Hyperbolic PDEs

by   Sethupathy Subramanian, et al.

GPU computing is expected to play an integral part in all modern Exascale supercomputers. It is also expected that higher order Godunov schemes will make up about a significant fraction of the application mix on such supercomputers. It is, therefore, very important to prepare the community of users of higher order schemes for hyperbolic PDEs for this emerging opportunity. We focus on three broad and high-impact areas where higher order Godunov schemes are used. The first area is computational fluid dynamics (CFD). The second is computational magnetohydrodynamics (MHD) which has an involution constraint that has to be mimetically preserved. The third is computational electrodynamics (CED) which has involution constraints and also extremely stiff source terms. Together, these three diverse uses of higher order Godunov methodology, cover many of the most important applications areas. In all three cases, we show that the optimal use of algorithms, techniques and tricks, along with the use of OpenACC, yields superlative speedups on GPUs! As a bonus, we find a most remarkable and desirable result: some higher order schemes, with their larger operations count per zone, show better speedup than lower order schemes on GPUs. In other words, the GPU is an optimal stratagem for overcoming the higher computational complexities of higher order schemes! Several avenues for future improvement have also been identified. A scalability study is presented for a real-world application using GPUs and comparable numbers of high-end multicore CPUs. It is found that GPUs offer a substantial performance benefit over comparable number of CPUs, especially when all the methods designed in this paper are used.


page 1

page 2

page 3

page 4


Fast sparse grid simulations of fifth order WENO scheme for high dimensional hyperbolic PDEs

The weighted essentially non-oscillatory (WENO) schemes are a popular cl...

Higher-order magnetohydrodynamic numerics

In this chapter, we aim at presenting the basic techniques necessary to ...

Higher Order Probabilities

A number of writers have supposed that for the full specification of bel...

On Higher-Order Cryptography (Long Version)

Type-two constructions abound in cryptography: adversaries for encryptio...

Higher-order Segmentation via Multicuts

Multicuts enable to conveniently represent discrete graphical models for...

The sequence of higher order Mersenne numbers and associated binomial transforms

In this article, we introduce and study a new integer sequence referred ...

Iso-geometric Integral Equation Solvers and their Compression via Manifold Harmonics

The state of art of electromagnetic integral equations has seen signific...

Please sign up or login with your details

Forgot password? Click here to reset