[Scipy-tickets] [SciPy] #553: sparse solvers segfault with singular matrices

SciPy scipy-tickets@scipy....
Tue Feb 19 14:06:05 CST 2008


#553: sparse solvers segfault with singular matrices
---------------------------------------+------------------------------------
 Reporter:  insyzygy                   |        Owner:  somebody
     Type:  defect                     |       Status:  new     
 Priority:  normal                     |    Milestone:  0.7     
Component:  scipy.linsolve             |      Version:          
 Severity:  normal                     |   Resolution:          
 Keywords:  SuperLU segfault singular  |  
---------------------------------------+------------------------------------
Comment (by nils):

 The current behavior is as follows
 {{{
 gdb /usr/bin/python
 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 "i586-suse-linux"...Using host libthread_db
 library "/lib/tls/libthread_db.so.1".

 (gdb) run test_superlu.py
 Starting program: /usr/bin/python test_superlu.py
 [Thread debugging using libthread_db enabled]
 [New Thread 1077176000 (LWP 8758)]
 /usr/lib/python2.4/site-packages/scipy/splinalg/dsolve/linsolve.py:64:
 SparseEfficiencyWarning: spsolve requires CSC or CSR matrix format
   warn('spsolve requires CSC or CSR matrix format',
 SparseEfficiencyWarning)
 warning: singular matrix
 zeroing nan and inf entries...
 warning: (almost) singular matrix! (estimated cond. number: inf)
 warning: singular matrix
 zeroing nan and inf entries...
 warning: (almost) singular matrix! (estimated cond. number: inf)

 Program received signal SIGSEGV, Segmentation fault.
 [Switching to Thread 1077176000 (LWP 8758)]
 0x424888ce in dpivotL (jcol=0, u=1, usepr=0xbfffe3d0, perm_r=0x824f658,
 iperm_r=0x40343840, iperm_c=0x40343840,
     pivrow=0xbfffe3d4, Glu=0x429f9a80, stat=0x40343840) at dpivotL.c:120
 120             perm_r[*pivrow] = jcol;
 Current language:  auto; currently c
 (gdb) bt
 #0  0x424888ce in dpivotL (jcol=0, u=1, usepr=0xbfffe3d0,
 perm_r=0x824f658, iperm_r=0x40343840, iperm_c=0x40343840,
     pivrow=0xbfffe3d4, Glu=0x429f9a80, stat=0x40343840) at dpivotL.c:120
 #1  0x4247e5b2 in dgstrf (options=0xbfffe570, A=0xbfffe490, drop_tol=0,
 relax=5, panel_size=10, etree=0x820f748,
     work=0x40343840, lwork=1077164096, perm_c=0x8275d50, perm_r=0x824f658,
 L=0xbfffe5d0, U=0xbfffe5b0,
     stat=0xbfffe550, info=0xbfffe52c) at dgstrf.c:310
 #2  0x42461fe3 in dgssv (options=0xbfffe570, A=0xbfffe610,
 perm_c=0x8275d50, perm_r=0x824f658, L=0xbfffe5d0,
     U=0xbfffe5b0, B=0xbfffe5f0, stat=0xbfffe550, info=0xbfffe52c) at
 dgssv.c:204
 #3  0x4245eec9 in Py_dgssv (self=0x0, args=0x40343840, kwdict=0x40343840)
 at _dsuperlumodule.c:94
 #4  0x4007b5d9 in PyCFunction_Call () from /usr/lib/libpython2.4.so.1.0
 #5  0x400b2a91 in PyEval_EvalFrame () from /usr/lib/libpython2.4.so.1.0
 #6  0x400b4bc1 in PyEval_EvalCodeEx () from /usr/lib/libpython2.4.so.1.0
 #7  0x400b2b65 in PyEval_EvalFrame () from /usr/lib/libpython2.4.so.1.0
 #8  0x400b4bc1 in PyEval_EvalCodeEx () from /usr/lib/libpython2.4.so.1.0
 #9  0x400b4e95 in PyEval_EvalCode () from /usr/lib/libpython2.4.so.1.0
 #10 0x400cf618 in run_node () from /usr/lib/libpython2.4.so.1.0
 #11 0x400d0db3 in PyRun_SimpleFileExFlags () from
 /usr/lib/libpython2.4.so.1.0
 #12 0x400d137a in PyRun_AnyFileExFlags () from
 /usr/lib/libpython2.4.so.1.0
 #13 0x400d750a in Py_Main () from /usr/lib/libpython2.4.so.1.0
 #14 0x0804871a in main (argc=1077164096, argv=0x40343840) at
 ccpython.cc:10
 }}}

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


More information about the Scipy-tickets mailing list