Extended Krylov Subspaces: Approximation of the Matrix Square Root and Related Functions
- 1 July 1998
- journal article
- Published by Society for Industrial & Applied Mathematics (SIAM) in SIAM Journal on Matrix Analysis and Applications
- Vol. 19 (3), 755-771
- https://doi.org/10.1137/s0895479895292400
Abstract
We introduce an economical Gram--Schmidt orthogonalization on the extended Krylov subspace originated by actions of a symmetric matrix and its inverse. An error bound for a family of problems arising from the elliptic method of lines is derived. The bound shows that, for the same approximation quality, the diagonal variant of the extended subspaces requires about the square root of the dimension of the standard Krylov subspaces using only positive or negative matrix powers. An example of an application to the solution of a 2.5-D elliptic problem attests to the computational efficiency of the method for large-scale problems. We introduce an economical Gram--Schmidt orthogonalization on the extended Krylov subspace originated by actions of a symmetric matrix and its inverse. An error bound for a family of problems arising from the elliptic method of lines is derived. The bound shows that, for the same approximation quality, the diagonal variant of the extended subspaces requires about the square root of the dimension of the standard Krylov subspaces using only positive or negative matrix powers. An example of an application to the solution of a 2.5-D elliptic problem attests to the computational efficiency of the method for large-scale problems.Keywords
This publication has 17 references indexed in Scilit:
- Using Nonorthogonal Lanczos Vectors in the Computation of Matrix FunctionsSIAM Journal on Scientific Computing, 1998
- Solution of 2.5-dimensional problems using the Lanczos decompositionRadio Science, 1994
- Semivectorial Helmholtz beam propagation by Lanczos reductionIEEE Journal of Quantum Electronics, 1993
- Efficient Solution of Parabolic Equations by Krylov Approximation MethodsSIAM Journal on Scientific and Statistical Computing, 1992
- A comparison of Lanczos electric field propagation methodsJournal of Lightwave Technology, 1992
- Rational Approximations and OrthogonalityTranslations of Mathematical Monographs, 1991
- A method for exponential propagation of large systems of stiff nonlinear differential equationsJournal of Scientific Computing, 1989
- Behavior of slightly perturbed Lanczos and conjugate-gradient recurrencesLinear Algebra and its Applications, 1989
- New Approach to Many-State Quantum Dynamics: The Recursive-Residue-Generation MethodPhysical Review Letters, 1983
- On George’s Nested Dissection MethodSIAM Journal on Numerical Analysis, 1976