Design and analysis of Toeplitz preconditioners

A general approach is presented for constructing preconditioners P for a symmetric positive-definite Toeplitz matrix A so that the Toeplitz system Ax=b can be solved efficiently by the preconditioned conjugate-gradient (PCG) method. The work is a generalization of recent work by G. Strang (1986) and R. Chan (1988). It is shown that the new preconditioners can be inverted directly by means of various fast transform methods with O(N log N) operations and that the eigenvalues of the preconditioned matrices P/sup -1/ A are clustered around one so that the PCG method converges very quickly.<<ETX>>

[1]  G. Strang,et al.  Toeplitz equations by conjugate gradients with circulant preconditioner , 1989 .

[2]  Hideo Kitajima.  A Symmetric Cosine Transform , 1980, IEEE Transactions on Computers.

[3]  T. Chan.  An Optimal Circulant Preconditioner for Toeplitz Systems , 1988 .

[4]  B. Anderson,et al.  Asymptotically fast solution of toeplitz and related systems of linear equations , 1980 .

[5]  David Y. Y. Yun,et al.  Fast Solution of Toeplitz Systems of Equations and Computation of Padé Approximants , 1980, J. Algorithms.

[6]  Gene H. Golub,et al.  Matrix computations , 1983 .

[7]  J. Makhoul.  On the eigenvectors of symmetric Toeplitz matrices , 1981 .

[8]  Philippe Delsarte,et al.  The split Levinson algorithm , 1986, IEEE Trans. Acoust. Speech Signal Process..

[9]  F. Hoog.  A new algorithm for solving Toeplitz systems of equations , 1987 .

[10]  J. L. Hock,et al.  An exact recursion for the composite nearest‐neighbor degeneracy for a 2×N lattice space , 1984 .

[11]  W. Gragg,et al.  Superfast solution of real positive definite toeplitz systems , 1988 .

[12]  R. H. Chan.  The spectrum of a family of circulant preconditioned Toeplitz systems , 1989 .

[13]  J. Cullum,et al.  Lanczos algorithms for large symmetric eigenvalue computations , 1985 .

[14]  J. H. Wilkinson.  The algebraic eigenvalue problem , 1966 .

[15]  B. Dickinson.  Solution of linear equations with rational Toeplitz matrices , 1980 .

[16]  K. R. Rao,et al.  A Fast Computational Algorithm for the Discrete Sine Transform , 1980, IEEE Trans. Commun..

[17]  B. Dickinson.  Efficient solutions of linear equations with banded Toeplitz matrices , 1979 .

[18]  Henrique S. Malvar.  Fast computation of the discrete cosine transform and the discrete Hartley transform , 1987, IEEE Trans. Acoust. Speech Signal Process..

[19]  A. Cantoni,et al.  Eigenvalues and eigenvectors of symmetric centrosymmetric matrices , 1976 .

[20]  N. Levinson.  The Wiener (Root Mean Square) Error Criterion in Filter Design and Prediction , 1946 .

[21]  David G. Luenberger,et al.  Linear and nonlinear programming , 1984 .

[22]  G. Strang.  A proposal for toeplitz matrix calculations , 1986 .