A threaded SPIKE algorithm for solving general banded systems
A parallel SPIKE algorithm for solving banded diagonally-dominant and non-diagonally dominant linear systems, has been implemented for shared-memory architectures. The numerical experiments demonstrate that the SPIKE-OpenMP solvers offer a highly competitive threaded alternative strategy to the LAPACK BLAS-threaded LU model, for solving banded linear systems on current multicore architectures.
spike algorithm non-diagonally dominant multicore banded system solver diagonally dominant