mirror of
https://sourceware.org/git/glibc.git
synced 2024-12-14 23:30:06 +00:00
b4012b7542
1997-04-21 13:25 Ulrich Drepper <drepper@cygnus.com> * manual/arith.texi: Add description for INFINITY, _Imaginary_I, fpclassify & friends, and complex number operations. Update various other math functions for ISO C 9X. * manual/math.texi: Update various entries for ISO C 9X. Add description for complex number functions. Add description of rand48 function family. * manual/string.h: Add description of a64l and l64a. * math/cmathcalls.h: Fix typo. * stdlib/a64l.c: Pretty printing. * stdlib/seed48_r.c: Also reset `a' and `c' to default values. * stdlib/srand48_r.c: Likewise. * stdlib/stdlib.h: Pretty printing. * sysdeps/i386/fpu/__math.h: Fix typo. * sysdeps/libm-ieee754/s_nearbyintf.c: Correctly name function. * sysdeps/libm-ieee754/s_nearbyintl.c: Likewise. 1997-04-19 22:16 Andreas Schwab <schwab@issan.informatik.uni-dortmund.de> * sysdeps/m68k/fpu/e_pow.c: Rewrite handling of integral exponent. 1997-04-18 19:34 Andreas Schwab <schwab@issan.informatik.uni-dortmund.de> * sysdeps/m68k/fpu/__math.h: Define optimized versions of isgreater, isgreaterequal, isless, islessequal, islessgreater, and isunordered. 1997-04-20 01:28 Richard Henderson <rth@tamu.edu> * rellns-sh: Handle files in the same directory correctly. 1997-04-20 11:22 Ulrich Drepper <drepper@cygnus.com> * csu/initfini.c: Place ALIGN instruction at correct positions. Patch by Richard Henderson <richard@twiddle.rth.home>. 1997-04-19 17:12 Ulrich Drepper <drepper@cygnus.com> * Make-dist: Don't automatically ignore .c files if the .S or .s file is ignored. * csu/Makefile (distribute): Add defs.awk. 1997-04-19 15:39 Ulrich Drepper <drepper@cygnus.com> * sysdeps/stub/shmat.c: Update to XPG4.2 interface. * sysdeps/stub/shmdt.c: Likewise. Reported by Thomas Bushnell, n/BSG. 1997-04-19 13:22 Ulrich Drepper <drepper@cygnus.com> * manual/stdio.texi: Add description of printf_size and printf_size_info. Partly based on the documentation by Larry McVoy. 1997-04-19 02:21 Ulrich Drepper <drepper@cygnus.com> * stdio-common/printf_size.c (printf_size): Correct values for `units'. Report by Larry McVoy <lm@neteng.engr.sgi.com>. * stdio-common/tst-printfsz.c: New file. * stdio-common/Makefile (tests): Add tst-printfsz.c. (CFLAGS-tst-printfsz.c): Define to prevent warnings about format strings. 1997-04-18 15:48 Ulrich Drepper <drepper@cygnus.com> * login/utmp.h: Add prototype for updwtmp. * login/logwtmp.c: Add new function updwtmp which allows to write a complete record to the wtmp file. Patch by Miquel van Smoorenburg <miquels@cistron.nl>. 1997-04-17 17:57 Andreas Schwab <schwab@issan.informatik.uni-dortmund.de> * math/Makefile (headers): Add mathbits.h. 1997-04-16 21:20 Andreas Schwab <schwab@issan.informatik.uni-dortmund.de> * sysdeps/m68k/fpu/__math.h: Add inlined sincos{,l,f}. * sysdeps/m68k/fpu/s_sincos.c: New file. * sysdeps/m68k/fpu/s_sincosf.c: New file. * sysdeps/m68k/fpu/s_sincosl.c: New file. * sysdeps/libm-ieee754/e_scalb.c: Use internal names of the functions. * sysdeps/libm-ieee754/e_scalbl.c: Likewise. * sysdeps/libm-ieee754/s_ctanh.c: Use isfinite instead of finite. * sysdeps/libm-ieee754/s_ctanhf.c: Likewise. * sysdeps/libm-ieee754/s_ctanhl.c: Likewise. * sysdeps/libm-ieee754/s_ctan.c: Likewise. * sysdeps/libm-ieee754/s_ctanf.c: Likewise. * sysdeps/libm-ieee754/s_ctanl.c: Likewise. Fix type of `res'. 1997-04-18 11:21 Ulrich Drepper <drepper@cygnus.com> * shadow/fgetspent_r.c: Set *RESULT to NULL before returning error. Patch by Thorsten Kukuk <kukuk@vt.uni-paderborn.de>.
78 lines
1.7 KiB
C
78 lines
1.7 KiB
C
/* s_rintf.c -- float version of s_rint.c.
|
|
* Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com.
|
|
*/
|
|
/* Adapted for use as nearbyint by Ulrich Drepper <drepper@cygnus.com>. */
|
|
|
|
/*
|
|
* ====================================================
|
|
* Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
|
|
*
|
|
* Developed at SunPro, a Sun Microsystems, Inc. business.
|
|
* Permission to use, copy, modify, and distribute this
|
|
* software is freely granted, provided that this notice
|
|
* is preserved.
|
|
* ====================================================
|
|
*/
|
|
|
|
|
|
#include <fenv.h>
|
|
#include "math.h"
|
|
#include "math_private.h"
|
|
|
|
#ifdef __STDC__
|
|
static const float
|
|
#else
|
|
static float
|
|
#endif
|
|
TWO23[2]={
|
|
8.3886080000e+06, /* 0x4b000000 */
|
|
-8.3886080000e+06, /* 0xcb000000 */
|
|
};
|
|
|
|
#ifdef __STDC__
|
|
float __nearbyintf(float x)
|
|
#else
|
|
float __nearbyintf(x)
|
|
float x;
|
|
#endif
|
|
{
|
|
fenv_t env;
|
|
int32_t i0,j0,sx;
|
|
u_int32_t i,i1;
|
|
float w,t;
|
|
GET_FLOAT_WORD(i0,x);
|
|
sx = (i0>>31)&1;
|
|
j0 = ((i0>>23)&0xff)-0x7f;
|
|
if(j0<23) {
|
|
if(j0<0) {
|
|
if((i0&0x7fffffff)==0) return x;
|
|
i1 = (i0&0x07fffff);
|
|
i0 &= 0xfff00000;
|
|
i0 |= ((i1|-i1)>>9)&0x400000;
|
|
SET_FLOAT_WORD(x,i0);
|
|
feholdexcept (&env);
|
|
w = TWO23[sx]+x;
|
|
t = w-TWO23[sx];
|
|
fesetenv (&env);
|
|
GET_FLOAT_WORD(i0,t);
|
|
SET_FLOAT_WORD(t,(i0&0x7fffffff)|(sx<<31));
|
|
return t;
|
|
} else {
|
|
i = (0x007fffff)>>j0;
|
|
if((i0&i)==0) return x; /* x is integral */
|
|
i>>=1;
|
|
if((i0&i)!=0) i0 = (i0&(~i))|((0x100000)>>j0);
|
|
}
|
|
} else {
|
|
if(j0==0x80) return x+x; /* inf or NaN */
|
|
else return x; /* x is integral */
|
|
}
|
|
SET_FLOAT_WORD(x,i0);
|
|
feholdexcept (&env);
|
|
w = TWO23[sx]+x;
|
|
t = w-TWO23[sx];
|
|
fesetenv (&env);
|
|
return t;
|
|
}
|
|
weak_alias (__nearbyintf, nearbyintf)
|