A COMPARISON OF PARALLEL SOLVERS FOR DIAGONALLY DOMINANT
AND GENERAL NARROW BANDED LINEAR SYSTEMS
y z x
PETER ARBENZ ANDREW CLEARY JACK DONGARRA AND MARKUS HEGLAND
Abstract We investigate and compare stable parallel algorithms for solving diagonally dominant and general
narrow banded linear systems of equations Narrow banded means that the bandwidth is very small compared with
the matrix order and is typically b etween and The solvers compared are the banded system solvers of ScaLA
PACK and those investigated by Arb enz and Hegland For the diagonally dominant case the algorithms
are analogs of the well known tridiagonal cyclic reduction algorithm while the inspiration for the general case is
the lesser known bidiagonal cyclic reduction which allows a clean parallel implementation of partial pivoting These
divide and conquer type algorithms complement ne grained algorithms which p erform well only for wide banded ma
trices with each family of algorithms having a range of problem sizes for which it is sup erior We present theoretical
analyses as well as numerical exp eriments conducted on the Intel Paragon
Key words narrow banded linear systems stable factorization parallel solution cyclic reduction ScaLAPACK