A GPU-Parallelized Interpolation-Based Fast Multipole Method for the Relativistic Space-Charge Field Calculation
arXiv:2302.07567 · doi:10.1016/j.cpc.2023.108825
Abstract
The fast multipole method (FMM) has received growing attention in the beam physics simulation. In this study, we formulate an interpolation-based FMM for the computation of the relativistic space-charge field. Different to the quasi-electrostatic model, our FMM is formulated in the lab-frame and can be applied without the assistance of the Lorentz transformation. In particular, we derive a modified admissibility condition which can effectively control the interpolation error of the proposed FMM. The algorithms and their GPU parallelization are discussed in detail. A package containing serial and GPU-parallelized solvers is implemented in the Julia programming language. The GPU-parallelized solver can reach a speedup of more than a hundred compared to the execution on a single CPU core.
arXiv comment30 pages, 10 figures
References in corpus (7)
- A Hierarchical O(N) Force Calculation Algorithm
- Effective Extensible Programming: Unleashing Julia on GPUs
- Biomolecular electrostatics using a fast multipole BEM on up to 512 GPUs and a billion unknowns
- Petascale turbulence simulation using a highly parallel fast multipole method on GPUs
- A GPU-Accelerated Fast Summation Method Based on Barycentric Lagrange Interpolation and Dual Tree Traversal
- Point-to-Point Coulomb Effects in High Brightness Photoelectron Beamlines for Ultrafast Electron Diffraction
- Relativistic Space-Charge Field Calculation by Interpolation-Based Treecode