Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

An asynchronous direct solver for banded linear systems

An asynchronous direct solver for banded linear systems 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 q×q 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. http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Numerical Algorithms Springer Journals

An asynchronous direct solver for banded linear systems

Loading next page...
1
 
/lp/springer_journal/an-asynchronous-direct-solver-for-banded-linear-systems-uW2OOih5CU

References (14)

Publisher
Springer Journals
Copyright
Copyright © 2017 by Springer Science+Business Media New York (outside the USA)
Subject
Computer Science; Numeric Computing; Algorithms; Algebra; Theory of Computation; Numerical Analysis
ISSN
1017-1398
eISSN
1572-9265
DOI
10.1007/s11075-016-0251-3
Publisher site
See Article on Publisher Site

Abstract

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 q×q 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.

Journal

Numerical AlgorithmsSpringer Journals

Published: Jan 4, 2017

There are no references for this article.