<p dir="auto">I tried to solve the eigenvalues and eigenvectors of a Hermitian matrix, the code used to work nicely and still works in Cocalc, but not on my windows 10 (the system got some updates) now, and the error is very confusing, I can not figure out what is happening now, maybe someone also has a similar issue.</p>
<p dir="auto">Python version: 3.7.2<br>
NumPy version: 1.16<br>
windows version:  Windows Pro 19041.508<br>
Intel: i7-3930K, 64 GB memory.</p>
<h3 dir="auto">Reproducing code example:</h3>
<p dir="auto">I am not sure if this can be reproduced on another system, but it can be reproduced in my computer, I recorded a video showing how that happens:</p>
<p dir="auto"><a href="https://www.youtube.com/watch?v=ZA_0kyXiXco" rel="nofollow">https://www.youtube.com/watch?v=ZA_0kyXiXco</a></p>
<p dir="auto">The matrix is 2048 by 2048, and if I call the   <code class="notranslate">linalg.eigh  </code>twice, it works. That is funny!</p>
<p dir="auto">Add support for the <a href="http://en.wikipedia.org/wiki/Cholesky#LDL_decomposition_2" rel="nofollow">LDL decomposition</a>, which is a variant of the Cholesky decomposition that doesn't take any square roots (faster). LAPACK has a function for this called DPTTRF, so I'm guessing that supporting this is just a matter of adding a wrapper for this function. Also, if this is something that would fit better in Scipy, let me know.</p>