[Scipy-tickets] [SciPy] #1459: Crash in Arpack

SciPy Trac scipy-tickets@scipy....
Mon Jun 13 19:53:41 CDT 2011

#1459: Crash in Arpack
 Reporter:  pv                   |       Owner:  wnbell     
     Type:  defect               |      Status:  new        
 Priority:  high                 |   Milestone:  Unscheduled
Component:  scipy.sparse.linalg  |     Version:  0.9.0      
 Keywords:  arpack segfault      |  
 The following code (courtesy of jakevdp) triggers out-of-bounds memory
 access inside ARPACK:
 import numpy as np
 from scipy.sparse.linalg import eigs

 N = 6
 k = 2

 # with this random seed, I get a memory error on the third
 # iteration below
 A = np.random.random((N,N))

 for i in range(3):
     print i
     evals2, evecs2 = eigs(A, k)
 The memory corruption does not always seem to lead directly to a SEGV, but
 I've seen segfaults in other cases, with probably the same cause.

 Valgrind log is attached. It seems probable that either the work array
 sizes in current Git master (tested with e23e4621) are invalid, or there
 is a bug in ARPACK itself.

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

More information about the Scipy-tickets mailing list