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>.
101 lines
2.5 KiB
C
101 lines
2.5 KiB
C
/* s_rintl.c -- long double version of s_rint.c.
|
|
* Conversion to long double by Ulrich Drepper,
|
|
* Cygnus Support, drepper@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.
|
|
* ====================================================
|
|
*/
|
|
|
|
/*
|
|
* rintl(x)
|
|
* Return x rounded to integral value according to the prevailing
|
|
* rounding mode.
|
|
* Method:
|
|
* Using floating addition.
|
|
* Exception:
|
|
* Inexact flag raised if x not equal to rintl(x).
|
|
*/
|
|
|
|
#include <fenv.h>
|
|
#include "math.h"
|
|
#include "math_private.h"
|
|
|
|
#ifdef __STDC__
|
|
static const long double
|
|
#else
|
|
static long double
|
|
#endif
|
|
TWO63[2]={
|
|
9.223372036854775808000000e+18, /* 0x403E, 0x00000000, 0x00000000 */
|
|
-9.223372036854775808000000e+18 /* 0xC03E, 0x00000000, 0x00000000 */
|
|
};
|
|
|
|
#ifdef __STDC__
|
|
long double __nearbyintl(long double x)
|
|
#else
|
|
long double __nearbyintl(x)
|
|
long double x;
|
|
#endif
|
|
{
|
|
fenv_t env;
|
|
int32_t se,j0,sx;
|
|
u_int32_t i,i0,i1;
|
|
long double w,t;
|
|
GET_LDOUBLE_WORDS(se,i0,i1,x);
|
|
sx = (se>>15)&1;
|
|
j0 = (se&0x7fff)-0x3fff;
|
|
if(j0<31) {
|
|
if(j0<0) {
|
|
if(((se&0x7fff)|i0|i1)==0) return x;
|
|
i1 |= i0;
|
|
i0 &= 0xe0000000;
|
|
i0 |= (i1|-i1)&0x80000000;
|
|
SET_LDOUBLE_MSW(x,i0);
|
|
feholdexcept (&env);
|
|
w = TWO63[sx]+x;
|
|
t = w-TWO63[sx];
|
|
fesetenv (&env);
|
|
GET_LDOUBLE_EXP(i0,t);
|
|
SET_LDOUBLE_EXP(t,(i0&0x7fff)|(sx<<15));
|
|
return t;
|
|
} else {
|
|
i = (0x7fffffff)>>j0;
|
|
if(((i0&i)|i1)==0) return x; /* x is integral */
|
|
i>>=1;
|
|
if(((i0&i)|i1)!=0) {
|
|
if(j0==31) i1 = 0x40000000; else
|
|
i0 = (i0&(~i))|((0x20000000)>>j0);
|
|
/* Shouldn't this be
|
|
if (j0 >= 30) i1 = 0x80000000 >> (j0 - 30);
|
|
i0 = (i0&(~i))|((0x20000000)>>j0);
|
|
If yes, this should be correct in s_rint and
|
|
s_rintf, too. -- drepper@cygnus.com */
|
|
}
|
|
}
|
|
} else if (j0>62) {
|
|
if(j0==0x4000) return x+x; /* inf or NaN */
|
|
else return x; /* x is integral */
|
|
} else {
|
|
i = ((u_int32_t)(0xffffffff))>>(j0-31);
|
|
if((i1&i)==0) return x; /* x is integral */
|
|
i>>=1;
|
|
if((i1&i)!=0) i1 = (i1&(~i))|((0x40000000)>>(j0-31));
|
|
}
|
|
SET_LDOUBLE_WORDS(x,se,i0,i1);
|
|
feholdexcept (&env);
|
|
w = TWO63[sx]+x;
|
|
t = w-TWO63[sx];
|
|
fesetenv (&env);
|
|
return t;
|
|
}
|
|
weak_alias (__nearbyintl, nearbyintl)
|