[Scipy-tickets] [SciPy] #1472: Recent Arpack changes segfaulting in test suite

SciPy Trac scipy-tickets@scipy....
Sun Jul 3 12:00:58 CDT 2011


#1472: Recent Arpack changes segfaulting in test suite
---------------------------------+------------------------------------------
 Reporter:  rgommers             |       Owner:  pv    
     Type:  defect               |      Status:  new   
 Priority:  high                 |   Milestone:  0.10.0
Component:  scipy.sparse.linalg  |     Version:  0.9.0 
 Keywords:  arpack               |  
---------------------------------+------------------------------------------

Comment(by rgommers):

 Snippet to trigger segfault:
 {{{
 import numpy as np
 from scipy.sparse.linalg.isolve import gmres

 np.random.seed(1234)
 data = np.random.rand(4, 4) + 1j*np.random.rand(4, 4)
 b  = np.arange(data.shape[0], dtype=float)
 x, info = gmres(data, b, x0=0*b, tol=1e-8)
 }}}

 Routine is named {{{zgmresrevcom_}}} in _iterative.so (use "otool -tV" on
 OS X instead of objdump.

 Setting the breakpoint and stepping from there ends with "warning:
 Couldn't find previous frame to set return breakpoint, stopping now."

 {{{
 $ gdb python
 ...
 (gdb) run -c "from scipy.sparse.linalg.isolve import _iterative;
 getattr(_iterative, 'zgmresrevcom')"
 Starting program:
 /Library/Frameworks/Python.framework/Versions/2.6/bin/python -c "from
 scipy.sparse.linalg.isolve import _iterative; getattr(_iterative,
 'zgmresrevcom')"
 Reading symbols for shared libraries +. done

 Program received signal SIGTRAP, Trace/breakpoint trap.
 0x8fe01030 in __dyld__dyld_start ()
 (gdb) c
 Continuing.
 Reading symbols for shared libraries warning: Could not find object file
 "/Volumes/Rivendell/Users/ronald/Projects/python/r264/Mac/BuildScript/_py/_bld/python/libpython2.6.a(getbuildinfo.o)"
 - no debug information available for
 "/Volumes/Rivendell/Users/ronald/Projects/python/r264/Modules/getbuildinfo.c".

 ...<lots more of the same noise>...

 (gdb) break zgmresrevcom_
 Breakpoint 1 at 0x189bcfb: file
 build/scons/scipy/sparse/linalg/isolve/GMRESREVCOM.f, line 325.
 (gdb) run -c "import numpy as np; from scipy.sparse.linalg.isolve import
 gmres; np.random.seed(1234); data = np.random.rand(4, 4) +
 1j*np.random.rand(4, 4); b = np.arange(data.shape[0], dtype=float); x,
 info = gmres(data, b, x0=0*b, tol=1e-8)"
 Starting program:
 /Library/Frameworks/Python.framework/Versions/2.6/bin/python -c "import
 numpy as np; from scipy.sparse.linalg.isolve import gmres;
 np.random.seed(1234); data = np.random.rand(4, 4) + 1j*np.random.rand(4,
 4); b = np.arange(data.shape[0], dtype=float); x, info = gmres(data, b,
 x0=0*b, tol=1e-8)"

 Program received signal SIGTRAP, Trace/breakpoint trap.
 0x8fe01030 in __dyld__dyld_start ()
 (gdb) s
 Single stepping until exit from function __dyld__dyld_start,
 which has no line number information.
 0x8fe015a6 in __dyld__ZN13dyldbootstrap5startEPK12macho_headeriPPKcl ()
 (gdb) s
 Single stepping until exit from function
 __dyld__ZN13dyldbootstrap5startEPK12macho_headeriPPKcl,
 which has no line number information.
 0x8fe01057 in __dyld__dyld_start ()
 (gdb) s
 Single stepping until exit from function __dyld__dyld_start,
 which has no line number information.
 warning: Couldn't find previous frame to set return breakpoint, stopping
 now.
 0x00001e80 in ?? ()
 (gdb) s
 Cannot find bounds of current function
 }}}

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


More information about the Scipy-tickets mailing list