GPU porting of ASCOT5 code for Monte Carlo simulations in fusion plasmas

Loading...
Thumbnail Image

Access rights

openAccess
CC BY
publishedVersion

URL

Journal Title

Journal ISSN

Volume Title

A1 Alkuperäisartikkeli tieteellisessä aikakauslehdessä

Major/Subject

Mcode

Degree programme

Language

en

Pages

Series

Plasma Physics and Controlled Fusion, Volume 67, issue 10

Abstract

We present the graphics processing unit (GPU) porting of a Monte Carlo particle-following code, called ASCOT5. The code solves the distribution function of minority species in fusion plasmas. Originally developed with an MPI-OpenMP hybrid parallelism taking full advantage of single instruction, multiple data operations, the code has been ported to GPU architecture using the OpenACC programming model. Subsequently, modifications were made to implement three distinct algorithmic strategies: history-based, event-based, and event-based-packing. In the first implementation, each GPU processing unit deals with the entire history of one or more particles, while the event-based algorithm operates on the principle of executing a single low-level event type at a time for all particles still alive (with or without packing particles). Performance results on NVIDIA GPUs are presented to showcase the effectiveness and efficiency of the code adaptations for GPU execution. These results provide insights into the comparative performance of the implemented approaches on the specified hardware architecture. Portability across other architectures such as INTEL and AMD GPUs with OpenMP Offload is also presented.

Description

Publisher Copyright: © 2025 The Author(s). Published by IOP Publishing Ltd.

Other note

Citation

Peybernes, M, Fourestey, G, Äkäslompolo, S, Särkimäki, K & Spiga, F 2025, 'GPU porting of ASCOT5 code for Monte Carlo simulations in fusion plasmas', Plasma Physics and Controlled Fusion, vol. 67, no. 10, 105020. https://doi.org/10.1088/1361-6587/ae0bf3