cuGMEC: A high-performance code for gyrokinetic-MHD hybrid simulation on GPUs with CUDA C++

To numerically simulate EP-driven Alfvén instabilities, we have developed GMEC[1,2] (Gyrokinetic-Magnetohydrodynamics Energetic-particle Code) on the central processing unit (CPU) platform. However, long-time and large-scale simulations require higher grid resolution and more particles, demanding substantial computational resources.

Therefore, applying more efficient numerical algorithms and optimization strategies on advanced computing architectures becomes critical. Computer scientists have designed the graphics processing unit (GPU) to implement single instruction multiple thread (SIMT) with fewer control units and higher throughput than CPUs. The GPU is optimized for data-parallel tasks and has emerged as a powerful platform for scientific computing.