Volume 28, pp. 206-222, 2007-2008.

A generalization of the steepest descent method for matrix functions

M. Afanasjew, M. Eiermann, O. G. Ernst, and S. Güttel

Abstract

We consider the special case of the restarted Arnoldi method for approximating the product of a function of a Hermitian matrix with a vector which results when the restart length is set to one. When applied to the solution of a linear system of equations, this approach coincides with the method of steepest descent. We show that the method is equivalent to an interpolation process in which the node sequence has at most two points of accumulation. This knowledge is used to quantify the asymptotic convergence rate.

Full Text (PDF) [670 KB], BibTeX

Key words

Matrix function, Krylov subspace approximation, restarted Krylov subspace method, restarted Arnoldi/Lanczos method, linear system of equations, steepest descent, polynomial interpolation

AMS subject classifications

65F10, 65F99, 65M20

< Back