EEMCS

Home > Publications
Home University of Twente
Education
Research
Prospective Students
Jobs
Publications
Intranet (internal)
 
 Nederlands
 Contact
 Search
 Organisation

EEMCS EPrints Service


18832 Residual, restarting and Richardson iteration for the matrix exponential
Home Policy Brochure Browse Search User Area Contact Help

Botchev, M.A. (2010) Residual, restarting and Richardson iteration for the matrix exponential. Memorandum 1928, Department of Applied Mathematics, University of Twente, Enschede. ISSN 1874-4850

Full text available as:

PDF

321 Kb
PDF (Revised version)

339 Kb
Open Access



Official URL: http://www.math.utwente.nl/publications

Exported to Metis

Abstract

A well-known problem in computing some matrix functions iteratively is a lack of a clear, commonly accepted residual notion. An important matrix function for which this is the case is the matrix exponential. Assume, the matrix exponential of a given matrix times a given vector has to be computed. We interpret the sought after vector as a value of a vector function satisfying the linear system of ordinary differential equations (ODE), whose coefficients form the given matrix. The residual is then defined with respect to the initial-value problem for this ODE system. The residual introduced in this way can be seen as a backward error. We show how the residual can efficiently be computed within several iterative methods for the matrix exponential. This completely resolves the question of reliable stopping criteria for these methods. Furthermore, we show that the residual concept can be used to construct new residual-based iterative methods. In particular, a variant of the Richardson method for the new residual appears to provide an efficient way to restart Krylov subspace methods for evaluating the matrix exponential.

Item Type:Internal Report (Memorandum)
Research Group:EWI-MACS: Mathematics of Computational Science
Research Program:IMPACT-General
Additional Information:Author's family name can also be spelled as "Bochev"
Uncontrolled Keywords:Matrix exponential, Residual, Krylov subspace methods, Restarting, Chebyshev polynomials, Stopping criterion, Richardson iteration, Backward stability, Matrix cosine
ID Code:18832
Deposited On:26 November 2010
More Information:statisticsmetis

Export this item as:

To correct this item please ask your editor

Repository Staff Only: edit this item