A Parallel Approach for a Block Tri-Diagonal Solver in Fusion Simulation

摘要

Simulations of fusion experiments provide important support for fusion experiments. Up to date, the study of fusion experiments has been focused on numerical methods by means of computing power. However, the accurate numerical method has always been a major issue. Based on a Scalable Island Equilibrium Solver for Toroidal Application (SIESTA), an algorithm for solving a huge tri-diagonal system is introduced in this paper. The main effort of the paper is placed on studying the computation of a huge tri-diagonal matrix in parallel. A parallel approach is investigated, in which a mapping is defined. The approach guarantees each different matrix block is associated with a unique processor. By means of Basic Linear Algebra Communication Subprograms (BLACS) and Message Passing Interface (MPI), the parallel approach is implemented on the Jaguar supercomputer. The numerical results and efficiency analysis are provided.