Comparison of software for computing the action of the matrix exponential

The implementation of exponential integrators requires the action of the matrix exponential and related functions of a large matrix. There are various methods in the literature for carrying out this task. In this paper we describe a new implementation of a method based on interpolation at Leja points. We numerically compare this method with others from the literature. As we are interested in exponential intergrators we choose the test examples from spatial discretization of time dependent partial differential equations in two and three space dimensions. The test matrices thus have large eigenvalues and can be nonnormal.