[SciPy-dev] FFT with DJBFFT failing

John Byrnes byrnes at bu.edu
Wed Mar 29 18:28:31 EST 2006


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)




More information about the SciPy-Dev mailing list