Evaluation of Single-Node Performance of Parallel Algorithms for Multigroup Monte Carlo Particle Transport Methods

Abstract

Monte Carlo (MC) methods have been widely used to solve the particle transport equation due to their high accuracy and capability of processing complex geometries. History-based and event-based algorithms that are applicable to different architectures are two methods for parallelizing the MC code. There is a large work on evaluating and optimizing parallel algorithms with continuous-energy schemes. In this work, we evaluate the single-node performance of history-based and event-based algorithms for multigroup MC methods on both CPUs and GPUs with Quicksilver, a multigroup MC transport code that has already implemented the history-based algorithms. We first implement and optimize the event-based algorithm based on Quicksilver and then perform the evaluation work extensively on the Coral2 benchmark. Numerical results indicate that contrary to continuous-energy schemes, the history-based approach with multigroup schemes outperforms the event-based algorithm on both architectures in all cases. We summarize that the performance loss of the event-based algorithm is mainly due to: 1) extra operations to reorganize particles, 2) batched atomic operations, and 3) poor particle data locality. Despite the poor performance, the event-based algorithm achieves higher memory bandwidth utilization. We further discuss the impact of memory access patterns and calculation of cross sections (xs) on the performance of the GPU. Built on the analytics, and shed light on the algorithm choice and optimizations for paralleling the MC transport code on different architectures.

Monte Carlo neutron transport using low power mobile GPU devices

“…There are many GPU implementation of Monte Carlo methods [1] [10]. Here discusses one important optimization.…”

“…Therefore, the more powerful graphics processing units (GPU) gradually becomes a powerful tool to accelerate the Monte Carlo method. Existing researches [1] [10] have optimized the GPU algorithms in many aspects, and achieved significant progresses. One of the important method worthy of mention is the optimization method of particle sorting.…”

Exploring Numba and CuPy for GPU-Accelerated Monte Carlo Radiation Transport

“…Since the MCRT utilizes a large number of MC particles that are evolved independently from each other (at least, with a timestep), the many-core architecture can evolve particles in parallel, leading to a significant speed-up compared to serial calculations [57][58][59]. There have been many studies and applications of GPU-accelerated MCRT [60][61][62][63]. The generation of random numbers on GPUs was investigated by, e.g., [64][65][66].…”

Positron Interactions with Some Human Body Organs Using Monte Carlo Probability Method

“…Microscopic collisions of charged particles as many-body problems cannot be solved accurately [10], and enough of the theories apply to several assumptions and approximations. The Monte Carlo estimates are an excellent method for analyzing the particle transport in matter [11].…”