Any recommendations for libraries to implement calculation of the matrix exponential?
Expokit is written in fortran, but can be embedded in C++. It works great and includes an optimized algorithm for sparse matrices.