A Cyclic Reduction Algorithm for Solving Block Tridiagonal Systems of Arbitrary Dimension
Roland A. Sweet · SIAM Journal on Numerical Analysis · 1977
A generalization of the Buneman variant of cyclic odd–even reduction algorithm for solving finite difference approximations to Poisson’s equations is presented. This generalization places no restriction on the block size, n, of the system and computes the solution in $O(n^2 \log _2 n)$ operations.