[SciPy-dev] 0.7.1: new schedule

[Pauli Virtanen]

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