Hi,<br><br>Did you try numpy.dot(sqrtm(a), sqrtm(a)) ?<br><br>Matthieu<br><br><div><span class="gmail_quote">2007/7/16, Kevin Jacobs <<a href="mailto:jacobs@bioinformed.com">jacobs@bioinformed.com</a>> <<a href="mailto:bioinformed@gmail.com">
bioinformed@gmail.com</a>>:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Hi all,<br><br>This is a bit of a SciPy question, but I thought I'd ask here since I'm already subscribed.  I'd like to add some new LAPACK bindings to SciPy and was wondering if there was a minimum version requirement for LAPACK, since it would be ideal if I could use some of the newer 
3.0 features.  In addition to using some block methods only added in 3.0, it is very convenient to use the WORK=-1 for space queries instead of reimplementing the underlying logic in the calc_work module.<br><br>The routines of most interest to me are:
<br>DGELSD<br>DGGGLM<br>DGGLSE<br><br>I've also found that SciPy's sqrtm is broken:<br><br>>>> a=matrix([[59, 64, 69],<br>        [64, 72, 80],<br>        [69, 80, 91]])<br>>>> sqrtm(a)<br>array([[ 
5.0084801  +1.03420519e-08j,  4.40747825 -2.06841037e-08j,<br>         3.8064764  +1.03420519e-08j],<br>       [ 4.40747825 -2.06841037e-08j,  4.88353492 +4.13682075e-08j,<br>         5.3595916  -2.06841037e-08j],<br>       [ 
3.8064764  +1.03420519e-08j,  5.3595916  -2.06841037e-08j,<br>         6.9127068  +1.03420519e-08j]])<br>>>> sqrtm(a)*sqrtm(a)<br>array([[ 25.08487289 +1.03595922e-07j,  19.42586452 -1.82329475e-07j,<br>         
14.48926259 +7.87335527e-08j],<br>       [ 19.42586452 -1.82329475e-07j,  23.84891336 +4.04046172e-07j,<br>         28.72522212 -2.21716697e-07j],<br>       [ 14.48926259 +7.87335527e-08j,  28.72522212 -2.21716697e-07j,<br>

         47.78551529 +1.42983144e-07j]])<br><br>So not even close...  (and yes, it does deserve a bug report if one doesn't already exist)<br><span class="sg"><br>-Kevin<br><br>
</span><br>_______________________________________________<br>Numpy-discussion mailing list<br><a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:Numpy-discussion@scipy.org">Numpy-discussion@scipy.org
</a><br><a onclick="return top.js.OpenExtLink(window,event,this)" href="http://projects.scipy.org/mailman/listinfo/numpy-discussion" target="_blank">http://projects.scipy.org/mailman/listinfo/numpy-discussion</a><br><br></blockquote>
</div><br>