The method of odd/even reduction and factorization with application to Poisson''s equation, part II

Several algorithms are presented for solving block tridiagonal systems of linear algebraic equations when the matrices on the diagonal are equal to each other and the matrices on the subdiagonals are all equal to each other. It is shown that these matrices arise from the finite difference approximation to certain elliptic partial differential equations on rectangular regions. Generalizations are derived for higher order equations and non-rectangular regions.