Computation-Communication Overlap Techniques for Parallel Spectral Calculations in Gyrokinetic Vlasov Simulations
نویسندگان
چکیده
One of the important phenomena in magnetically-confined fusion plasma is plasma turbulence, which causes particle and heat transport and degrades plasma confinement. To address multi-scale turbulence including temporal and spatial scales of electrons and ions, we extend our gyrokinetic Vlasov simulation code GKV to run efficiently on peta-scale supercomputers. A key numerical technique is the parallel Fast Fourier Transform (FFT) required for parallel spectral calculations, where masking of the cost of inter-node transpose communications is essential to improve strong scaling. To mask communication costs, computation-communication overlap techniques are applied for FFTs and transpose with the help of the hybrid parallelization of message passing interface and open multi-processing. Integrated overlaps including whole spectral calculation procedures show better scaling than simple overlaps of FFTs and transpose. The masking of communication costs significantly improves strong scaling of the GKV code, and makes substantial speed-up toward multi-scale turbulence simulations.
منابع مشابه
Optimizations of the gyrokinetic Vlasov simulation code GKV for the K computer
We report the optimizations of the gyrokinetic Vlasov simulation code GKV, which performs five-dimensional CFD calculation by using spectral and finite difference methods, on for the K computer. First, the segmented rank mapping on the three-dimensional torus network is advantageous to utilize the bi-section bandwidth and to avoid collisions of the messages, and, therefore, to reduce communicat...
متن کاملA Numerical Method for Parallel Particle Motions in Gyrokinetic Vlasov Simulations∗)
A semi-Lagrangian scheme is applied for the first time to computations of charged particle motions along magnetic field lines, to numerically solve the δf gyrokinetic equations in a flux tube geometry. This new solver adopted in the gyrokinetic Vlasov simulations has an advantage over the conventional Eulerian codes in calculating the parallel dynamics, because semi-Lagrangian schemes are free ...
متن کاملOptimizations of a spectral/finite difference gyrokinetic code for improved strong scaling toward million cores
No optimizations Map. Map.+Simul.Comm. Map.+Simul.Comm. +Overlap(spectral) Map.+Simul.Comm. +Overlap(spectral+finite diff.) Problem size: (nx, ny, nz, nv, nμ, ns) = (256, 256, 32, 32, 32, 2) Parallelization: (Nxy, Nz, Nv, Nμ, Ns, Nthreads) = (4, 4, 4, 8, 2, 4) on BX900, (4, 4, 4, 4, 2, 8) on K Optimizations of a spectral/finite difference gyrokinetic code for improved strong scaling toward mill...
متن کاملAn efficient numerical method for gyrokinetic δ f Vlasov simulations in helical plasmas
Three-dimensional configurations of helical plasmas bring numerical difficulties to gyrokinetic simulations of plasma turbulent transport in toroidal fusion devices. Watanabe et al. have investigated the ion-temperature-gradient (ITG)-driven turbulence in helical plasmas by means of the gyrokinetic Vlasov simulation code GVK [1]. Because of a large number of grid points along a magnetic field l...
متن کاملOverlap-based feature weighting: The feature extraction of Hyperspectral remote sensing imagery
Hyperspectral sensors provide a large number of spectral bands. This massive and complex data structure of hyperspectral images presents a challenge to traditional data processing techniques. Therefore, reducing the dimensionality of hyperspectral images without losing important information is a very important issue for the remote sensing community. We propose to use overlap-based feature weigh...
متن کامل