What is in this document? Mathematical and algorithmic explanation of this process. I will put a naive Python implementation of this algorithm to make it more clear.

Abstract: A Toeplitz-block-Toeplitz (TBT) matrix is block Toeplitz with Toeplitz blocks. TBT systems of equations arise in 2D interpolation, 2D linear prediction and 2D least-squares deconvolution ...

Abstract: Using the structure of the inverse of a Toeplitz matrix and some simplifications resulting when the Toeplitz matrix is banded, algorithms for the solution of a set of n linear equations ...