aboutsummaryrefslogtreecommitdiff
path: root/sysdeps/ieee754/flt-32
diff options
context:
space:
mode:
authorAndreas Jaeger <aj@suse.de>2002-08-23 09:29:27 +0000
committerAndreas Jaeger <aj@suse.de>2002-08-23 09:29:27 +0000
commitc0ad824e402c7afde1b46be1620c70107aa3eda9 (patch)
tree193845a8400499e3d01009631761793183a04bfe /sysdeps/ieee754/flt-32
parenta7f5255ad25b5b7b637e3f2c1a648da2df0f1efc (diff)
downloadglibc-c0ad824e402c7afde1b46be1620c70107aa3eda9.tar
glibc-c0ad824e402c7afde1b46be1620c70107aa3eda9.tar.gz
glibc-c0ad824e402c7afde1b46be1620c70107aa3eda9.tar.bz2
glibc-c0ad824e402c7afde1b46be1620c70107aa3eda9.zip
Update.
2002-08-23 Andreas Jaeger <aj@suse.de> * locale/loadarchive.c: Include <stdlib.h> for prototypes of alloca, free, qsort. * dlfcn/bug-dlsym1-lib1.c (foo): Add prototype to avoid warning. * sysdeps/ieee754/flt-32/e_jnf.c: Remove unused variable. * sysdeps/ieee754/ldbl-96/e_j0l.c: Likewise. * sysdeps/generic/fgetxattr.c: New file. * sysdeps/generic/flistxattr.c: New file. * sysdeps/generic/fremovexattr.c: New file. * sysdeps/generic/fsetxattr.c: New file. * sysdeps/generic/getxattr.c: New file. * sysdeps/generic/listxattr.c: New file. * sysdeps/generic/lgetxattr.c: New file. * sysdeps/generic/llistxattr.c: New file. * sysdeps/generic/lremovexattr.c: New file. * sysdeps/generic/lsetxattr.c: New file. * sysdeps/generic/removexattr.c: New file. * sysdeps/generic/setxattr.c: New file. * misc/sys/xattr.h: New file, contributed by Christoph Hellwig <hch@lst.de>. * misc/Makefile (headers): Add sys/xattr.h. (routines): Add *xattr calls. * misc/Versions: Add *xattr calls with version 2.3. * sysdeps/unix/sysv/linux/m68k/Versions: Remove *xattr syscalls. * sysdeps/unix/sysv/linux/s390/s390-32/Versions: Likewise. * sysdeps/unix/sysv/linux/sparc/sparc32/Versions: Likewise. * sysdeps/unix/sysv/linux/sparc/sparc64/Versions: Likewise. * sysdeps/unix/sysv/linux/x86_64/Versions: Likewise. * sysdeps/unix/sysv/linux/ia64/Versions: Likewise. * sysdeps/unix/sysv/linux/alpha/Versions: Likewise. * sysdeps/unix/sysv/linux/arm/Versions: Likewise. * sysdeps/powerpc/elf/libc-start.c * sysdeps/unix/sysv/linux/powerpc/dl-sysdep.c
Diffstat (limited to 'sysdeps/ieee754/flt-32')
-rw-r--r--sysdeps/ieee754/flt-32/e_jnf.c19
1 files changed, 9 insertions, 10 deletions
diff --git a/sysdeps/ieee754/flt-32/e_jnf.c b/sysdeps/ieee754/flt-32/e_jnf.c
index 9e5279c30a..34c4d95a32 100644
--- a/sysdeps/ieee754/flt-32/e_jnf.c
+++ b/sysdeps/ieee754/flt-32/e_jnf.c
@@ -25,7 +25,6 @@ static const float
#else
static float
#endif
-invsqrtpi= 5.6418961287e-01, /* 0x3f106ebb */
two = 2.0000000000e+00, /* 0x40000000 */
one = 1.0000000000e+00; /* 0x3F800000 */
@@ -62,7 +61,7 @@ static float zero = 0.0000000000e+00;
if(n==1) return(__ieee754_j1f(x));
sgn = (n&1)&(hx>>31); /* even n -- 0, odd n -- sign(x) */
x = fabsf(x);
- if(ix==0||ix>=0x7f800000) /* if x is 0 or inf */
+ if(ix==0||ix>=0x7f800000) /* if x is 0 or inf */
b = zero;
else if((float)n<=x) {
/* Safe to use J(n+1,x)=2n/x *J(n,x)-J(n-1,x) */
@@ -90,11 +89,11 @@ static float zero = 0.0000000000e+00;
}
} else {
/* use backward recurrence */
- /* x x^2 x^2
+ /* x x^2 x^2
* J(n,x)/J(n-1,x) = ---- ------ ------ .....
* 2n - 2(n+1) - 2(n+2)
*
- * 1 1 1
+ * 1 1 1
* (for large x) = ---- ------ ------ .....
* 2n 2(n+1) 2(n+2)
* -- - ------ - ------ -
@@ -107,7 +106,7 @@ static float zero = 0.0000000000e+00;
* 1
* w - -----------------
* 1
- * w+h - ---------
+ * w+h - ---------
* w+2h - ...
*
* To determine how many terms needed, let
@@ -144,15 +143,15 @@ static float zero = 0.0000000000e+00;
v = two/x;
tmp = tmp*__ieee754_logf(fabsf(v*tmp));
if(tmp<(float)8.8721679688e+01) {
- for(i=n-1,di=(float)(i+i);i>0;i--){
+ for(i=n-1,di=(float)(i+i);i>0;i--){
temp = b;
b *= di;
b = b/x - a;
a = temp;
di -= two;
- }
+ }
} else {
- for(i=n-1,di=(float)(i+i);i>0;i--){
+ for(i=n-1,di=(float)(i+i);i>0;i--){
temp = b;
b *= di;
b = b/x - a;
@@ -164,9 +163,9 @@ static float zero = 0.0000000000e+00;
t /= b;
b = one;
}
- }
+ }
}
- b = (t*__ieee754_j0f(x)/b);
+ b = (t*__ieee754_j0f(x)/b);
}
}
if(sgn==1) return -b; else return b;