Implementation and Optimization of a 1D2V PIC Method for Nonlinear Kinetic Models on GPUs

Title data

Korch, Matthias ; Raithel, Philipp ; Werner, Tim:
Implementation and Optimization of a 1D2V PIC Method for Nonlinear Kinetic Models on GPUs.
In: 2020 28th Euromicro International Conference on Parallel, Distributed and Network-Based Processing : PDP 2020 ; Proceedings. - Piscataway, NJ : IEEE , 2020 . - pp. 30-37
ISBN 978-1-7281-6582-0
DOI: https://doi.org/10.1109/PDP50117.2020.00012

Project information

Abstract in another language

This paper considers the parallel numerical simulation of the time evolution of galaxies and globular clusters on GPUs. The model used is the Einstein-Vlasov system, which is designed, in particular, to study the formation of black holes and spacetime singularities in a general relativistic framework. First, a reference implementation is derived using NVIDIA CUDA as programming model, which is then optimized in several steps. Bottlenecks are identified by profiling, and different approaches, namely particle sort, improved treatment of atomic operations, and kernel fusion are investigated to overcome these bottlenecks. Each optimized variant is evaluated in relation to the other variants using detailed runtime experiments and profiling results. Using in the order of 10^7 to 10^8 particles, speedups between 1.84 and 2.38 w.r.t. the reference implementation have been observed.