Thu, 28 May 2009 16:29:38 +0900, David Cournapeau wrote: [clip]
Ok, there are 2 failures on mac os x in scipy.special. test_yn_zeros and test_ynp_zeros are (way) below the tolerance you set up, I put the output of the test in #953.
I think I found a single-line fix: --- a.f 2009-05-29 00:12:14.000000000 +0300 +++ b.f 2009-05-29 21:27:18.246736714 +0300 @@ -159,7 +159,7 @@ BS=BS+2.0D0*F SU=SU+(-1)**(K/2)*F/K ELSE IF (K.GT.1) THEN - SV=SV+(-1)**(K/2)*K/(K*K-1.0)*F + SV=SV+(-1)**(K/2)*K*F/(K*K-1.0) ENDIF F2=F1 15 F1=F Fixes the issue at least with Gfortran ffloat-store. Who said multiplication is commutative :) -- Pauli Virtanen