I&#39;m bumping into the old &quot;Eigenvalues did not converge&quot; error using numpy.linalg.eigh() on several different linux builds of numpy (1.4.1). The matrix is 166x166. I can compute the eigenvalues on a Macintosh build of numpy, and I can confirm that there aren&#39;t degenerate eigenvalues, and that the matrix appears to be negative definite.<br>
<br>I&#39;ve seen this before (though not for several years), and what I normally do is to build lapack with -O0. This trick did not work in the current instance. Does anyone have any tricks to getting eigh to work?<br><br>
Other weird things that I&#39;ve noticed about this case: I can compute the eigenvalues using eigvals and eigvalsh, and can compute the eigenvals/vecs using eig(). The matrix is real symmetric, and I&#39;ve tested that it&#39;s symmetric enough by forcibly symmetrizing it.<br>
<br>Thanks in advance for any help you can offer.<br clear="all"><br>-- <br>Rick Muller<br><a href="mailto:rpmuller@gmail.com">rpmuller@gmail.com</a><br>505-750-7557<br><br>