On 7/20/07, <b class="gmail_sendername">Nils Wagner</b> <<a href="mailto:nwagner@iam.uni-stuttgart.de">nwagner@iam.uni-stuttgart.de</a>> wrote:<div><span class="gmail_quote"></span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Your sqrtm_eig(x) function won't work if x is defective.<br>See test_defective.py for details.<br></blockquote></div><br><br>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.
<br><br>-Kevin<br><br>