[Numpy-discussion] expm

Kevin Jacobs <jacobs@bioinformed.com> bioinformed at gmail.com
Fri Jul 20 15:12:28 EDT 2007

On 7/20/07, Nils Wagner <nwagner at iam.uni-stuttgart.de> wrote:
> Your sqrtm_eig(x) function won't work if x is defective.
> See test_defective.py for details.

I've added several defective matrices to my test cases and the SVD method
doesn't work well as I'd thought (which is obvious in retrospect).  I'm
going to circle back and see what I can do to speed up Nick Higham's
algorithm to the point where it is useful for me.  Otherwise, my need is for
a fast sqrtm for symmetric positive definite inputs, so I may still end up
using the SVD approach and contributing a complementary sqrtm_nondefective
back to scipy.

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/numpy-discussion/attachments/20070720/5ba805cd/attachment.html>

More information about the NumPy-Discussion mailing list