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.
Liste complète des métadonnées

Littérature citée [12 références]  Voir  Masquer  Télécharger

Contributeur : Fabien Rozar <>
Soumis le : vendredi 5 février 2016 - 15:48:36
Dernière modification le : jeudi 7 février 2019 - 14:37:33
Document(s) archivé(s) le : vendredi 11 novembre 2016 - 16:48:42


Fichiers produits par l'(les) auteur(s)


  • HAL Id : hal-01261427, version 1
  • ARXIV : 1602.02886


Fabien Rozar, Christophe Steiner, Guillaume Latu, Michel Mehrenberger, Virginie Grandgirard, et al.. Optimization of the gyroaverage operator based on hermite interpolation. Jul 2014, Luminy, France. 2014, 〈http://smai.emath.fr/cemracs/cemracs14/〉. 〈hal-01261427〉



Consultations de la notice


Téléchargements de fichiers