[SciPy-dev] FFT with DJBFFT failing
John Byrnes
byrnes at bu.edu
Thu Mar 30 00:47:18 CST 2006
I suppose I should give some system details.
Pentium 4 - 32 bit, Linux, Python 2.4, DJBFFT version 0.76 compiled from source.
GCC 4.0 for C code, GCC 3.4.5 for fortran.
Thanks!
John
On Wed, Mar 29, 2006 at 06:28:31PM -0500, John Byrnes wrote:
> Hello all,
>
> I've found that DJBFFT causes scipy to fail several tests. These
> failures only occur when DJBFFT is found by distutils when building.
>
> I'm running the current SVN copies of both scipy and numpy.
>
> Regards,
> John
>
>
> Report is below:
>
> ======================================================================
> FAIL: check_normal (scipy.stats.tests.test_morestats.test_anderson)
> ----------------------------------------------------------------------
> Traceback (most recent call last):
> File
> "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/stats/tests/test_morestats.py",
> line 51, in check_normal
> assert_array_less(A, crit[-2:])
> File
> "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py",
> line 255, in assert_array_less
> assert cond,\
> AssertionError:
> Arrays are not less-ordered (mismatch 50.0%):
> Array 1: 0.90160265659743288
> Array 2: [ 0.858 1.0209999999999999]
>
>
> ======================================================================
> FAIL: check_definition (scipy.fftpack.tests.test_pseudo_diffs.test_diff)
> ----------------------------------------------------------------------
> Traceback (most recent call last):
> File
> "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_pseudo_diffs.py",
> line 85, in check_definition
> assert_array_almost_equal(diff(sin(x)),direct_diff(sin(x)))
> File
> "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py",
> line 231, in assert_array_almost_equal
> assert cond,\
> AssertionError:
> Arrays are not almost equal (mismatch 87.5%):
> Array 1: [ 1.0000000e+00 9.2387953e-01 7.0710678e-01
> 3.8268343e-01
> 8.8009815e-17 -3.8268343e-01 -7.0710678e-01 -9.23...
> Array 2: [ -0.0000000e+00 5.0000000e-01 7.7662794e-17
> 1.2967736e-16
> 1.2490009e-16 2.1612893e-16 1.4972165e-16 -9.71...
>
>
> ======================================================================
> FAIL: check_definition
> (scipy.fftpack.tests.test_pseudo_diffs.test_hilbert)
> ----------------------------------------------------------------------
> Traceback (most recent call last):
> File
> "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_pseudo_diffs.py",
> line 301, in check_definition
> assert_array_almost_equal (y,y1)
> File
> "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py",
> line 231, in assert_array_almost_equal
> assert cond,\
> AssertionError:
> Arrays are not almost equal (mismatch 87.5%):
> Array 1: [ 1.0000000e+00 9.2387953e-01 7.0710678e-01
> 3.8268343e-01
> 8.2304423e-17 -3.8268343e-01 -7.0710678e-01 -9.23...
> Array 2: [ -0.0000000e+00 +0.0000000e+00j 5.0000000e-01
> -7.3105720e-17j
> 3.8831397e-17 -2.4238713e-17j 4.3225785e-17 -6.813...
>
>
> ======================================================================
> FAIL: check_random_even
> (scipy.fftpack.tests.test_pseudo_diffs.test_hilbert)
> ----------------------------------------------------------------------
> Traceback (most recent call last):
> File
> "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_pseudo_diffs.py",
> line 332, in check_random_even
> assert_array_almost_equal(direct_hilbert(direct_ihilbert(f)),f)
> File
> "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py",
> line 231, in assert_array_almost_equal
> assert cond,\
> AssertionError:
> Arrays are not almost equal (mismatch 100.0%):
> Array 1: [ 0. -0.j -0.0003392+0.0003636j
> -0.0010602-0.0002132j
> -0.0005238+0.0002887j 0.0003091-0.0020875j 0.00076...
> Array 2: [-0.1106071 -0.497691 0.4783423 -0.337373
> 0.2458101 0.0210909
> -0.4631084 0.192853 0.4792201 -0.3906347 0.22212...
>
>
> ======================================================================
> FAIL: check_tilbert_relation
> (scipy.fftpack.tests.test_pseudo_diffs.test_hilbert)
> ----------------------------------------------------------------------
> Traceback (most recent call last):
> File
> "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_pseudo_diffs.py",
> line 311, in check_tilbert_relation
> assert_array_almost_equal (y,y1)
> File
> "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py",
> line 231, in assert_array_almost_equal
> assert cond,\
> AssertionError:
> Arrays are not almost equal (mismatch 62.5%):
> Array 1: [ 1.0000000e+00 6.5328148e-01 2.1280927e-16
> -2.7059805e-01
> -1.0626020e-16 2.7059805e-01 1.5393080e-16 -6.53...
> Array 2: [ -0.0000000e+00 +0.0000000e+00j 2.5000000e-01
> -3.5284004e-17j
> -1.1905260e-18 +4.4746813e-17j 2.5000000e-01 -9.941...
>
>
> ======================================================================
> FAIL: check_definition
> (scipy.fftpack.tests.test_pseudo_diffs.test_ihilbert)
> ----------------------------------------------------------------------
> Traceback (most recent call last):
> File
> "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_pseudo_diffs.py",
> line 370, in check_definition
> assert_array_almost_equal (y,y1)
> File
> "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py",
> line 231, in assert_array_almost_equal
> assert cond,\
> AssertionError:
> Arrays are not almost equal (mismatch 87.5%):
> Array 1: [ -1.0000000e+00 -9.2387953e-01 -7.0710678e-01
> -3.8268343e-01
> -8.2304423e-17 3.8268343e-01 7.0710678e-01 9.23...
> Array 2: [ 0.0000000e+00 -0.0000000e+00j -5.0000000e-01
> +7.3105720e-17j
> -3.8831397e-17 +2.4238713e-17j -4.3225785e-17 +6.813...
>
>
> ======================================================================
> FAIL: check_itilbert_relation
> (scipy.fftpack.tests.test_pseudo_diffs.test_ihilbert)
> ----------------------------------------------------------------------
> Traceback (most recent call last):
> File
> "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_pseudo_diffs.py",
> line 380, in check_itilbert_relation
> assert_array_almost_equal (y,y1)
> File
> "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py",
> line 231, in assert_array_almost_equal
> assert cond,\
> AssertionError:
> Arrays are not almost equal (mismatch 62.5%):
> Array 1: [ -1.0000000e+00 -6.5328148e-01 -2.1280927e-16
> 2.7059805e-01
> 1.0626020e-16 -2.7059805e-01 -1.5393080e-16 6.53...
> Array 2: [ 0.0000000e+00 -0.0000000e+00j -2.5000000e-01
> +3.5284004e-17j
> 1.1905260e-18 -4.4746813e-17j -2.5000000e-01 +9.941...
>
>
> ======================================================================
> FAIL: check_definition
> (scipy.fftpack.tests.test_pseudo_diffs.test_itilbert)
> ----------------------------------------------------------------------
> Traceback (most recent call last):
> File
> "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_pseudo_diffs.py",
> line 288, in check_definition
> assert_array_almost_equal (y,y1)
> File
> "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py",
> line 231, in assert_array_almost_equal
> assert cond,\
> AssertionError:
> Arrays are not almost equal (mismatch 87.5%):
> Array 1: [ -9.9667995e-02 -9.2081220e-02 -7.0475915e-02
> -3.8141290e-02
> -7.1261911e-18 3.8141290e-02 7.0475915e-02 9.20...
> Array 2: [ -0.0000000e+00 +0.0000000e+00j -4.9833997e-02
> +7.2863005e-18j
> -7.6643594e-18 +4.7841238e-18j -1.2592216e-17 +1.984...
>
>
> ======================================================================
> FAIL: check_definition
> (scipy.fftpack.tests.test_pseudo_diffs.test_shift)
> ----------------------------------------------------------------------
> Traceback (most recent call last):
> File
> "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_pseudo_diffs.py",
> line 390, in check_definition
> assert_array_almost_equal(shift(sin(x),a),direct_shift(sin(x),a))
> File
> "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py",
> line 231, in assert_array_almost_equal
> assert cond,\
> AssertionError:
> Arrays are not almost equal (mismatch 100.0%):
> Array 1: [ 0.0998334 0.1968802 0.2920308 0.3843691
> 0.4730057 0.5570869
> 0.6358031 0.7083962 0.7741671 0.8324823 0.88278...
> Array 2: [ 1.2321788e-17 4.9916708e-02 2.4648710e-19
> 8.8465810e-18
> 5.9412517e-18 2.0994071e-17 1.4776750e-17 1.22...
>
>
> ======================================================================
> FAIL: check_definition
> (scipy.fftpack.tests.test_pseudo_diffs.test_tilbert)
> ----------------------------------------------------------------------
> Traceback (most recent call last):
> File
> "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_pseudo_diffs.py",
> line 227, in check_definition
> assert_array_almost_equal (y,y1)
> File
> "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py",
> line 231, in assert_array_almost_equal
> assert cond,\
> AssertionError:
> Arrays are not almost equal (mismatch 87.5%):
> Array 1: [ 1.0033311e+01 9.2695708e+00 7.0946223e+00
> 3.8395819e+00
> 1.1013865e-15 -3.8395819e+00 -7.0946223e+00 -9.26...
> Array 2: [ -0.0000000e+00 +0.0000000e+00j 5.0166556e+00
> -7.3349243e-16j
> 1.9673887e-16 -1.2280519e-16j 1.4838281e-16 -2.338...
>
>
> ======================================================================
> FAIL: check_random_even
> (scipy.fftpack.tests.test_pseudo_diffs.test_tilbert)
> ----------------------------------------------------------------------
> Traceback (most recent call last):
> File
> "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_pseudo_diffs.py",
> line 240, in check_random_even
> assert_array_almost_equal(direct_tilbert(direct_itilbert(f,h),h),f)
> File
> "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py",
> line 231, in assert_array_almost_equal
> assert cond,\
> AssertionError:
> Arrays are not almost equal (mismatch 100.0%):
> Array 1: [ -0.0000000e+00 +0.0000000e+00j 1.4398428e-03
> +1.8373263e-03j
> -8.7680356e-04 +9.9658995e-04j -7.3017115e-04 +1.287...
> Array 2: [ 0.5075109 -0.2380892 -0.2038641 -0.3047494
> -0.2114436 0.4541402
> -0.1193784 -0.0768124 -0.1561373 0.1053 0.03414...
>
>
> ======================================================================
> FAIL: check_definition (scipy.fftpack.tests.test_basic.test_fft)
> ----------------------------------------------------------------------
> Traceback (most recent call last):
> File
> "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_basic.py",
> line 98, in check_definition
> assert_array_almost_equal(y,y1)
> File
> "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py",
> line 231, in assert_array_almost_equal
> assert cond,\
> AssertionError:
> Arrays are not almost equal (mismatch 100.0%):
> Array 1: [ 1.+0.j 2.+0.j 3.+0.j 4.+1.j 1.+0.j 2.+0.j
> 3.+0.j 4.+2.j]
> Array 2: [ 20. +3.j -0.7071068+0.7071068j -7.
> +4.j
> -0.7071068-0.7071068j -4. -3.j 0.707106...
>
>
> ======================================================================
> FAIL: check_definition (scipy.fftpack.tests.test_basic.test_fftn)
> ----------------------------------------------------------------------
> Traceback (most recent call last):
> File
> "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_basic.py",
> line 431, in check_definition
> assert_array_almost_equal(fftn(x),direct_dftn(x))
> File
> "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py",
> line 231, in assert_array_almost_equal
> assert cond,\
> AssertionError:
> Arrays are not almost equal (mismatch 100.0%):
> Array 1: [[[[ 5.9716569e+02 +0.j -4.7701973e+00
> -6.248201j
> 7.2274468e+00 +2.123387j ..., -1.8242164e+01 +4.20633...
> Array 2: [[[[ 1.4378537e+02+0.j
> 3.8003180e-01-0.9697448j
> 6.8443765e+00+3.6280536j ..., -2.9207977e+00+2.7763892j...
>
>
> ======================================================================
> FAIL: check_definition (scipy.fftpack.tests.test_basic.test_ifft)
> ----------------------------------------------------------------------
> Traceback (most recent call last):
> File
> "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_basic.py",
> line 183, in check_definition
> assert_array_almost_equal(y,y1)
> File
> "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py",
> line 231, in assert_array_almost_equal
> assert cond,\
> AssertionError:
> Arrays are not almost equal (mismatch 100.0%):
> Array 1: [ 0.125+0.j 0.25 +0.j 0.375+0.j 0.5
> +0.125j 0.125+0.j
> 0.25 +0.j 0.375+0.j 0.5 +0.25j ]
> Array 2: [ 2.5 +0.375j 0.0883883+0.0883883j -0.125
> -0.5j
> 0.0883883-0.0883883j -0.5 -0.375j -0.0883883-0.0...
>
>
> ======================================================================
> FAIL: check_random_complex (scipy.fftpack.tests.test_basic.test_ifft)
> ----------------------------------------------------------------------
> Traceback (most recent call last):
> File
> "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_basic.py",
> line 211, in check_random_complex
> assert_array_almost_equal (ifft(fft(x)),x)
> File
> "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py",
> line 231, in assert_array_almost_equal
> assert cond,\
> AssertionError:
> Arrays are not almost equal (mismatch 100.0%):
> Array 1: [ 0.000438 +1.2422571e-03j 0.0081279 +1.3153027e-03j
> 0.0058322 +2.5472850e-03j 0.0098841 +1.5452180e-02j
> 0.007895...
> Array 2: [ 0.0280313+0.0795045j 0.5201845+0.0841794j
> 0.3732589+0.1630262j
> 0.6325855+0.9889395j 0.5052878+0.7382801j 0.12457...
>
>
> ======================================================================
> FAIL: check_random_real (scipy.fftpack.tests.test_basic.test_ifft)
> ----------------------------------------------------------------------
> Traceback (most recent call last):
> File
> "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_basic.py",
> line 217, in check_random_real
> assert_array_almost_equal (ifft(fft(x)),x)
> File
> "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py",
> line 231, in assert_array_almost_equal
> assert cond,\
> AssertionError:
> Arrays are not almost equal (mismatch 100.0%):
> Array 1: [ 0.4665868+0.j 0.0158205-0.0104689j
> -0.0024869-0.0038622j
> 0.0255389-0.0128942j 0.0246149-0.0121173j 0.00113...
> Array 2: [ 0.5673542 0.8727196 0.7445887 0.5776637 0.910135
> 0.2396252
> 0.4235496 0.0873539 0.9231279 0.0897232 0.15871...
>
>
> ======================================================================
> FAIL: check_definition (scipy.fftpack.tests.test_basic.test_ifftn)
> ----------------------------------------------------------------------
> Traceback (most recent call last):
> File
> "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_basic.py",
> line 601, in check_definition
> assert_array_almost_equal(ifftn(x),direct_idftn(x))
> File
> "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py",
> line 231, in assert_array_almost_equal
> assert cond,\
> AssertionError:
> Arrays are not almost equal (mismatch 100.0%):
> Array 1: [[[[ 4.9883390e-01 +0.0000000e+00j 7.9039717e-03
> +1.7999071e-03j
> -2.5307390e-03 -2.3994048e-03j ..., 5.8921717...
> Array 2: [[[[ 1.2726888e-01 +0.0000000e+00j 8.6321110e-04
> +2.0907331e-03j
> 2.2428627e-03 +4.2203592e-04j ..., -2.7428294...
>
>
> ======================================================================
> FAIL: check_definition (scipy.fftpack.tests.test_basic.test_irfft)
> ----------------------------------------------------------------------
> Traceback (most recent call last):
> File
> "/usr/lib/python2.4/site-packages/scipy-0.4.9.1788-py2.4-linux-i686.egg/scipy/fftpack/tests/test_basic.py",
> line 341, in check_definition
> assert_array_almost_equal(y,ifft(x1))
> File
> "/usr/lib/python2.4/site-packages/numpy-0.9.7.2301-py2.4-linux-i686.egg/numpy/testing/utils.py",
> line 231, in assert_array_almost_equal
> assert cond,\
> AssertionError:
> Arrays are not almost equal (mismatch 100.0%):
> Array 1: [ 2.625 -1.6856602 -0.375 -1.1856602 0.625
> 0.4356602 -0.375
> 0.9356602]
> Array 2: [ 0.125+0.j 0.25 +0.375j 0.5 +0.125j 0.25
> +0.375j 0.5 +0.j
> 0.25 -0.375j 0.5 -0.125j 0.25 -0.375j]
>
>
> ----------------------------------------------------------------------
> Ran 1508 tests in 2.893s
>
> FAILED (failures=18)
>
> _______________________________________________
> Scipy-dev mailing list
> Scipy-dev at scipy.net
> http://www.scipy.net/mailman/listinfo/scipy-dev
--
Anyone who is capable of getting themselves made President should on no
account be allowed to do the job.
-- Douglas Adams, "The Hitchhiker's Guide to the Galaxy"
More information about the Scipy-dev
mailing list