[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