An asynchronous direct solver for banded linear systems

被引:2
作者
Jandron, Michael A. [1 ]
Ruffa, Anthony A. [1 ]
Baglama, James [2 ]
机构
[1] Naval Undersea Warfare Ctr, Newport, RI 02841 USA
[2] Univ Rhode Isl, Dept Math, Kingston, RI 02881 USA
关键词
LAPACK; Banded; Linear systems; Parallel; Direct solver; Asynchronous; TOEPLITZ; INSTABILITY; ALGORITHM;
D O I
10.1007/s11075-016-0251-3
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
Banded linear systems occur frequently in mathematics and physics. However, direct solvers for large systems cannot be performed in parallel without communication. The aim of this paper is to develop a general asymmetric banded solver with a direct approach that scales across many processors efficiently. The key mechanism behind this is that reduction to a row-echelon form is not required by the solver. The method requires more floating point calculations than a standard solver such as LU decomposition, but by leveraging multiple processors the overall solution time is reduced. We present a solver using a superposition approach that decomposes the original linear system into q subsystems, where q is the number of superdiagonals. These methods show optimal computational cost when q processors are available because each system can be solved in parallel asynchronously. This is followed by a qxq dense constraint matrix problem that is solved before a final vectorized superposition is performed. Reduction to row echelon form is not required by the solver, and hence the method avoids fill-in. The algorithm is first developed for tridiagonal systems followed by an extension to arbitrary banded systems. Accuracy and performance is compared with existing solvers and software is provided in the supplementary material.
引用
收藏
页码:211 / 235
页数:25
相关论文
共 31 条
[1]   Hybrid scheduling for the parallel solution of linear systems [J].
Amestoy, PR ;
Guermouche, A ;
L'Excellent, JY ;
Pralet, S .
PARALLEL COMPUTING, 2006, 32 (02) :136-156
[2]   Multifrontal parallel distributed symmetric and unsymmetric solvers [J].
Amestoy, PR ;
Duff, IS ;
L'Excellent, JY .
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 2000, 184 (2-4) :501-520
[3]   A fully asynchronous multifrontal solver using distributed dynamic scheduling [J].
Amestoy, PR ;
Duff, IS ;
L'Excellent, JY ;
Koster, J .
SIAM JOURNAL ON MATRIX ANALYSIS AND APPLICATIONS, 2001, 23 (01) :15-41
[4]  
Anderson E., 1999, SOC IND APPL MATH
[5]  
[Anonymous], 2013, MATRIX COMPUTATIONS
[6]  
[Anonymous], 1997, Applied numerical linear algebra
[7]  
Bartels R.H., 1998, INTRO SPLINES USE CO, P9
[8]  
Demmel J. W., 1999, SIAM J MATRIX ANAL A, V20
[9]   A survey of recent developments in parallel implementations of Gaussian elimination [J].
Donfack, Simplice ;
Dongarra, Jack ;
Faverge, Mathieu ;
Gates, Mark ;
Kurzak, Jakub ;
Luszczek, Piotr ;
Yamazaki, Ichitaro .
CONCURRENCY AND COMPUTATION-PRACTICE & EXPERIENCE, 2015, 27 (05) :1292-1309
[10]   A review of frontal methods for solving linear systems [J].
Duff, IS .
COMPUTER PHYSICS COMMUNICATIONS, 1996, 97 (1-2) :45-52