I provide an improvement on the state of the art technique for solving tridiagonal Toeplitz systems, such as as the Poisson equation.
My method outperforms the Thomas algorithm, LAPACK, and BLAS for this task, in terms of speed, accuracy, and stability.