Trying to improve performance on Matrix Exponentials in Python
Julia package for adaptive matrix exponential based on tolerance