Thu May 2 22:18:28 1996 Andreas Schwab <schwab@issan.informatik.uni-dortmund.de>
authorroland <roland>
Fri, 3 May 1996 16:50:14 +0000 (16:50 +0000)
committerroland <roland>
Fri, 3 May 1996 16:50:14 +0000 (16:50 +0000)
* sysdeps/libm-ieee754/e_j0f.c (__ieee754_j0f, __ieee754_y0f):
Replace 0x80000000 by 0x48000000.
* sysdeps/libm-ieee754/e_j1f.c (__ieee754_j1f): Likewise.

sysdeps/libm-ieee754/e_j0f.c
sysdeps/libm-ieee754/e_j1f.c

index b4e97a2..eed171c 100644 (file)
@@ -79,7 +79,7 @@ static float zero = 0.0;
         * j0(x) = 1/sqrt(pi) * (P(0,x)*cc - Q(0,x)*ss) / sqrt(x)
         * y0(x) = 1/sqrt(pi) * (P(0,x)*ss + Q(0,x)*cc) / sqrt(x)
         */
-               if(ix>0x80000000) z = (invsqrtpi*cc)/__sqrtf(x);
+               if(ix>0x48000000) z = (invsqrtpi*cc)/__sqrtf(x);
                else {
                    u = pzerof(x); v = qzerof(x);
                    z = invsqrtpi*(u*cc-v*ss)/__sqrtf(x);
@@ -161,7 +161,7 @@ v04  =  4.4111031494e-10; /* 0x2ff280c2 */
                     if ((s*c)<zero) cc = z/ss;
                     else            ss = z/cc;
                 }
-                if(ix>0x80000000) z = (invsqrtpi*ss)/__sqrtf(x);
+                if(ix>0x48000000) z = (invsqrtpi*ss)/__sqrtf(x);
                 else {
                     u = pzerof(x); v = qzerof(x);
                     z = invsqrtpi*(u*ss+v*cc)/__sqrtf(x);
index 06384bd..e6f14a1 100644 (file)
@@ -80,7 +80,7 @@ static float zero    = 0.0;
         * j1(x) = 1/sqrt(pi) * (P(1,x)*cc - Q(1,x)*ss) / sqrt(x)
         * y1(x) = 1/sqrt(pi) * (P(1,x)*ss + Q(1,x)*cc) / sqrt(x)
         */
-               if(ix>0x80000000) z = (invsqrtpi*cc)/__sqrtf(y);
+               if(ix>0x48000000) z = (invsqrtpi*cc)/__sqrtf(y);
                else {
                    u = ponef(y); v = qonef(y);
                    z = invsqrtpi*(u*cc-v*ss)/__sqrtf(y);