A parallel method for near-circulant penta-diagonal systems

Summary form only given. A fast algorithm for solving a large system with an essentially Toeplitz five-band coefficient matrix is presented. The first two and last two rows are influenced by boundary conditions. The five band core of this matrix is factored as the product of tridiagonal matrices for the purpose of obtaining a solution of a linear system more efficiently. An error term for the approximate solution is presented following the work by Yan and Chung. An algorithm is developed for solving the two systems and is tested using two multiprocessor machines with different architectures.