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

Abstract

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.

arXiv comment49 pages

References in corpus (10)

• Second Order Accurate Schemes for Magnetohydrodynamics With Divergence-Free Reconstruction
• Divergence-Free Adaptive Mesh Refinement for Magnetohydrodynamics
• Efficient, High Accuracy ADER-WENO Schemes for Hydrodynamics and Divergence-Free Magnetohydrodynamics
• A Divergence-Free Upwind Code for Multidimensional Magnetohydrodynamic Flows
• Multidimensional HLLE Riemann solver; Application to Euler and Magnetohydrodynamic Flows
• ADER-WENO Finite Volume Schemes with Space-Time Adaptive Mesh Refinement
• A Two-dimensional HLLC Riemann Solver for Conservation Laws : Application to Euler and MHD Flows
• Divergence-Free Reconstruction of Magnetic Fields and WENO Schemes for Magnetohydrodynamics
• A High-Order Relativistic Two-Fluid Electrodynamic Scheme with Consistent Reconstruction of Electromagnetic Fields and a Multidimensional Riemann Solver for Electromagnetism
• Modeling Magnetically Channeled Winds in 3D: I. Isothermal Simulations of a Magnetic O Supergiant