Abstract
The four vector extrapolation methods, minimal polynomial extrapolation, reduced rank extrapolation, modified minimal polynomial extrapolation and the topological epsilon algorithm, when applied to linearly generated vector sequences are Krylov subspace methods and it is known that they are equivalent to some well-known conjugate gradient type methods. However, the vector ε-algorithm is an extrapolation method, older than the four extrapolation methods above, and no similar results are known for it. In this paper, a determinantal formula for the vector ε-algorithm is given. Then it is shown that, when applied to a linearly generated vector sequence, the algorithm is also a Krylov subspace method and for a class of matrices the method is equivalent to a preconditioned Lanczos method. A new determinantal formula for the CGS is given, and an algebraic comparison between the vector ε-algorithm for linear systems and CGS is also given.
Similar content being viewed by others
References
W.E. Arnoldi, The principle of minimized iterations in the solution of the matrix eigenvalue problem, Quart. Appl. Math. 9 (1951) 17–29.
O. Axelsson, Conjugate gradient type methods for unsymmetric and inconsistent systems of linear equations, Linear Algebra Appl. 29 (1980) 1–16.
C. Brezinski, Some results in the theory of the vector e-algorithm, Linear Algebra Appl. 8 (1974) 77–86.
C. Brezinski, Généralisations de la transformation de Shanks, de la table de Padé et de l' e-algorithme, Calcolo 12 (1975) 317–360.
C. Brezinski, Padé-Type Approximation and General Orthogonal Polynomials (Birkhäuser, Basel, 1980).
C. Brezinski and M. Redivo-Zaglia, Treatment of near-breakdown in the CGS algorithms, Numer. Algorithms 7 (1994) 33–73.
C. Brezinski and M. Redivo-Zaglia, Look-ahead in BiCGSTAB and other product-type methods for linear systems, BIT 35 (1995) 169–201.
C. Brezinski and M. Redivo-Zaglia, Extrapolation Methods. Theory and Practice (North-Holland, Amsterdam, 1996).
C. Brezinski and H. Sadok, Avoiding breakdown in the CGS algorithm, Numer. Algorithms 1 (1991) 199–206.
C. Brezinski and H. Sadok, Lanczos-type algorithms for solving systems of linear equations, Appl. Numer. Math. 11 (1993) 443–473.
S. Cabay and L.W. Jackson, A polynomial extrapolation method for finding limits and antilimits of vector sequences, SIAM J. Numer. Anal. 13 (1976) 734–752.
R.P. Eddy, Extrapolating to the limit of a vector sequence, in: Information Linkage between Applied Mathematics and Industry, ed. P.C.C. Wang (Academic Press, New York, 1979) pp. 387–396.
S.C. Eisenstat, H.C. Elman and M.H. Schultz, Variational iterative methods for nonsymmetric systems of linear equations, SIAM J. Numer. Anal. 20 (1983) 345–357.
W. Gander, G.H. Golub and D. Gruntz, Solving linear equations by extrapolation, in: Supercomputing, Trondheim, 1989, Computer Systems Science, Vol. 62 (Springer, Berlin, 1989) pp. 279–293.
P.R. Graves-Morris, Vector valued rational interpolants I, Numer. Math. 42 (1983) 331–348.
P.R. Graves-Morris, Vector valued rational interpolants II, IMA J. Numer. Anal. 4 (1984) 209–224.
P.R. Graves-Morris, G.A. Baker Jr. and C.F. Woodcock, Cayley's theorem and its application in the theory of vector Padé approximants, J. Comput. Appl. Math. 66 (1996) 255–265.
P.R. Graves-Morris and C.D. Jenkins, Vector-valued, rational interpolants III, Constr. Approx. 2 (1986) 263–289.
C. Lanczos, Solution of systems of linear equations by minimized iteration, J. Res. N.B.S. 49 (1952) 33–53.
J.B. McLeod, A note on the e-algorithm, Computing 7 (1971) 17–24.
M. Mešina, Convergence acceleration for the iterative solution of the equations X= AX+ f, Comput. Methods Appl. Mech. Engrg. 10 (1977) 165–173.
J. Nuttall, Convergence of Padé approximants of meromorphic functions, J. Math. Anal. Appl. 31 (1970) 147–153.
D.E. Roberts, Clifford algebras and vector-valued rational forms I, Proc. Roy. Soc. London A 431 (1990) 285–300.
D.E. Roberts, Clifford algebras and vector-valued rational forms II, Numer. Algorithms 3 (1992) 371–381.
D.E. Roberts, On a q-d algorithm, Adv. Comput. Math. 8 (1998) 193–219.
Y. Saad, Krylov subspace methods for solving large unsymmetric linear systems, Math. Comp. 37 (1981) 105–126.
Y. Saad, Iterative Methods for Sparse Linear Systems (PWS, Boston, 1996).
A. Salam, Formal vector orthogonal polynomials, Adv. Comput. Math. 8 (1988) 267–289.
A. Salam, Vector Padé-type approximants and vector Padé approximants, J. Approx. Theory 97 (1999) 92–112.
A. Salam, On vector Hankel determinant, Linear Algebra Appl. 313 (2000) 127–139.
D. Shanks, Non-linear transformations of divergent and slowly convergent sequences, J. Math. Phys. 34 (1955) 1–42.
A. Sidi, Convergence and stability properties of minimal polynomial and reduced rank extrapolation algorithms, SIAM J. Numer. Anal. 23 (1986) 197–209.
A. Sidi, Extrapolation vs. projection methods for linear systems of equations, J. Comput. Appl. Math. 22 (1988) 71–88.
A. Sidi, W.F. Ford and D.A. Smith, Acceleration of convergence of vectors sequences, SIAM J. Numer. Anal. 23 (1986) 178–196.
D.A. Smith, W.F. Ford and A. Sidi, Extrapolation methods for vector sequences, SIAM J. Numer Anal. 23 (1986) 178–196.
P. Sonneveld, CGS, A fast Lanczos-type solver for nonsymmetric linear systems, SIAM J. Sci. Statist. Comput. 10 (1989) 35–52.
H.A. Van Der Vorst, BiCGSTAB: A fast and smoothly converging variant of the Bi-CG for the solution of nonsymmetric linear systems, SIAM J. Sci. Statist. Comput. 13 (1992) 631–644.
P. Wynn, On a device for computing the e m (S n )transformation, Math. Tables Autom. Comp. 10 (1956) 91–96.
P. Wynn, Acceleration techniques for iterated vector and matrix problems, Math. Comput. 16 (1962) 301–322.
P. Wynn, Vector continued fractions, Linear Algebra Appl. 1 (1968) 357–395.
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Salam, A., Graves-Morris, P. On the Vector ε-Algorithm for Solving Linear Systems of Equations. Numerical Algorithms 29, 229–247 (2002). https://doi.org/10.1023/A:1014832627766
Issue Date:
DOI: https://doi.org/10.1023/A:1014832627766