[SciPy-User] Faster approach to calculate the matrix exponential

[Gregor Thalhammer]

Am 26.06.2014 um 10:31 schrieb Troels Emtekær Linnet <tlinnet at gmail.com>:

> Dear NMR wizards.
> 
> Do any of you know any faster approach to calculate the matrix exponential ?
> 
> I currently do it via eigenvalue decomposition approach.
> But my profiling scripts tells me, that my bottleneck is the eig function.
> Stealing 86% of the time.
> 

scipy.linalg.expm provides a matrix exponential using Padé approximation. Is this sufficiently accurate for your application?

Gregor