Massively parallel computation of atmospheric neutrino oscillations on CUDA enabled accelerators
Abstract The computation of neutrino flavor transition amplitudes through inhomogeneous matter is a time consuming step and thus could benefit from optimization and parallelization. Next to reliable parameter estimation of intrinsic physical quantities such as neutrino masses and mixing angles, these transition amplitudes are important in hypothesis testing of potential extensions of the standard model of elementary particle physics, such as additional neutrino flavors. Hence, fast yet precise implementations are of high importance to research. In the recent past, massively parallel accelerators such as CUDA enabled GPUs featuring thousands of compute units have been widely adopted due to their superior memory bandwidth, vast compute capability, and highly competitive compute to energy ratio in comparison to traditional multi core architectures with a few tens of monolithic cores. In this paper, we introduce two scalable multi GPU extensions of common neutrino oscillation frameworks namely Prob3+ and n SQuIDS allowing for the acceleration of oscillation dynamics computation by one to three orders of magnitude while preserving numerical accuracy. Our software is licensed under LGPLv3 and can be accessed at https:/github.com/fkallen/CUDAProb3 and https:/github.com/fkallen/CUDAnuSQuIDS Program summary Program Title: CUDAProb3+ and CUDA n SQuIDS Program Files doi: http:/dx.doi.org/10.17632/j54yymmg5h.1 Licensing provisions: LGPLv3 Programming language: CUDA, C and C+ Nature of problem: Fast computation of neutrino oscillation probabilities in inhomogeneous matter using constant coefficient or variable coefficient ODE solvers. Solution method: Solving ordinary differential equations on CUDA enabled GPUs.
Comput. Phys. Commun.