[BOK][B] Functions of matrices: theory and computation

NJ Higham - 2008 - SIAM
Functions of matrices have been studied for as long as matrix algebra itself. Indeed, in his
seminal A Memoir on the Theory of Matrices (1858), Cayley investigated the square root of a …

Nineteen dubious ways to compute the exponential of a matrix, twenty-five years later

C Moler, C Van Loan - SIAM review, 2003 - SIAM
In principle, the exponential of a matrix could be computed in many ways. Methods involving
approximation theory, differential equations, the matrix eigenvalues, and the matrix …

On Krylov subspace approximations to the matrix exponential operator

M Hochbruck, C Lubich - SIAM Journal on Numerical Analysis, 1997 - SIAM
Krylov subspace methods for approximating the action of matrix exponentials are analyzed
in this paper. We derive error bounds via a functional calculus of Arnoldi and Lanczos …

Efficient integration of large stiff systems of ODEs with exponential propagation iterative (EPI) methods

M Tokman - Journal of Computational Physics, 2006 - Elsevier
A new class of exponential propagation techniques which we call exponential propagation
iterative (EPI) methods is introduced in this paper. It is demonstrated how for large stiff …

RD-rational approximations of the matrix exponential

I Moret, P Novati - BIT Numerical Mathematics, 2004 - Springer
Restricted Denominator (RD) rational approximations to the matrix exponential operator are
constructed by interpolation in points related to Krylov subspaces associated to a rational …

Error estimates and evaluation of matrix functions via the Faber transform

B Beckermann, L Reichel - SIAM Journal on Numerical Analysis, 2009 - SIAM
The need to evaluate expressions of the form f(A) or f(A)b, where f is a nonlinear function, A
is a large sparse n*n matrix, and b is an n-vector, arises in many applications. This paper …

PARAEXP: A parallel integrator for linear initial-value problems

MJ Gander, S Güttel - SIAM Journal on Scientific Computing, 2013 - SIAM
A novel parallel algorithm for the integration of linear initial-value problems is proposed. This
algorithm is based on the simple observation that homogeneous problems can typically be …

Rational Krylov methods for operator functions

S Güttel - 2010 - eprints.maths.manchester.ac.uk
We present a unified and self-contained treatment of rational Krylov methods for
approximating the product of a function of a linear operator with a vector. With the help of …

Balanced truncation model order reduction in limited time intervals for large systems

P Kürschner - Advances in Computational Mathematics, 2018 - Springer
In this article we investigate model order reduction of large-scale systems using time-limited
balanced truncation, which restricts the well known balanced truncation framework to …

Efficient and stable Arnoldi restarts for matrix functions based on quadrature

A Frommer, S Güttel, M Schweitzer - SIAM Journal on Matrix Analysis and …, 2014 - SIAM
When using the Arnoldi method for approximating f(A)\mathbfb, the action of a matrix
function on a vector, the maximum number of iterations that can be performed is often limited …