Skip to Main content Skip to Navigation
Journal articles

Optimization of the gyroaverage operator based on hermite interpolation

Abstract : Gyrokinetic modeling is appropriate for describing Tokamak plasma turbulence, and the gyroaverage operator is a cornerstone of this approach. In a gyrokinetic code, the gyroaveraging scheme needs to be accurate enough to avoid spoiling the data but also requires a low computation cost because it is applied often on the main unknown, the 5D guiding-center distribution function, and on the 3D electric potentials. In the present paper, we improve a gyroaverage scheme based on Hermite interpolation used in the Gysela code. This initial implementation represents a too large fraction of the total execution time. The gyroaverage operator has been reformulated and is now expressed as a matrix-vector product and a cache-friendly algorithm has been setup. Different techniques have been investigated to quicken the computations by more than a factor two. Description of the algorithms is given, together with an analysis of the achieved performance.
Complete list of metadatas

Cited literature [12 references]  Display  Hide  Download

https://hal.inria.fr/hal-01261427
Contributor : Fabien Rozar <>
Submitted on : Friday, February 5, 2016 - 3:48:36 PM
Last modification on : Friday, June 19, 2020 - 9:22:05 AM
Long-term archiving on: : Friday, November 11, 2016 - 4:48:42 PM

Files

gyroaverage-esaim.pdf
Files produced by the author(s)

Identifiers

Citation

Fabien Rozar, Christophe Steiner, Guillaume Latu, Michel Mehrenberger, Virginie Grandgirard, et al.. Optimization of the gyroaverage operator based on hermite interpolation. ESAIM: Proceedings and Surveys, EDP Sciences, 2016, ⟨10.1051/proc/201653012⟩. ⟨hal-01261427⟩

Share

Metrics

Record views

1210

Files downloads

323