Approximating the matrix exponential of an advection-diffusion operator using the incomplete orthogonalization method
2015 (English)In: Lecture Notes in Computational Science and Engineering, ISSN 1439-7358, 345-353 p.Article in journal (Refereed) Published
In this paper we give first results for the approximation of eAb, i.e. the matrix exponential times a vector, using the incomplete orthogonalization method. The benefits compared to the Arnoldi iteration are clear: shorter orthogonalization lengths make the algorithm faster and a large memory saving is also possible. For the case of three term orthogonalization recursions, simple error bounds are derived using the norm and the field of values of the projected operator. In addition, an a posteriori error estimate is given which in numerical examples is shown to work well for the approximation. In the numerical examples we particularly consider the case where the operator A arises from spatial discretization of an advection-diffusion operator.
Place, publisher, year, edition, pages
2015. 345-353 p.
Advection, Error analysis, A-posteriori error estimates, Advection-diffusion operators, Arnoldi iterations, Field of values, Matrix exponentials, Memory savings, Orthogonalization, Spatial discretizations, Iterative methods
IdentifiersURN: urn:nbn:se:kth:diva-167774DOI: 10.1007/978-3-319-10705-9__34ScopusID: 2-s2.0-84919797963OAI: oai:DiVA.org:kth-167774DiVA: diva2:814277
QC 201505262015-05-262015-05-222015-05-26Bibliographically approved