Semi-Lagrangian Vlasov simulation on GPUs
arXiv:1907.08316 · doi:10.1016/j.cpc.2020.107351
Abstract
In this paper, our goal is to efficiently solve the Vlasov equation on GPUs. A semi-Lagrangian discontinuous Galerkin scheme is used for the discretization. Such kinetic computations are extremely expensive due to the high-dimensional phase space. The SLDG code, which is publicly available under the MIT license abstracts the number of dimensions and uses a shared codebase for both GPU and CPU based simulations. We investigate the performance of the implementation on a range of both Tesla (V100, Titan V, K80) and consumer (GTX 1080 Ti) GPUs. Our implementation is typically able to achieve a performance of approximately 470 GB/s on a single GPU and 1600 GB/s on four V100 GPUs connected via NVLink. This results in a speedup of about a factor of ten (comparing a single GPU with a dual socket Intel Xeon Gold node) and approximately a factor of 35 (comparing a single node with and without GPUs). In addition, we investigate the effect of single precision computation on the performance of the SLDG code and demonstrate that a template based dimension independent implementation can achieve good performance regardless of the dimensionality of the problem.
References in corpus (9)
- A positivity-preserving high-order semi-Lagrangian discontinuous Galerkin scheme for the Vlasov-Poisson equations
- Vlasov methods in space physics and astrophysics
- A semi-Lagrangian Vlasov solver in tensor train format
- A low-rank projector-splitting integrator for the Vlasov--Maxwell equations with divergence correction
- A comparison of semi-Lagrangian discontinuous Galerkin and spline based Vlasov solvers in four dimensions
- Block-Structured Adaptive Mesh Refinement Algorithms for Vlasov Simulation
- Convergence analysis of a discontinuous Galerkin/Strang splitting approximation for the Vlasov--Poisson equations
- High performance computing aspects of a dimension independent semi-Lagrangian discontinuous Galerkin code
- A study on conserving invariants of the Vlasov equation in semi-Lagrangian computer simulations
Cited by in corpus (5)
- A quantum-inspired method for solving the Vlasov-Poisson equations
- Exponential methods for solving hyperbolic problems with application to kinetic equations
- An accurate and time-parallel rational exponential integrator for hyperbolic and oscillatory PDEs
- A -mode integrator for solving evolution equations in Kronecker form
- A performance portable implementation of the semi-Lagrangian algorithm in six dimensions