A new investigation of the extended Krylov subspace method for matrix function evaluations
Authors:
- 18 July 2010
- journal article
- research article
- Published by Wiley in Numerical Linear Algebra with Applications
- Vol. 17 (4) , 615-638
Abstract
For large square matrices A and functions f, the numerical approximation of the action of f(A) to a vector v has received considerable attention in the last two decades. In this paper we investigate theextended Krylov subspace method, a technique that was recently proposed to approximate f(A)v for A symmetric. We provide a new theoretical analysis of the method, which improves the original result for A symmetric, and gives a new estimate for A nonsymmetric. Numerical experiments confirm that the new error estimates correctly capture the linear asymptotic convergence rate of the approximation. By using recent algorithmic improvements, we also show that the method is computationally competitive with respect to other enhancement techniques. Copyright © 2009 John Wiley & Sons, Ltd.Keywords
This publication has 40 references indexed in Scilit:
- Talbot quadratures and rational approximationsBIT Numerical Mathematics, 2006
- Image numérique, GMRES et polynômes de FaberComptes Rendus Mathematique, 2005
- Computational methods for UV-suppressed fermionsJournal of Computational Physics, 2003
- Numerical methods for the QCDd overlap operator. I. Sign-function and error boundsComputer Physics Communications, 2002
- Magnus-factorized method for numerical solving the time-dependent Schrödinger equationComputer Physics Communications, 2000
- Algorithm 756: a MATLAB toolbox for Schwarz-Christoffel mappingACM Transactions on Mathematical Software, 1996
- Krylov subspace approximation of eigenpairs and matrix functions in exact and computer arithmeticNumerical Linear Algebra with Applications, 1995
- A method for exponential propagation of large systems of stiff nonlinear differential equationsJournal of Scientific Computing, 1989
- Two polynomial methods of calculating functions of symmetric matricesUSSR Computational Mathematics and Mathematical Physics, 1989
- New Approach to Many-State Quantum Dynamics: The Recursive-Residue-Generation MethodPhysical Review Letters, 1983