GPUにおける核融合プラズマ乱流コードの最適化

Optimization of fusion plasma turbulence code on GPU

朝比 祐一   ; 井戸村 泰宏   ; 伊奈 拓也

Asahi, Yuichi; Idomura, Yasuhiro; Ina, Takuya

高温の核融合プラズマの特徴は無衝突性であり、速度空間に構造を持つ。そのため、核融合プラズマの閉じ込め性能劣化に繋がるプラズマ乱流の解析においては、通常の3次元流体モデルではなく、5次元の運動論的モデルが用いられる。本講演では、差分法に基づく核融合プラズマ乱流コードGT5DのGPU化を通じ、高次元のステンシル計算における最適化手法について論じる。

Because of its collision-less characteristics, fusion plasmas have fine structures in velocity space. Thus, for the analysis of fusion plasma turbulence which degrades plasma confinement, the five dimensional kinetic models are often employed rather than the usual three dimensional fluid model. In this study, we present the optimization of a fusion plasma code, GT5D, which employs a finite difference method. We then discuss the optimization techniques effective for high dimensional stencil based calculations.