Segmantation fault sparse.lu_factor(A).solve(r) (Trac #13)

[scipy-gitbot]

Original ticket http://projects.scipy.org/scipy/ticket/13 on 2006-03-02 by @nilswagner01, assigned to unknown.

from scipy import *

from scipy.sparse import *

n = 20

A = csc_matrix((n,n))

x = rand(n)

y = rand(n-1)+1j*rand(n-1)

r = rand(n)

for i in range(len(x)):

A[i,i] = x[i]

for i in range(len(y)):

A[i,i+1] = y[i]

A[i+1,i] = conjugate(y[i])

xx = sparse.lu_factor(A).solve(r)

Program received signal SIGSEGV, Segmentation fault.
[Switching to Thread 16384 (LWP 15138)]
0x00002aaab2d8f8a5 in dpivotL (jcol=0, u=1, usepr=0x7fffffffd3e4, perm_r=0x6cf550, iperm_r=,
iperm_c=, pivrow=0x7fffffffd3e8, Glu=0x2aaab34bb360, stat=0xffffffffab7dfa48) at dpivotL.c:120
120 perm_r[*pivrow] = jcol;
Current language: auto; currently c
(gdb) bt
#0 0x00002aaab2d8f8a5 in dpivotL (jcol=0, u=1, usepr=0x7fffffffd3e4, perm_r=0x6cf550, iperm_r=,

iperm_c=<value optimized out>, pivrow=0x7fffffffd3e8, Glu=0x2aaab34bb360, stat=0xffffffffab7dfa48) at dpivotL.c:120

#1 0x00002aaab2d85c49 in dgstrf (options=, A=0x7fffffffd570, drop_tol=,

relax=<value optimized out>, panel_size=10, etree=<value optimized out>, work=<value optimized out>, lwork=7140688,
perm_c=0xa17510, perm_r=0x6cf550, L=0x2aaab55658c0, U=0x2aaab55658e0, stat=0x7fffffffd510, info=0x7fffffffd508)
at dgstrf.c:310

#2 0x00002aaab2d6886d in newSciPyLUObject (A=0x7fffffffd670, diag_pivot_thresh=1, drop_tol=0, relax=1, panel_size=10,

permc_spec=2, intype=<value optimized out>) at _superluobject.c:372

#3 0x00002aaab2d67bac in Py_dgstrf (self=, args=,

keywds=<value optimized out>) at _dsuperlumodule.c:187

#4 0x00002aaaaac1c17e in PyCFunction_Call () from /usr/lib64/libpython2.4.so.1.0
#5 0x00002aaaaac4f661 in PyEval_EvalFrame () from /usr/lib64/libpython2.4.so.1.0
#6 0x00002aaaaac51705 in PyEval_EvalCodeEx () from /usr/lib64/libpython2.4.so.1.0
#7 0x00002aaaaac4f789 in PyEval_EvalFrame () from /usr/lib64/libpython2.4.so.1.0
#8 0x00002aaaaac51705 in PyEval_EvalCodeEx () from /usr/lib64/libpython2.4.so.1.0
#9 0x00002aaaaac51992 in PyEval_EvalCode () from /usr/lib64/libpython2.4.so.1.0
#10 0x00002aaaaac6aeb9 in run_node () from /usr/lib64/libpython2.4.so.1.0
#11 0x00002aaaaac6c3c9 in PyRun_SimpleFileExFlags () from /usr/lib64/libpython2.4.so.1.0
#12 0x00002aaaaac6c921 in PyRun_AnyFileExFlags () from /usr/lib64/libpython2.4.so.1.0
#13 0x00002aaaaac72023 in Py_Main () from /usr/lib64/libpython2.4.so.1.0
#14 0x00000000004008d9 in main (argc=, argv=) at ccpython.cc:10