On Fri, Jan 14, 2011 at 12:48 PM, Pauli Virtanen <pav@iki.fi> wrote:
Fri, 14 Jan 2011 12:09:54 -0500, josef.pktd wrote:
compiled against numpy 1.5.1
Josef
====================================================================== FAIL: Test singular pair ---------------------------------------------------------------------- [clip] "...\scipy\linalg\tests\test_decomp.py", line 189, in _check_gen_eig err_msg=msg)
This is actually a bit strange. The test does this:
w, vr = eig(A,B) wt = eigvals(A,B) # ... some stuff that does not touch w and wt # The failing assertion: assert_array_almost_equal(sort(w[isfinite(w)]), sort(wt[isfinite(wt)]), err_msg=msg)
The code paths in `eig` and `eigvals` are identical as far as the computation of `w` and `wt` is concerned.
So perhaps this is another manifestation of the no-determinism of linear algebra which I recall occurred on your computer? It seems this particular is supposed to be ill-conditioned and susceptible to rounding error...
***
It would be nice to find out what exactly is causing the non-determinism, since it seems you are not the only one seeing it (other machines where this occurs were reported on the Numpy ML).
http://permalink.gmane.org/gmane.comp.python.numeric.general/42082
Are you able to reproduce it using only the `dot` function? What happens if you then use the non-BLAS `numpy.core.multiarray.dot` instead?
I can try a bit later, I'm trying to figure out a stats bug and don't want to mess with my scipy since it takes me more than an hour to build. Looking at the test and try the non atlas dot should be fine. Josef
-- Pauli Virtanen
_______________________________________________ SciPy-Dev mailing list SciPy-Dev@scipy.org http://mail.scipy.org/mailman/listinfo/scipy-dev