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

[Troels Emtekær Linnet]

Dear Gregor.

Using the scripy method does not sound viable.

It was discussed here,  post 2

http://thread.gmane.org/gmane.science.nmr.relax.devel/6446/match=eliminating+86+bottleneck+numeric+r1rho

The problem is also, that this function expects square matrices.
I need to tweak to data of dimension a, b, c, d, e, X, X where X, X is the
square matrices.
A universal method for the data array.
Best
Troels
On 26 Jun 2014 12:52, "Gregor Thalhammer" <gregor.thalhammer at gmail.com>
wrote:

>
> 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
>
> _______________________________________________
> SciPy-User mailing list
> SciPy-User at scipy.org
> http://mail.scipy.org/mailman/listinfo/scipy-user
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20140626/937ea305/attachment.html>