In this paper we present an algorithm for the construction of the superoptimal circulant preconditioner for a two-level Toeplitz linear system. The algorithm is fast, in the sense that it operates in FFT time. Numerical results are given to assess its performance when applied to the solution of two-level Toeplitz systems by the conjugate gradient method, compared with the Strang and optimal circulant preconditioners.
Similar content being viewed by others
References
O. Axelsson and G. Lindskog, On the rate of convergence of the preconditioned conjugate gradient method, Numer. Math. 48(5) (1986) 499–523.
M.R. Banham and A.K. Katsaggelos, Digital image restoration, IEEE Signal Process. Mag. 14 (1997) 24–41.
R.H. Chan, The spectrum of a family of circulant preconditioned Toeplitz systems, SIAM J. Numer. Anal. 26 (1989) 503–506.
R.H. Chan, X.Q. Jin and M.C. Yeung, The circulant operator in the Banach algebra of matrices, Linear Algebra Appl. 149 (1991) 41–53.
R.H. Chan, X.Q. Jin and M.C. Yeung, The spectra of super-optimal circulant preconditioned Toeplitz systems, SIAM J. Numer. Anal. 28(3) (1991) 871–879.
R.H. Chan and M.K. Ng, Conjugate gradient methods for Toeplitz systems, SIAM Rev. 38(3) (September 1996) 427–482.
R.H. Chan and G. Strang, Toeplitz equations by conjugate gradients with circulant preconditioner, SIAM J. Sci. Statist. Comput. 10(1) (1989) 104–119.
R.H. Chan and M.C. Yeung, Circulant preconditioners for Toeplitz matrices with positive continuous generating functions, Math. Comput. 58(197) (1992) 233–240.
T.F. Chan, An optimal circulant preconditioner for Toeplitz systems, SIAM J. Sci. Statist. Comput. 9(4) (1988) 766–771.
P.J. Davis, Circulant Matrices, Wiley, New York, 1979.
F. Di Benedetto, C. Estatico and S. Serra Capizzano, Superoptimal preconditioned conjugate gradient iteration for image deblurring, SIAM J. Sci. Comput. 26(3) (2005) 1012–1035.
M. Donatelli, C. Estatico and S. Serra Capizzano, Regularization of image restoration problems with anti-reflective boundary conditions, Submitted, (2004).
P.C. Hansen, Deconvolution and regularization with Toeplitz matrices, Numer. Algorithms 29(4) (2002) 323–378.
M. Kac, W.L. Murdock and G. Szegö, On the eigenvalues of certain Hermitian forms, J. Rat. Mech. Anal. 2 (1953) 767–800.
The MathWorks, Inc., Natick, Massachusetts. Matlab ver. 7.0, 2004.
C.V.M. van der Mee, G. Rodriguez and S. Seatzu, Fast superoptimal preconditioning of multiindex Toeplitz matrices. Submitted, (2005).
C.V.M. van der Mee and S. Seatzu, A method for generating infinite positive self-adjoint test matrices and Riesz bases, SIAM J. Matrix Anal. Appl. 26(4) (2005) 1132–1149.
S. Serra Capizzano and E.E. Tyrtyshnikov, Any circulant-like preconditioner for multilevel matrices is not superlinear, SIAM J. Matrix Anal. Appl. 21(2) (1999) 431–439.
A. van der Sluis and H.A. van der Vorst, The rate of convergence of conjugate gradient. Numer. Math. 48(5) (1986) 543–560.
Software in the Public Interest, Inc. Debian/GNU Linux 3.1 (Sarge), 2005. http://www.debian.org/.
G. Strang, A proposal for Toeplitz matrix calculations, Stud. Appl. Math. 74 (1986) 171–176.
V.V. Strela and E.E. Tyrtyshnikov, Which circulant preconditioner is better? Math. Comput. 65(213) (1996) 137–150.
M. Tismenetsky, A decomposition of Toeplitz matrices and optimal circulant preconditioning, Linear Algebra Appl. 154/156 (1991) 105–121.
W.F. Trench, Numerical solution of the eigenvalue problem for Hermitian Toeplitz matrices, SIAM J. Matrix Anal. Appl. 10 (1989) 135–146.
E.E. Tyrtyshnikov, Optimal and superoptimal circulant preconditioners, SIAM J. Matrix Anal. Appl. 13(2) (1992) 459–473.
Author information
Authors and Affiliations
Corresponding author
Additional information
Communicated by C. Brezinski
Rights and permissions
About this article
Cite this article
van der Mee, C., Rodriguez, G. & Seatzu, S. Fast computation of two-level circulant preconditioners. Numer Algor 41, 275–295 (2006). https://doi.org/10.1007/s11075-005-9011-5
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11075-005-9011-5