[Scipy-tickets] [SciPy] #1723: scipy.test() crashes part way through

SciPy Trac scipy-tickets@scipy....
Sun Aug 26 19:05:18 CDT 2012


#1723: scipy.test() crashes part way through
---------------------+------------------------------------------------------
 Reporter:  chrisfs  |       Owner:  somebody   
     Type:  defect   |      Status:  new        
 Priority:  normal   |   Milestone:  Unscheduled
Component:  Other    |     Version:  0.10.0     
 Keywords:           |  
---------------------+------------------------------------------------------

Comment(by chrisfs):

 Replying to [comment:3 rgommers]:
 > This is a stand-alone version of the test that crashes
 > {{{
 >
 > import numpy as np
 > from numpy.testing import assert_array_almost_equal
 > from scipy.linalg import qr
 >
 > np.random.seed(1234)
 >
 > m = 200
 > n = 100
 > for k in range(2):
 >     a = np.random.random([m, n])
 >     q,r = qr(a)
 >     assert_array_almost_equal(np.dot(np.transpose(q),q), np.identity(m))
 >     assert_array_almost_equal(np.dot(q,r), a)
 > }}}
 > Does that crash too? If so, can you confirm that it crashes on the call
 to {{{qr}}}? You can do so by executing this line by line in IPython for
 example.

 When I go through line by line (and take out the FOR statement, so I can
 run through that manually), it crashes at
 assert_array_almost_equal(np.dot(np.transpose(q),q), np.identity(m))
 Keeping the For Statement and stepping through the code using the debugger
 in IDLE gives a crash at the same place.
 The debugger says:
 '__main__'.<module>(),line 12:
 assert_array_almost_equal(np.dot(np.transpose(q),q), np.identity(m))
 > 'numpy,core,fromnumeric'transpose(), line 487: return transpose(axes)

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


More information about the Scipy-tickets mailing list