Abstract:Lagrangian particle dispersion models are indispensable tools for the study of atmospheric transport processes. However, Lagrangian transport simulations can become numerically expensive when large numbers of air parcels are involved. To accelerate these simulations, we made considerable efforts to port the Massive-Parallel Trajectory Calculations (MPTRAC) model to graphics processing units (GPUs). Here we discuss performance optimizations of the major bottleneck of the GPU code of MPTRAC, the advection kernel. Timeline, roofline, and memory analyses of the baseline GPU code revealed that the application is memory-bound, and performance suffers from near-random memory access patterns. By changing the data structure of the horizontal wind and vertical velocity fields of the global meteorological data driving the simulations from structure of arrays (SoAs) to array of structures (AoSs) and by introducing a sorting method for better memory alignment of the particle data, performance was greatly improved. We evaluated the performance on NVIDIA A100 GPUs of the Jülich Wizard for European Leadership Science (JUWELS) Booster module at the Jülich Supercomputing Center, Germany. For our largest test case, transport simulations with 10 8 particles driven by the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA5 reanalysis, we found that the runtime for the full set of physics computations was reduced by 75 %, including a reduction of 85 % for the advection kernel. In addition to demonstrating the benefits of code optimization for GPUs, we show that the runtime of central processing unit (CPU-)only simulations is also improved. For our largest test case, we found a runtime reduction of 34 % for the physics computations, including a reduction of 65 % for the advection kernel. The code optimizations discussed here bring the MPTRAC model closer to applications on upcoming exascale high-performance computing systems and will also be of interest for optimizing the performance of other models using particle methods.

Speeding up MadGraph5_aMC@NLO

Speeding up Madgraph5 aMC@NLO through CPU vectorization and GPU offloading: towards a first alpha release

Acceleration beyond lowest order event generation: An outlook on further parallelism within MadGraph5_aMC@NLO

Madgraph5_aMC@NLO on GPUs and vector CPUs Experience with the first alpha release

aMCfast: automation of fast NLO computations for PDF fits

The MadNIS Reloaded

Simplifying QCD event generation with chirality flow, reference vectors and spin directions

Fast Computation of MadGraph Amplitudes on Graphics Processing Unit (GPU)

Asymmetric collisions in MadGraph5_aMC@NLO

Accelerating Lagrangian transport simulations on graphics processing units: performance optimizations of Massive-Parallel Trajectory Calculations (MPTRAC) v2.6

Optimizations of the Accelerated Axial Polynomial Method of Characteristics of the APOLLO3® Deterministic Neutron Transport Code

Optimizing Near Field Computation in the MLFMA Algorithm with Data Redundancy and Performance Modeling on a Single GPU

Quadratic Speed-up in Infinite Variance Quantum Monte Carlo

Multi-GPU RI-HF Energies and Analytic Gradients $-$ Towards High Throughput Ab Initio Molecular Dynamics

Embracing a new era of highly efficient and productive quantum Monte Carlo simulations

Accelerating Lattice QCD Simulations using GPUs

General Method for Speeding Up Kinetic Monte Carlo Simulations

Accelerating VASP electronic structure calculations using graphic processing units

Multi-GPU RI-HF Energies and Analytic Gradients─Toward High-Throughput Ab Initio Molecular Dynamics

Resolved photoproduction in {\tt MadGraph5\_aMC@NLO}