<div dir="ltr"><br><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Oct 9, 2013 at 8:04 AM, Charanpal Dhanjal <span dir="ltr"><<a href="mailto:charanpal@gmail.com" target="_blank">charanpal@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">I get a segmentation fault upon running the following:<br>
<br>
import numpy<br>
A = numpy.ones((7000000, 8))<br>
Q, R = numpy.linalg.qr(A)<br>
<br>
on Python 2.7.3, Linux 64-bit using numpy 1.9.0.dev-ec3603f linked<br>
against OpenBLAS. If A is a smaller matrix then the QR decomposition<br>
works (for example A has shape (4000000, 8)). I haven't quite narrowed<br>
down the exact threshold where the crash occurs, however I know that the<br>
above A is 448 MB (Q and R are no bigger), and the machine in question<br>
has 32GB of RAM. I also tested scipy.linalg.qr (version<br>
0.14.0.dev-ced994c) with the same results.<br>
<br>
I don't get the same problem on my laptop which is Python 2.7.3, Linux<br>
64-bit but with numpy 1.8.0rc1 linked to OpenBLAS. Both machines have<br>
OpenBLAS 0.2.6. Does anyone have some insight into why this problem is<br>
occuring?<br>
<br></blockquote><div><br></div><div>Works fine here with '1.9.0.dev-7457f15' linked against ATLAS. I suspect problem is in OpenBLAS. What architecture/os do the two machines have? Is OpenBLAS using more than one core?<br>
<br></div><div>Chuck<br></div><br></div></div></div>