T. Jawecki, W. Auzinger, O. Koch:

"Computable upper error bounds for Krylov approximations to matrix exponentials and associated phi-functions";

accepted for publication in BIT Numerical Mathematics (2019), ISBN: 978-3-902627-11-7.

An a posteriori estimate for the error of a standard Krylov approximation

to the matrix exponential is derived.

The estimate is based on the defect (residual) of the Krylov approximation and

is proven to constitute a rigorous upper bound on the error,

in contrast to existing asymptotical approximations.

It can be computed economically in the underlying Krylov space.

In view of time-stepping applications, assuming that

the given matrix is scaled by a time step, it is shown that

the bound is asymptotically correct (with an order related to the

dimension of the Krylov space) for the time step tending to zero.

This means that the deviation of the error estimate from the true error

tends to zero faster than the error itself.

Furthermore, this result is extended to

Krylov approximations of phi-functions

and to improved versions of such approximations.

The accuracy of the derived bounds is demonstrated by examples

and compared with different variants known from the literature,

which are also investigated more closely.

Alternative error bounds are tested on examples, in particular

a version based on the concept of effective order.

For the case where the matrix exponential is used

in time integration algorithms, a step size selectio

matrix exponential · Krylov approximation · a posteriori error estimation · strict upper bound · Schršodinger equation

Created from the Publication Database of the Vienna University of Technology.