A module named poly has been developed for matfaust and pyfaust (it is still in beta status though).
It allows among other functionalities to compute the action of the exponential of a matrix on a vector (using Chebyshev polynomial series).
This module is also GPU compatible.
A Jupyter Notebook is available to introduce how to use this module.
Reference: Sibylle Marcotte, Amélie Barbe, Rémi Gribonval, Titouan Vayer, Marc Sebban, et al.. Fast Multiscale Diffusion on Graphs. 2021. hal-03212764