G. Heinig et F. Hellinger, MOORE-PENROSE INVERSION OF SQUARE TOEPLITZ MATRICES, SIAM journal on matrix analysis and applications, 15(2), 1994, pp. 418-450
Fast algorithms for the computation of the Moore-Penrose inverse A+ of
a square Toeplitz matrix A are constructed based on Bezoutian represe
ntations of A+. Two approaches are presented. The first approach is a
recursion of the nested submatrices, the second approach uses generali
zed inverses of A that are Bezoutians and a global formula.