Abstract
Given a multigrid procedure for linear systems with coefficient matrices \(A_n,\) we discuss the optimality of a related multigrid procedure with the same smoother and the same projector, when applied to properly related algebraic problems with coefficient matrices \(B_n\): we assume that both \(A_n\) and \(B_n\) are Hermitian positive definite with \(A_n\le \vartheta B_n,\) for some positive \(\vartheta \) independent of \(n.\) In this context we prove the Two-Grid Method optimality. We apply this elementary strategy for designing a multigrid solution for modifications of multilevel structured linear systems, in which the Hermitian positive definite coefficient matrix is banded in a multilevel sense. As structured matrices, Toeplitz, circulants, Hartley, sine (\(\tau \) class) and cosine algebras are considered. In such a way, several linear systems arising from the approximation of integro–differential equations with various boundary conditions can be efficiently solved in linear time (with respect to the size of the algebraic problem). Some numerical experiments are presented and discussed, both with respect to Two-Grid and multigrid procedures.
Similar content being viewed by others
References
Aricò, A., Donatelli, M.: A V-cycle multigrid for multilevel matrix algebras: proof of optimality. Numer. Math. 105(4), 511–547 (2007)
Aricò, A., Donatelli, M., Serra-Capizzano, S.: V-cycle optimal convergence for certain (multilevel) structured linear system. SIAM J. Matrix Anal. Appl. 26–1, 186–214 (2004)
Axelsson, O.: Iterative Solution Methods. Cambridge University Press, Cambridge (1994)
Axelsson, O., Neytcheva, M.: The algebraic multilevel iteration methods—theory and applications. In: Bainov D. (ed.) Proceedings of the 2nd International Coll. on Numerical Analysis, Plovdiv (Bulgaria), pp. 13–23 (1993)
Bertaccini, D., Golub, G.H., Serra-Capizzano, S., Tablino-Possio, C.: Preconditioned HSS methods for the solution of non-Hermitian positive definite linear systems and applications to the discrete convection-diffusion equation. Numer. Math. 99, 441–484 (2005)
Bertero, M., Boccacci, P.: Introduction to Inverse Problems in Imaging. Inst. of Physics Publ. Bristol and Philadelphia, London (1998)
Bhatia, R.: Matrix Analysis. Springer, New York (1997)
Bini, D., Capovani, M.: Spectral and computational properties of band symmetric Toeplitz matrices. Linear Algebra Appl. 52(53), 99–125 (1983)
Chan, R.H., Chang, Q., Sun, H.: Multigrid method for ill-conditioned symmetric Toeplitz systems. SIAM J. Sci. Comput. 19–2, 516–529 (1998)
Chan, R.H., Ng, M.K.: Conjugate gradient methods for Toeplitz systems. SIAM Rev. 38, 427–482 (1996)
Chan, R.H., Serra-Capizzano, S., Tablino-Possio, C.: Two-Grid methods for banded linear systems from DCT III algebra. Numer. Linear Algebra Appl. 12–2(3), 241–249 (2005)
Chang, Q., Jin, X., Sun, H.: Convergence of the multigrid method for ill-conditioned block toeplitz systems. BIT 41–1, 179–190 (2001)
Elsner, U., Mehrmann, V., Milde, F., Römer, R.A., Schreiber, M.: The Anderson model of localization: A challenge for modern eigenvalue methods. SIAM J. Sci. Comput. 20, 2089–2102 (1999)
Falgout, R., Vassilevski, P., Zikatanov, L.: On two-grid convergence estimates. Numer. Linear Algebra Appl. 12, 471–494 (2005)
Fiorentino, G., Serra-Capizzano, S.: Multigrid methods for Toeplitz matrices. Calcolo 28, 283–305 (1991)
Fiorentino, G., Serra-Capizzano, S.: Multigrid methods for symmetric positive definite block Toeplitz matrices with nonnegative generating functions. SIAM J. Sci. Comput. 17–4, 1068–1081 (1996)
Hackbusch, W.: Multigrid Methods and Applications. Springer, Berlin (1985)
Ho, M.K., Ng, M.K.: Splitting iterations for circulant-plus-diagonal systems. Numer. Linear Algebra Appl. 12, 779–792 (2005)
Huckle, T., Staudacher, J.: Multigrid preconditioning and Toeplitz matrices. Electr. Trans. Numer. Anal. 13, 81–105 (2002)
Lund, J., Bowers, K.: Sinc Methods for Quadrature and Differential Equations. SIAM Publications, Philadelphia (1992)
Ng, M.K., Chan, R.H., Chan, T., Yip, A.: Cosine transform preconditioners for high resolution image reconstruction. Linear Algebra Appl. 316, 89–104 (2000)
Ng, M.K., Chan, R.H., Tang, W.C.: A fast algorithm for deblurring models with Neumann boundary conditions. SIAM J. Sci. Comput. 21, 851–866 (1999)
Ng, M.K., Serra-Capizzano, S., Tablino-Possio, C.: Numerical behaviour of multigrid methods for symmetric Sinc–Galerkin systems. Numer. Linear Algebra Appl. 12–2(3), 261–269 (2005)
Noutsos, D., Serra-Capizzano, S., Vassalos, P.: Matrix algebra preconditioners for multilevel Toeplitz systems do not insure optimal convergence rate. Theor. Comput. Sci. 315, 557–579 (2004)
Ruge, J., Stüben, K.: Algebraic multigrid. In: McCormick, S. (ed.) Frontiers in Applied Mathematics: Multigrid Methods, pp. 73–130. SIAM, Philadelphis (1987)
Serra-Capizzano, S.: Multi-iterative methods. Comput. Math. Appl. 26–4, 65–87 (1993)
Serra-Capizzano, S.: Convergence analysis of two grid methods for elliptic Toeplitz and PDEs matrix sequences. Numer. Math. 92–3, 433–465 (2002)
Serra-Capizzano, S.: A note on anti-reflective boundary conditions and fast deblurring models. SIAM J. Sci. Comput. 25–3, 1307–1325 (2003)
Serra-Capizzano, S., Tablino-Possio, C.: Multigrid methods for multilevel circulant matrices. SIAM J. Sci. Comput. 26–1, 55–85 (2004)
Serra-Capizzano, S., Tablino-Possio, C.: A note on algebraic multigrid methods for the discrete weighted Laplacian. Comput. Math. Appl. 60–5, 1290–1298 (2010)
Serra-Capizzano, S., Tablino-Possio, C.: A proposal of multigrid methods for Hermitian positive definite linear systems enjoying an order relation, TR-2012-15, Dipartimento di Matematica e Applicazioni, Università di Milano-Bicocca. http://home.matapp.unimib.it/quaderni_di_dipartimento/2012-15
C. Tablino-Possio, V-cycle optimal convergence for DCT-III matrices. In: Bini, D., Mehrmann, V., Olshevsky, V., Tyrtyshnikov, E., Van Barel, M. (eds.) Operator Theory: Advances and Applications., 199 —Numerical Methods for Structured Matrices and Applications: Georg Heinig memorial volume, pp. 377–396. Birkhauser Verlag (2010)
Tyrtyshnikov, E.: Circulant preconditioners with unbounded inverses. Linear Algebra Appl. 216, 1–23 (1995)
Trottenberg, U., Oosterlee, C.W., Schüller, A.: Multigrid. Academic Press, London (2001)
Acknowledgments
The authors thank the Referees for their very pertinent and useful remarks.
Author information
Authors and Affiliations
Corresponding author
Additional information
The work of the authors was partially supported by MIUR 20083KLJEZ.
Rights and permissions
About this article
Cite this article
Serra-Capizzano, S., Tablino-Possio, C. Two-Grid Methods for Hermitian positive definite linear systems connected with an order relation. Calcolo 51, 261–285 (2014). https://doi.org/10.1007/s10092-013-0081-9
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10092-013-0081-9