English

Krylov subspace restarting for matrix Laplace transforms

Numerical Analysis 2023-11-17 v3 Numerical Analysis

Abstract

A common way to approximate F(A)bF(A)b -- the action of a matrix function on a vector -- is to use the Arnoldi approximation. Since a new vector needs to be generated and stored in every iteration, one is often forced to rely on restart algorithms which are either not efficient, not stable or only applicable to restricted classes of functions. We present a new representation of the error of the Arnoldi iterates if the function FF is given as a Laplace transform. Based on this representation we build an efficient and stable restart algorithm. In doing so we extend earlier work for the class of Stieltjes functions which are special Laplace transforms. We report several numerical experiments including comparisons with the restart method for Stieltjes functions.

Keywords

Cite

@article{arxiv.2205.13842,
  title  = {Krylov subspace restarting for matrix Laplace transforms},
  author = {Andreas Frommer and Karsten Kahl and Marcel Schweitzer and Manuel Tsolakis},
  journal= {arXiv preprint arXiv:2205.13842},
  year   = {2023}
}
R2 v1 2026-06-24T11:30:40.294Z