An L-banded approximation to the inverse of symmetric Toeplitz matrices

We apply the banded matrix inversion theorem given by Kavcic and Moura [IEEE Trans. Inf. Theory 46: 1495-1509, 2000] to symmetric Toeplitz matrices. If the inverse is banded with bandwidth smaller than its size, there is a gain in arithmetic complexity compared to the current methods for Toeplitz matrix inversion. Our algorithm can also be used to find an approximation of the inverse matrix even though it is not exactly banded, but only well localized around its diagonal