Exploiting Kronecker structure in exponential integrators: fast approximation of the action of φ-functions of matrices via quadrature

by   Matteo Croci, et al.

In this article, we propose an algorithm for approximating the action of φ-functions of matrices against vectors, which is a key operation in exponential time integrators. In particular, we consider matrices with Kronecker sum structure, which arise from problems admitting a tensor product representation. The method is based on quadrature approximations of the integral form of the φ-functions combined with a scaling and modified squaring method. Owing to the Kronecker sum representation, only actions of 1D matrix exponentials are needed at each quadrature node and assembly of the full matrix can be avoided. Additionally, we derive a priori bounds for the quadrature error, which show that, as expected by classical theory, the rate of convergence of our method is supergeometric. Guided by our analysis, we construct a fast and robust method for estimating the optimal scaling factor and number of quadrature nodes that minimizes the total cost for a prescribed error tolerance. We investigate the performance of our algorithm by solving several linear and semilinear time-dependent problems in 2D and 3D. The results show that our method is accurate and orders of magnitude faster than the current state-of-the-art.


page 1

page 2

page 3

page 4


A μ-mode approach for exponential integrators: actions of φ-functions of Kronecker sums

We present a novel method for computing actions of the so-called φ-funct...

Exponential of tridiagonal Toeplitz matrices: applications and generalization

In this paper, an approximate method is presented for computing exponent...

Randomized sketching for Krylov approximations of large-scale matrix functions

The computation of f(A)b, the action of a matrix function on a vector, i...

Parallel Computation of functions of matrices and their action on vectors

We present a novel class of methods to compute functions of matrices or ...

Direction splitting of φ-functions in exponential integrators for d-dimensional problems in Kronecker form

In this manuscript, we propose an efficient, practical and easy-to-imple...

Efficient and accurate computation to the φ-function and its action on a vector

In this paper, we develop efficient and accurate algorithms for evaluati...

A Gauss Laguerre approach for the resolvent of fractional powers

This paper introduces a very fast method for the computation of the reso...

Please sign up or login with your details

Forgot password? Click here to reset