[SciPy-dev] Segmentation fault using linsolve

[Nils Wagner]

Hi all,

Can someone reproduce the segfault by running the following test

from scipy import *

n = 15
A = sparse.lil_matrix((n,n))
for i in arange(0,n):
  A[i,:n] = random.rand(n)
B = 2.*sparse.speye(n,n)
sigma =  1.0
sigma_solve = linsolve.splu(A - sigma*B).solve


Nils

I am using the latest svn versions of numpy/scipy.
This is the output of gdb

GNU gdb 6.3
Copyright 2004 Free Software Foundation, Inc.
GDB is free software, covered by the GNU General Public License, and you are
welcome to change it and/or distribute copies of it under certain
conditions.
Type "show copying" to see the conditions.
There is absolutely no warranty for GDB.  Type "show warranty" for details.
This GDB was configured as "x86_64-suse-linux"...(no debugging symbols
found)
Using host libthread_db library "/lib64/tls/libthread_db.so.1".

(gdb) run test_linsolve.py
Starting program: /usr/bin/python test_linsolve.py
(no debugging symbols found)
(no debugging symbols found)
[Thread debugging using libthread_db enabled]
[New Thread 46912509653888 (LWP 6613)]
(no debugging symbols found)
(no debugging symbols found)
(no debugging symbols found)
(no debugging symbols found)
(no debugging symbols found)
(no debugging symbols found)
(no debugging symbols found)
(no debugging symbols found)
(no debugging symbols found)
Use minimum degree ordering on A'+A.

Program received signal SIGSEGV, Segmentation fault.
[Switching to Thread 46912509653888 (LWP 6613)]
0x00002aaab0e6e9a1 in genmmd_ (neqns=0x7fffffa847ac, xadj=0x8d6d60,
adjncy=0x998550, invp=0x8e19ec, perm=0x88cb20, delta=0x7fffffa847a4,
    dhead=0x8e196c, qsize=0x8e6040, llist=0x92bcbc, marker=0x94c78c,
maxint=0x7fffffa847a0, nofsub=0x7fffffa8479c) at mmd.c:162
162             perm[nextmd] = -mdeg;