Implementation of matrix exponential times vector

Hi, thanks for this great software. I would like to request a new feature with regards to the calculation of a matrix exponential times a vector, which can be considerable faster for large matrices than calculating exp(M) and then taking the product. So, I would need exp(M)*v, where M is a square matrix (nxn) and v is a vector (nx1). A template to look at is given e.g. here:

https://github.com/higham/expmv

And some literature: A. H. Al-Mohy and N. J. Higham, "Computing the action of the matrix exponential, with an application to exponential integrators" SIAM J. Sci. Comput., 33(2):488--511, 2011.

I believe, that this would be great feature for the quantum dynamics field.

Thanks!

Alice