[Scipy-tickets] [SciPy] #794: linalg.eig segmentation fault on unpickled arrays

SciPy scipy-tickets@scipy....
Wed Nov 19 07:19:00 CST 2008

#794: linalg.eig segmentation fault on unpickled arrays
 Reporter:  tzito         |       Owner:  somebody
     Type:  defect        |      Status:  new     
 Priority:  high          |   Milestone:  0.7.0   
Component:  scipy.linalg  |     Version:          
 Severity:  critical      |    Keywords:          
 Working on the new linalg.eigh I noticed that feeding linalg.eig (and also
 linalg.eigh) with unpickled arrays leads to a segmentation fault, only if
 ATLAS with SSE2 optimization is installed. This happens on debian lenny,
 with system lapack and atlas libraries and scipy svn r5149.

 To reproduce try the following code:

 import scipy
 import scipy.linalg
 import pickle
 N = 100
 x = scipy.random.random((N, N))
 y = x.copy()
 out = scipy.linalg.eig(y.T, overwrite_a=True)
 print "OK"
 s = pickle.dumps(x)
 yp = pickle.loads(s)
 out = scipy.linalg.eig(yp.T, overwrite_a=True)
 print "OK"

 When libatlas3gf-sse2 is installed I get:

 Segmentation Fault

 When libatlas3gf-sse is installed I get:


 I suspect this problem has something to do with the alignment problems
 described in
 [http://projects.scipy.org/scipy/numpy/ticket/551 numpy ticket #551].

Ticket URL: <http://scipy.org/scipy/scipy/ticket/794>
SciPy <http://www.scipy.org/>
SciPy is open-source software for mathematics, science, and engineering.

More information about the Scipy-tickets mailing list