<div dir="ltr"><div><div><div>I've indeed opened an issue for this : <a href="https://github.com/numpy/numpy/issues/8090">https://github.com/numpy/numpy/issues/8090</a> . Recently, I've included the LAPACK routines into SciPy dev version that will come with version 0.19. Then you can use ?GECON, ?POCON and other ?XXCON routines for yourself or wait a bit more until I have time to implement it on the SciPy side. <br><br></div>@rkern told me that for NumPy, C translations are involved but I couldn't find an entrance point to contribute for yet. It's a bit above my abilities to fully grasp the way of working in NumPy. You can read more in <a href="https://github.com/numpy/numpy/issues/3755">https://github.com/numpy/numpy/issues/3755</a><br><br></div>Best, <br></div>ilhan<br><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jan 25, 2017 at 9:14 PM, Edward Richards <span dir="ltr"><<a href="mailto:edwardlrichards@gmail.com" target="_blank">edwardlrichards@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">What is the best way to make sure that a matrix inversion makes any sense before preforming it? I am currently struggling to understand some results from matrix inversions in my work, and I would like to see if I am dealing with an ill-conditioned problem. It is probably user error, but I don't like having the possibility hanging over my head.<br>
<br>
I naively put a call to np.linalg.cond into my code; all of my cores went to 100% and a few minutes later I got a number. To be fair A is 6400 elements square, but this takes ~20x more time than the inversion. This is not really practical for what I am doing, is there a better way?<br>
<br>
This is partly in response to Ilhan Polat's post about introducing the A\b operator to numpy. I also couldn't check the Numpy mailing list archives to see if this has been asked before, the numpy-discussion gmane link isn't working for me at all.<br>
<br>
Thanks for your time,<br>
Ned<br>
______________________________<wbr>_________________<br>
NumPy-Discussion mailing list<br>
<a href="mailto:NumPy-Discussion@scipy.org" target="_blank">NumPy-Discussion@scipy.org</a><br>
<a href="https://mail.scipy.org/mailman/listinfo/numpy-discussion" rel="noreferrer" target="_blank">https://mail.scipy.org/mailman<wbr>/listinfo/numpy-discussion</a><br>
</blockquote></div><br></div>