A direct Method for the Discrete Solution of Separable Elliptic Equations

This paper extends the direct method of cyclic reduction to linear systems which result from the discretization of separable elliptic equations with Dirichlet, Neumann, or periodic boundary conditions. For an $m \times n$ net, the operation count is proportional to $mn\log _2 n$ and $mn$ storage locations are required.