[SciPy-dev] fft segfault, 64 Bit Opteron

[Andrew Straw]

FYI, on my Athlon64 (Debian sarge, gcc 3.3, stock atlas), I'm not
getting segfaults using scipy.test(10,10) but only a single failure. 
This is a svn checkout last night, so it may lag slightly, but it seems
to indicate Arnd's segfaults may reflect a gcc vs. icc issue.

In [2]: scipy.__core_version__
Out[2]: '0.8.1.1611'

In [3]: scipy.__scipy_version__
Out[3]: '0.4.3.1479'

======================================================================
FAIL: check_odeint1 (scipy.integrate.test_integrate.test_odeint)
----------------------------------------------------------------------
Traceback (most recent call last):
  File
"/home/astraw/py24-amd64/lib/python2.4/site-packages/scipy/integrate/tests/test_integrate.py",
line 51, in check_odeint1
    assert res < 1.0e-6
AssertionError

----------------------------------------------------------------------
Ran 1238 tests in 60.816s

FAILED (failures=1)


Travis Oliphant wrote:

>Arnd Baecker wrote:
>
>  
>
>>Hi Travis,
>>
>>one more addition - the build log gives:
>>
>>In file included from build/src/scipy/base/src/umathmodule.c:8036:
>>scipy/base/src/ufuncobject.c: In function `PyUFunc_GenericFunction':
>>scipy/base/src/ufuncobject.c:1569: warning: passing arg 2 of pointer to
>>function from incompatible pointer type
>>gcc -pthread -shared
>>build/temp.linux-x86_64-2.4/build/src/scipy/base/src/umathm
>>
>> 
>>
>>    
>>
>Thanks much for this testing.
>
>Could you send the entire build log again.  Perhaps there is something 
>else.
>I've made a fix which should work better on 64-bit. 
>
>More 64-bit testing needed.
>
>One trick to test the buffered section of code using smaller arrays is 
>to set the buffer size to something very small (but a multiple of 16 --- 
>16 is the smallest).   For arrays smaller than the buffer size, array's 
>are just copied when a cast is needed.  But, for larger arrays, the 
>buffered code is exercised.
>
>For example:
>
>scipy.setbufsize(16)
>scipy.test(1,1)
>
>At some-point we might play with different buffer sizes to see if some 
>numbers are better than others.
>
>-Travis
>
>
>_______________________________________________
>Scipy-dev mailing list
>Scipy-dev at scipy.net
>http://www.scipy.net/mailman/listinfo/scipy-dev
>  
>