I narrowed down the segfault from the scipy tests on my machine to: OPENBLAS_CORETYPE='Barcelona' /c/Python34_x64/python -c"import numpy as np; print(np.linalg.svd(np.ones((129, 129), dtype=np.float64))"
Barcelona is the architecture detected by OpenBLAS. If I force Nehalem or if I reduce the matrix to shape (128, 128), it no longer segfaults. To debug it further I think I have to write a test program in C that uses the sgesdd function of the libopenblaspy.dll or against the official libopenblas.dll. Setting the dev environment under windows is painful though. I don't know if it's a bug in OpenBLAS' coretype detection or in the SGESDD kernel for the Barcelona architecture. -- Olivier _______________________________________________ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion
