The code is a versatile tool for numerical simulations of magnetohydrodynamic (MHD) processes in solar and stellar atmospheres. This paper describes the implementation of a CUDA version of the code, which enables it to run efficiently on modern supercomputers with GPUs. We carefully studied the numerical performance, with detailed comparisons to the CPU version, together with strong and weak scaling of the code. Low-level CUDA optimizations allowed us to reach an astonishing 600-fold higher performance with respect to the CPU version. Weak scaling, along with tuned communications between GPUs, enables the yield of more than 85% of the performance for up to 400 GPUs.
Modestov et al. (Tue,) studied this question.