We present preliminary performance results of gPLUTO, the new GPU-optimized implementation of the PLUTO code for computational plasma astrophysics. Like its predecessor, gPLUTO employs a finite-volume formulation to numerically solve the equations of magnetohydrodynamics (MHD) in multiple spatial dimensions. Still, this new implementation is a complete rewrite in C++ and leverages the OpenACC programming model to achieve acceleration on NVIDIA GPUs. While a more comprehensive description of the code and its several other modules will be presented in a future paper, here we focus on some preparatory results that demonstrate the code potential and performance on pre exa-scale parallel architectures.
@article{arxiv.2511.20337,
title = {The PLUTO Code on GPUs: A First Look at Eulerian MHD Methods},
author = {Marco Rossazza and Andrea Mignone and Matteo Bugli and Stefano Truzzi and Lubomir Riha and Tomas Panoc and Ondrej Vysocky and Nitin Shukla and Alessandro Romeo and Vittoria Berta},
journal= {arXiv preprint arXiv:2511.20337},
year = {2025}
}