Accurate conjugate gradient methods for families of shifted systems

Applied Numerical Mathematics

Volumn 49, Issue 1, 2004, Pages 17-37

  Van den Eshof, Jasper ^a   Sleijpen, Gerard L G ^b  

^a HEINRICH HEINE UNIVERSITY   (Germany)

^b UTRECHT UNIVERSITY   (Netherlands)

Author keywords

Accuracy; Finite precision arithmetic; Iterative methods; Shifted systems; Tikhonov regularization

Indexed keywords

ERROR ANALYSIS; LEAST SQUARES APPROXIMATIONS; LINEAR SYSTEMS; MATHEMATICAL TRANSFORMATIONS; MATRIX ALGEBRA; PERTURBATION TECHNIQUES; QUANTUM THEORY;

ACCURACY; FINITE PRECISION ARITHMETIC; SHIFTED SYSTEMS; TIKHONOV REGULARIZATION;

GRADIENT METHODS;

EID: 1542286086     PISSN: 01689274     EISSN: None     Source Type: Journal    
DOI: 10.1016/j.apnum.2003.11.010     Document Type: Conference Paper

References (25)

1. Bai Z., Freund R.W. A symmetric band Lanczos process based on coupled recurrences and some applications. SIAM J. Sci. Comput. 23:(2):2001;542-562. Copper Mountain Conference 2000, MR 1 861 264.
2. Björck Å. Numerical Methods for Least Squares Problems. 1996;SIAM, Philadelphia, PA.
3. Björck Å., Elfving T., Strakoš Z. Stability of conjugate gradient and Lanczos methods for linear least squares problems. SIAM J. Matrix Anal. Appl. 19:(3):1998;720-736.
4. Datta B.N., Saad Y. Arnoldi methods for large Sylvester-like observer matrix equations, and an associated algorithm for partial spectrum assignment. Linear Algebra Appl. 154/156:1991;225-244.
5. I. Dhillon, B. Parlett, Orthogonal eigenvectors and relative gaps, in preparation.
6. Freund R.W. Solution of shifted linear systems by quasi-minimal residual iterations. Reichel L., Ruttan A., Varga R.S. Numerical Linear Algebra, 1992. 1993;101-121 de Gruyter, Berlin.
7. Freund R.W., Nachtigal N.M. An implementation of the QMR method based on coupled two-term recurrences. SIAM J. Sci. Comput. 15:(2):1994;313-337. Iterative methods in numerical linear algebra, Copper Mountain Resort, CO, 1992.
8. A. Frommer, BICGSTAB(l) for families of shifted linear systems, Preprint BUGHW-SC 02/04, Bergische Universität GH Wuppertal, Wuppertal, Germany, November 2002.
9. Frommer A., Glässner U. Restarted GMRES for shifted linear systems. SIAM J. Sci. Comput. 19:(1):1998;15-26. Special issue on iterative methods, Copper Mountain, CO, 1996.
10. Frommer A., Maass P. Fast CG-based methods for Tikhonov-Phillips regularization. SIAM J. Sci. Comput. 20:1999;1831-1850.
11. Frommer A., Nöckel B., Güsken S., Lippert T., Schilling K. Many masses on one stroke: Economic computation of quark. Internat. J. Modern Phys. C. 6:1995;627-638.
12. Golub G.H., Van Loan C.F. Matrix Computations. third ed. 1996;The John Hopkins University Press, Baltimore, MD.
13. Greenbaum A. Behavior of slightly perturbed Lanczos and conjugate-gradient recurrences. Linear Algebra Appl. 113:1989;7-63.
14. Greenbaum A. Estimating the attainable accuracy of recursively computed residual methods. SIAM J. Matrix Anal. Appl. 18:(3):1997;535-551.
15. Gulliksson M.E., Wedin P.-Å., Wei Y. Perturbation identities for regularized Tikhonov inverses and weighted pseudoinverses. BIT. 40:(3):2000;513-523.
16. Hansen P.C. Regularization tools: A Matlab package for analysis and solution of discrete ill-posed problems. Numer. Algorithms. 6:(1-2):1994;1-35.
17. Hestenes M.R., Stiefel E. Methods of conjugate gradients for solving linear systems. J. Res. Nat. Bur. Standards. 49:1952;409-436. (1953).
18. Jegerlehner B. Krylov space solvers for shifted linear systems. http://arXiv.org/abs/hep-lat/9612014:1996;
19. Malyshev A.N. A unified theory of conditioning for linear least squares and Tikhonov regularization solutions. SIAM J. Matrix Anal. Appl. 24:(4):2003;1186-1196.
20. Neuberger H. Overlap Dirac operator. Frommer A., Lippert T., Medeke B., Schilling K. Numerical Challenges in Lattice Quantum Chromodynamics. 2000;Springer, Berlin.
21. Paige C.C. Error analysis of the Lanczos algorithm for tridiagonalizing a symmetric matrix. J. Inst. Math. Appl. 18:(3):1976;341-349.
22. Rutishauser H. Der Quotienten-Differenzen-Algorithmus. Mitt. Inst. Angew. Math. Zürich. 1957:(7):1957;74.
23. Simoncini V. Restarted full orthogonalization method for shifted linear systems. BIT. 43:(2):2003;459-466.
24. 2A+λB+C)x=b and application to structural dynamics. SIAM J. Sci. Comput. 23:(6):2002;1875-1897.
25. Strakoš Z., Tichý P. On error estimation in the conjugate gradient method and why it works in finite precision computations. ETNA. 13:2002;56-80.