Optimization of the gyroaverage operator based on Hermite interpolation
From MaRDI portal
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.
Recommendations
- Targeting realistic geometry in tokamak code Gysela
- Nearest-grid-point interpolation in gyrokinetic particle-in-cell simulation
- Towards optimal explicit time-stepping schemes for the gyrokinetic equations
- Development of semi-Lagrangian gyrokinetic code for full-f turbulence simulation in general tokamak geometry
- Field-aligned interpolation for semi-Lagrangian gyrokinetic simulations
Cites work
Cited in
(2)
This page was built for publication: Optimization of the gyroaverage operator based on Hermite interpolation
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q3181272)
