Exponential Integrators on Graphic Processing Units

Abstract

In this paper we revisit stencil methods on GPUs in the context of exponential integrators. We further discuss boundary conditions, in the same context, and show that simple boundary conditions (for example, homogeneous Dirichlet or homogeneous Neumann boundary conditions) do not affect the performance if implemented directly into the CUDA kernel. In addition, we show that stencil methods with position-dependent coefficients can be implemented efficiently as well. As an application, we discuss the implementation of exponential integrators for different classes of problems in a single and multi GPU setup (up to 4 GPUs). We further show that for stencil based methods such parallelization can be done very efficiently, while for some unstructured matrices the parallelization to multiple GPUs is severely limited by the throughput of the PCIe bus.

To appear in: Proceedings of the 2013 International Conference on High Performance Computing Simulation (HPCS 2013), IEEE (2013)

Cited by in corpus (4)

• On the performance of exponential integrators for problems in magnetohydrodynamics
• Exponential methods for solving hyperbolic problems with application to kinetic equations
• A conservative discontinuous Galerkin scheme for the 2D incompressible Navier--Stokes equations
• An accurate and time-parallel rational exponential integrator for hyperbolic and oscillatory PDEs