SciPost Submission Page
An Algorithm to Parallelise Parton Showers on a GPU
by Michael H. Seymour, Siddharth Sule
Submission summary
| Authors (as registered SciPost users): | Siddharth Sule |
| Submission information | |
|---|---|
| Preprint Link: | https://arxiv.org/abs/2403.08692v2 (pdf) |
| Code repository: | https://gitlab.com/siddharthsule/gaps |
| Date submitted: | 2024-04-25 16:15 |
| Submitted by: | Sule, Siddharth |
| Submitted to: | SciPost Physics Codebases |
| Ontological classification | |
|---|---|
| Academic field: | Physics |
| Specialties: |
|
| Approaches: | Computational, Phenomenological |
Abstract
The Single Instruction, Multiple Thread (SIMT) paradigm of GPU programming does not support the branching nature of a parton shower algorithm by definition. However, modern GPUs are designed to schedule threads with diverging processes independently, allowing them to handle such branches. With regular thread synchronization and careful treatment of the individual steps, one can simulate a parton shower on a GPU. We present a parallelized Sudakov veto algorithm designed to simulate parton branching on multiple events simultaneously. We also release a CUDA C++ program that generates matrix elements, showers partons, and computes jet rates and event shapes for LEP at 91.2 GeV on a GPU. To benchmark its performance, we also provide a near-identical C++ program designed to simulate events serially on a CPU. While the consequences of branching are not absent, we demonstrate that a GPU can provide the throughput of a many-core CPU. As an example, we show that the time taken to simulate 10^6 events on one NVIDIA TESLA V100 GPU is equivalent to that of 258 Intel Xeon E5-2620 CPUs.