[Scipy-tickets] [SciPy] #311: The test results in a segmentation fault

SciPy scipy-tickets at scipy.net
Tue Jan 30 14:38:26 CST 2007

#311: The test results in a segmentation fault
 Reporter:  nils            |        Owner:  timl    
     Type:  defect          |       Status:  assigned
 Priority:  normal          |    Milestone:          
Component:  scipy.linsolve  |      Version:  devel   
 Severity:  normal          |   Resolution:          
 Keywords:                  |  
Comment (by mauger):

 Ok...I also get the segfault. The code still has two bugs in it (one of
 which I unsuccessfully tried to fix with the previous patch). The bug
 causing the segfault seen above could be "fixed" (really circumvented) by
 setting the "delta" parameter in get_perm_c() to 0. Otherwise, mmd.c could
 be patched to correctly deal with non-zero delta parameters.

 I've attached a new patch that takes the first approach and sets delta to
 zero. The patch also fixes the second bug. I've tested the patch on the
 following code with no failures.

 from scipy import *
 from pylab import spy, show

 M = sparse.csr_matrix(arange(36).reshape((6,6)))
 lu = linsolve.splu(M)

 n = 20
 A = sparse.lil_matrix((n,n))
 for i in arange(n):
   A[i,:] = random.rand(n)
 A = A.tocsc()
 B = sparse.speye(n,n)
 linsolve.splu(A - B)


 I don't have ARPACK installed so I cannot test the second half of the
 eigs1.py example, but I don't get any errors running the first half.

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

More information about the Scipy-tickets mailing list