<br><br><div class="gmail_quote">On Sat, Jul 21, 2012 at 7:21 PM, John Hassler <span dir="ltr">&lt;<a href="mailto:hasslerjc@comcast.net" target="_blank">hasslerjc@comcast.net</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div class="im"><br>
On 7/20/2012 5:00 PM, Pauli Virtanen wrote:<br>
&gt; 20.07.2012 22:24, John Hassler kirjoitti:<br>
&gt;&gt; The offending line is:<br>
&gt;&gt; AA,BB,Q,Z,sdim = qz(A,B,sort=sort)<br>
&gt;&gt;<br>
&gt;&gt; I couldn&#39;t find &#39;sort&#39; defined anywhere, so I arbitrarily changed it to<br>
&gt;&gt; sort=&#39;lhp&#39;.  Then it runs, although the test fails.<br>
&gt;&gt;<br>
&gt;&gt; Is there something else I can try?<br>
&gt; Seems to be a problem with the callable sort function then.<br>
&gt;<br>
&gt; That it works with sort=&#39;lhp&#39; is strange, and probably means that there<br>
&gt; is not a problem with all callback functions, but something goes wrong<br>
&gt; in the algorithm. (Which would be expected, if *gees callbacks work<br>
&gt; without problems.)<br>
&gt;<br>
&gt; If you can/know how to recompile, try recompiling with<br>
&gt;<br>
&gt; set OPT=-g -DDEBUGCFUNCS<br>
&gt; python setup.py .........<br>
&gt;<br>
&gt; This should make the f2py wrappers spit out extra information on what is<br>
&gt; going on.<br>
&gt;<br>
</div>I&#39;m not set up to recompile.  </blockquote><div><br>Is there anyone who can reproduce the issue and investigate this further? Alternatively I can build some binaries with the above debug flags. I can&#39;t reproduce the issue though, so that may be a painful way to debug.<br>
<br>Ralf<br><br> <br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">I set up a little test program to<br>
reproduce the problem (below).  It leads to the Fortran call to gges in<br>
_decomp_qz.py.  This crashes if sort_t = 1 (which calls the lambda), but<br>
not for sort_t = 0.  I didn&#39;t play with it any further.<br>
john<br>
<br>
# Test qz<br>
import numpy as np<br>
from scipy.linalg import qz<br>
<br>
A =   np.array([[3.9,  12.5, -34.5,  -0.5],<br>
                 [ 4.3,  21.5, -47.5,   7.5],<br>
                 [ 4.3,  21.5, -43.5,   3.5],<br>
                 [ 4.4,  26.0, -46.0,   6.0 ]])<br>
<br>
B = np.array([[ 1.0,   2.0,  -3.0,   1.0],<br>
               [1.0,   3.0,  -5.0,   4.0],<br>
               [1.0,   3.0,  -4.0,   3.0],<br>
               [1.0,   3.0,  -4.0,   4.0]])<br>
<br>
<br>
#sort = lambda ar,ai,beta : ai == 0  ## crashes<br>
sort = None                         ## runs<br>
<br>
print(qz(A,B,sort=sort))<br>
<div class="HOEnZb"><div class="h5"><br>
_______________________________________________<br>
SciPy-User mailing list<br>
<a href="mailto:SciPy-User@scipy.org">SciPy-User@scipy.org</a><br>
<a href="http://mail.scipy.org/mailman/listinfo/scipy-user" target="_blank">http://mail.scipy.org/mailman/listinfo/scipy-user</a><br>
</div></div></blockquote></div><br>