[SciPy-user] SciPy with gcc4 and gfortran on OS X

[Zachary Pincus]

Robert -

Thanks for the feedback about the errors, etc.
Do you have any feel for the persistent errors on scipy.test() that  
I've been seeing with gfortran? You commented on the sporadic errors  
earlier, but the consistent ones worry me the most, given that it  
appears to be a simple dot product (on complex numbers, reproduced  
below for convenience) that is failing badly, and consistently. Is  
this indicative of gfortran problems, do you think, or something else  
that I can hunt down?

Zach

PS. Thanks for your time here, especially since this is all sort of  
tangential (because gfortran isn't a necessity on PPC Apple hardware,  
and apparently these problems aren't happening on Intel Apple scipy  
builds with gfortran (??) ).

That dot product problem:
FAIL: check_dot (scipy.lib.tests.test_blas.test_fblas1_simple)
----------------------------------------------------------------------
Traceback (most recent call last):
    File "/Library/Frameworks/Python.framework/Versions/2.4/lib/
python2.4/site-packages/scipy/lib/blas/tests/test_blas.py", line 76,
in check_dot
      assert_almost_equal(f([3j,-4,3-4j],[2,3,1]),-9+2j)
    File "/Library/Frameworks/Python.framework/Versions/2.4/lib/
python2.4/site-packages/numpy/testing/utils.py", line 148, in
assert_almost_equal
      assert round(abs(desired - actual),decimal) == 0, msg
AssertionError:
Items are not equal:
DESIRED: (-9+2j)
ACTUAL: (-1.9984917640686035-1.9984936714172363j)

>
>> SPORADIC FAILURES
>> ===================================================================== 
>> =
>> FAIL: check_expon (scipy.stats.tests.test_morestats.test_anderson)
>> --------------------------------------------------------------------- 
>> -
>> Traceback (most recent call last):
>>    File "/Library/Frameworks/Python.framework/Versions/2.4/lib/
>> python2.4/site-packages/scipy/stats/tests/test_morestats.py", line
>> 57, in check_expon
>>      assert_array_less(A, crit[-2:])
>>    File "/Library/Frameworks/Python.framework/Versions/2.4/lib/
>> python2.4/site-packages/numpy/testing/utils.py", line 255, in
>> assert_array_less
>>      assert cond,\
>> AssertionError:
>> Arrays are not less-ordered (mismatch 100.0%):
>>          Array 1: 2.1501866413808912
>>          Array 2: [ 1.587               1.9339999999999999]
>
> Well, this is a stochastic test; it is *supposed* to fail  
> sporadically. However,
> the mismatch in the array shapes is probably indicative of a real  
> bug. The
> anderson function is on the list for review, of course.
>
>   http://projects.scipy.org/scipy/scipy/ticket/159