Work

Highly parallel drift-kinetic semi-Lagrangian simulations in Python

Year: 2024

Type: article

Abstract: \nThis paper describes PyGyro, a parallel 4D drift-kinetic code which implements backward semi-Lagrangian methods combined with field-aligned interpolation and time-splitting. The purpose of PyGyro is twofold: to develop an advanced parallelisation scheme based on fast data transpositions, and to demonstrate the efficacy of a development workflow with code written in Python. Bottlenecks are accelerated with the Pyccel transpiler to achieve Fortran-like execution times. The resulting simulations are shown to correctly reproduce the expected physics and demonstrate good strong scaling, well beyond what was achieved by an existing Fortran code.\n

Source: Springer Link (Chiba Institute of Technology)

DOI: 10.1051/proc/202477176/pdf

Authors: Emily Bourne, Yaman Güçlü

Institutions: Commissariat à l'Énergie Atomique et aux Énergies Alternatives, CEA Cadarache, Institut de Recherche sur la Fusion par confinement Magnétique, École Polytechnique Fédérale de Lausanne, Max Planck Institute for Plasma Physics

Language: English

Topic: Parallel Computing and Optimization Techniques

Subfield: Hardware and Architecture

Field: Computer Science

Domain: Physical Sciences

Open Access status: green

Funders: European Commission, Horizon 2020 Framework Programme