
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? -- Pauli Virtanen