Gravitational Octree Code Performance Evaluation on Volta GPU

Abstract

In this study, the gravitational octree code originally optimized for the Fermi, Kepler, and Maxwell GPU architectures is adapted to the Volta architecture. The Volta architecture introduces independent thread scheduling requiring either the insertion of the explicit synchronizations at appropriate locations or the enforcement of the same implicit synchronizations as do the Pascal or earlier architectures by specifying -gencode arch=compute 60,code=sm 70. The performance measurements on Tesla V100, the current flagship GPU by NVIDIA, revealed that the N-body simulations of the Andromeda galaxy model with 2 23 = 8 388 608 particles took 3.8 × 10 −2 s or 3.3 × 10 −2 s per step for each case. Tesla V100 achieves a 1.4 to 2.2-fold acceleration in comparison with Tesla P100, the flagship GPU in the previous generation. The observed speed-up of 2.2 is greater than 1.5, which is the ratio of the theoretical peak performance of the two GPUs. The independence of the units for integer operations from those for floatingpoint number operations enables the overlapped execution of integer and floating-point number operations. It hides the execution time of the integer operations leading to the speed-up rate above the theoretical peak performance ratio. Tesla V100 can execute N-body simulation with up to 25 × 2 20 = 26 214 400 particles, and it took 2.0 × 10 −1 s per step. It corresponds to 3.5 TFlop/s, which is 22% of the single-precision theoretical peak performance.

B-spline-based material point method with dynamic load balancing technique for large-scale simulation

“…To address this issue and mitigate the imbalance in computation and communication as much as possible, especially in scenarios involving largescale and long-duration phenomena, we adopted the dynamic load balancing (DLB) technique proposed by [47] in this study. The implementation was initially developed for the gravitational octree code GOTHIC [48,49] and adjusted for the MPM in this study. Domain decomposition in particle methods generally does not necessarily concern the geometry of each domain.…”