A block-tridiagonal solver with two-level parallelization for finite element-spectral codesopen access
- Authors
- Lee, Jungpyo; Wright, John C.
- Issue Date
- Jun-2014
- Publisher
- ELSEVIER SCIENCE BV
- Keywords
- Two-level parallel computing; Block tridiagonal solver; Partitioned Thomas method; Cyclic odd even reduction; Plasma waves
- Citation
- COMPUTER PHYSICS COMMUNICATIONS, v.185, no.10, pp.2598 - 2608
- Indexed
- SCIE
SCOPUS
- Journal Title
- COMPUTER PHYSICS COMMUNICATIONS
- Volume
- 185
- Number
- 10
- Start Page
- 2598
- End Page
- 2608
- URI
- https://scholarworks.bwise.kr/hanyang/handle/2021.sw.hanyang/159702
- DOI
- 10.1016/j.cpc.2014.06.006
- ISSN
- 0010-4655
- Abstract
- Two-level parallelization is introduced to solve a massive block-tridiagonal matrix system. One-level is used for distributing blocks whose size is as large as the number of block rows due to the spectral basis, and the other level is used for parallelizing in the block row dimension. The purpose of the added parallelization dimension is to retard the saturation of the scaling due to communication overhead and inefficiencies in the single-level parallelization only distributing blocks. As a technique for parallelizing the tridiagonal matrix, the combined method of "Partitioned Thomas method" and "Cyclic Odd Even Reduction" is implemented in an MPI-Fortran90 based finite element-spectral code (TORIC) that calculates the propagation of electromagnetic waves in a tokamak. The two-level parallel solver using thousands of processors shows more than 5 times improved computation speed with the optimized processor grid compared to the single-level parallel solver under the same conditions. Three-dimensional RF field reconstructions in a tokamak are shown as examples of the physics simulations that have been enabled by this algorithmic advance.
- Files in This Item
-
- Appears in
Collections - 서울 공과대학 > 서울 원자력공학과 > 1. Journal Articles
Items in ScholarWorks are protected by copyright, with all rights reserved, unless otherwise indicated.