A parallel algorithm for solving block tridiagonal linear systems (Q1205899)

Iterative method: Decoupling of the system into two sets of totally \(q+1\) systems of order \(2p\) (\(n=p\cdot q\)) that can be solved independently on \(q+1\) processors, then mean value of the solution of the two sets (arithmetic mean method). Presentation of theory for convergence, proposition to use the method as preconditioner for the conjugate gradient method, examples on a 4 processor CRAY Y-MP (use of Gauss algorithm and cyclic reduction for individual systems).