[Scipy-tickets] [SciPy] #1113: scipy.linalg.eig crash on LinAlgError

SciPy Trac scipy-tickets@scipy....
Fri Feb 12 17:45:45 CST 2010


#1113: scipy.linalg.eig crash on LinAlgError
--------------------+-------------------------------------------------------
 Reporter:  anjiro  |       Owner:  somebody
     Type:  defect  |      Status:  new     
 Priority:  normal  |   Milestone:  0.8.0   
Component:  Other   |     Version:  0.7.0   
 Keywords:          |  
--------------------+-------------------------------------------------------
 The following code:


 {{{
 import scipy.linalg
 from scipy import arange, array, dot
 C = array([[0,0],[0,0]])
 D = array([[1,2,3], [1,2,3]])
 S = dot(D.T,D)
 scipy.linalg.eig(S,C)
 }}}

 gives this output:

 {{{
 Traceback (most recent call last):
   File "t.py", line 7, in <module>
     scipy.linalg.eig(S,C)
   File "/Library/Python/2.6/site-
 packages/scipy-0.7.1-py2.6-macosx-10.6-universal.egg/scipy/linalg/decomp.py",
 line 157, in eig
     return _geneig(a1,b,left,right,overwrite_a,overwrite_b)
   File "/Library/Python/2.6/site-
 packages/scipy-0.7.1-py2.6-macosx-10.6-universal.egg/scipy/linalg/decomp.py",
 line 84, in _geneig
     if info>0: raise LinAlgError,"generalized eig algorithm did not
 converge"
 numpy.linalg.linalg.LinAlgError: generalized eig algorithm did not
 converge
 }}}

 followed by a segmentation fault. Under gdb, I get the following
 traceback:

 {{{
 Program received signal EXC_BAD_ACCESS, Could not access memory.
 Reason: 13 at address: 0x0000000000000000
 0x000000010004e0ae in PyType_Modified ()
 (gdb) bt
 #0  0x000000010004e0ae in PyType_Modified ()
 #1  0x000000010004e0c0 in PyType_Modified ()
 #2  0x000000010004e0c0 in PyType_Modified ()
 #3  0x000000010004e0c0 in PyType_Modified ()
 #4  0x000000010004e15d in PyType_ClearCache ()
 #5  0x00000001000a30c5 in Py_Finalize ()
 #6  0x00000001000b03b1 in Py_Main ()
 #7  0x0000000100000e6c in ?? ()
 }}}

 I'm running scipy 0.7.1, under Python 2.6.1 on OSX (it's the Apple
 python).

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


More information about the Scipy-tickets mailing list