Affiliation:
1. Belarusian State University
Abstract
The topic of this paper refers to efficient parallel solvers of block-tridiagonal linear systems of equations. Such systems occur in numerous modeling problems and require usage of high-performance multicore computation systems. One of the widely used methods for solving block-tridiagonal linear systems in parallel is the original block-tridiagonal sweep method. We consider the algorithm based on the partitioning idea. Firstly, the initial matrix is split into parts and transformations are applied to each part independently to obtain equations of a reduced block-tridiagonal system. Secondly, the reduced system is solved sequentially using the classic Thomas algorithm. Finally, all the parts are solved in parallel using the solutions of a reduced system. We propose a modification of this method. It was justified that if known stability conditions for the matrix sweep method are satisfied, then the proposed modification is stable as well.
Publisher
Publishing House Belorusskaya Nauka
Subject
Computational Theory and Mathematics,General Physics and Astronomy,General Mathematics
Reference10 articles.
1. Samarskii A. A., Nikolaev E. S. Numerical Methods for Grid Equations. Vol. 1. Direct Methods. Birkhauser Verlag, 1989. 242 p.
2. Heller D. Some aspects of the cyclic reduction algorithm for block tridiagonal linear systems. SIAM Journal on Numerical Analysis. 1976, vol. 13, no. 4, pp. 484–496. https://doi.org/10.1137/0713042
3. Akimova E. N. Parallelization of the matrix sweep algorithm. Matematicheskoe modelirovanie = Mathematical Models, 1994, vol. 6, no. 9,pp. 61–67 (in Russian).
4. Akimova E. N, Belousov D. V. Parallel algorithms for solving the systems of equations with block-three-diagonal matrices on multiprocessors computer systems. Vestnik UGATU[Scientific Journal of Ufa State Aviation Technical University], 2011 vol. 15, no. 5,pp. 87–93 (in Russian).
5. Hirshman S. P., Perumalla K. S., Lynch V. E., Sanchez R. BCYCLIC: A parallel block tridiagonal matrix cyclic solver. Journal of Computational Physics, 2010, vol. 229, no. 18, pp. 6392–6404. https://doi.org/10.1016/j.jcp.2010.04.049