[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