On Thu, Mar 26, 2009 at 9:06 PM, Charles R Harris <charlesr.har...@gmail.com > wrote:
> > > 2009/3/26 Christian Marquardt <christ...@marquardt.sc> > >> Oh sorry - you are right (too late in the night here in Europe). >> >> The output is similar in all four cases - it looks like >> >> AssertionError: >> Arrays are not almost equal >> >> (mismatch 100.0%) >> x: array([ 4.60555124+0.j, -2.60555124+0.j], dtype=complex64) >> y: array([-2.60555124 +1.11022302e-16j, 4.60555124 -1.11022302e-16j], >> dtype=complex64) >> >> Are x and y the expected and actual results? That would just show that >> there >> are small rounding errors in the imaginary part, and that MKL returns the >> results >> in another order, no? >> > > Looks like a sorting error, the eigen values should be sorted. So it looks > like a buggy test from here. Having an imaginary part to the eigenvalues > returned by a routine that is supposed to process Hermitean matrices doesn't > look right, but the errors are in the double precision range, which is > pretty good for float32. > > I think we need a ticket to fix those tests. > Can you post the actual error messages? It will make it easier to find where the failure is. Chuck
_______________________________________________ Numpy-discussion mailing list Numpy-discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion
