On Sun, Aug 8, 2010 at 7:20 PM, Pauli Virtanen <p...@iki.fi> wrote: > Sun, 08 Aug 2010 15:09:00 +0800, Ralf Gommers wrote: > > On Mon, Aug 2, 2010 at 10:30 AM, Alan G Isaac <alan.is...@gmail.com> > [clip] > >> ====================================================================== > >> FAIL: test_umath_complex.TestCsqrt.test_special_values(<ufunc 'sqrt'>, > >> 1, inf, inf, inf) > >> test_umath_complex.TestCsqrt.test_special_values(<ufunc 'sqrt'>, 1, > >> inf, inf, inf) > >> ---------------------------------------------------------------------- > >> Traceback (most recent call last): > >> File "C:\Python27\lib\site-packages\nose\case.py", line 183, in > >> runTest > >> self.test(*self.arg) > >> File > >> "C:\Python27\lib\site-packages\numpy\core\tests\test_umath_complex.py", > >> line 539, i > >> n check_complex_value > >> assert_equal(f(z1), z2) > >> File "C:\Python27\lib\site-packages\numpy\testing\utils.py", line > >> 256, > >> in assert_equal > >> return assert_array_equal(actual, desired, err_msg, verbose) > >> File "C:\Python27\lib\site-packages\numpy\testing\utils.py", line > >> 686, > >> in assert_array_e > >> qual > >> verbose=verbose, header='Arrays are not equal') > >> File "C:\Python27\lib\site-packages\numpy\testing\utils.py", line > >> 596, > >> in assert_array_c > >> ompare > >> raise AssertionError(msg) > >> AssertionError: > >> Arrays are not equal > >> > >> (x and y nan location mismatch [ True], False mismatch) > >> x: array([ inf nanj]) > >> y: array((inf+infj)) > >> >> raise AssertionError('\nArrays are not equal\n\n(x and y nan > >> >> location > >> mismatch [ True] > >> , False mismatch)\n x: array([ inf nanj])\n y: array((inf+infj))') > [clip] > > Pauli or Charles, can you please have a look at these? Looks like > > that's related to your recent work on nans/infs. > > These errors probably come from the fact that the platform's C library > does not handle special nan/inf values in csqrt in the way C99 standard > requires. The recent changes in umath since 1.4.x are unrelated. > > I think the support for platform-provided complex functions was enabled > in r7986. > > To be sure, we'd need to see the build log. > > Ticket opened with 2.7 and 2.6 build logs attached, http://projects.scipy.org/numpy/ticket/1574. Looks like it tries to use C99: #define NPY_USE_C99_COMPLEX #define NPY_HAVE_COMPLEX_DOUBLE 1 #define NPY_HAVE_COMPLEX_FLOAT 1 #define NPY_HAVE_COMPLEX_LONG_DOUBLE 1 #define NPY_USE_C99_FORMATS 1
Ralf
_______________________________________________ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion