+2017-08-03 Joseph Myers <joseph@codesourcery.com>
+
+ * math/s_nextafter.c (__nextafter): Use uintN_t instead of
+ u_intN_t.
+ * math/s_nexttowardf.c (__nexttowardf): Likewise.
+ * sysdeps/generic/math_private.h (ieee_double_shape_type):
+ Likewise.
+ (ieee_float_shape_type): Likewise.
+ * sysdeps/i386/fpu/s_fpclassifyl.c (__fpclassifyl): Likewise.
+ * sysdeps/i386/fpu/s_isnanl.c (__isnanl): Likewise.
+ * sysdeps/i386/fpu/s_nextafterl.c (__nextafterl): Likewise.
+ * sysdeps/i386/fpu/s_nexttoward.c (__nexttoward): Likewise.
+ * sysdeps/i386/fpu/s_nexttowardf.c (__nexttowardf): Likewise.
+ * sysdeps/ieee754/dbl-64/e_acosh.c (__ieee754_acosh): Likewise.
+ * sysdeps/ieee754/dbl-64/e_cosh.c (__ieee754_cosh): Likewise.
+ * sysdeps/ieee754/dbl-64/e_fmod.c (__ieee754_fmod): Likewise.
+ * sysdeps/ieee754/dbl-64/e_gamma_r.c (__ieee754_gamma_r):
+ Likewise.
+ * sysdeps/ieee754/dbl-64/e_hypot.c (__ieee754_hypot): Likewise.
+ * sysdeps/ieee754/dbl-64/e_jn.c (__ieee754_jn): Likewise.
+ (__ieee754_yn): Likewise.
+ * sysdeps/ieee754/dbl-64/e_log10.c (__ieee754_log10): Likewise.
+ * sysdeps/ieee754/dbl-64/e_log2.c (__ieee754_log2): Likewise.
+ * sysdeps/ieee754/dbl-64/e_rem_pio2.c (__ieee754_rem_pio2):
+ Likewise.
+ * sysdeps/ieee754/dbl-64/e_sinh.c (__ieee754_sinh): Likewise.
+ * sysdeps/ieee754/dbl-64/s_ceil.c (__ceil): Likewise.
+ * sysdeps/ieee754/dbl-64/s_copysign.c (__copysign): Likewise.
+ * sysdeps/ieee754/dbl-64/s_erf.c (__erf): Likewise.
+ (__erfc): Likewise.
+ * sysdeps/ieee754/dbl-64/s_expm1.c (__expm1): Likewise.
+ * sysdeps/ieee754/dbl-64/s_finite.c (FINITE): Likewise.
+ * sysdeps/ieee754/dbl-64/s_floor.c (__floor): Likewise.
+ * sysdeps/ieee754/dbl-64/s_fpclassify.c (__fpclassify): Likewise.
+ * sysdeps/ieee754/dbl-64/s_isnan.c (__isnan): Likewise.
+ * sysdeps/ieee754/dbl-64/s_issignaling.c (__issignaling):
+ Likewise.
+ * sysdeps/ieee754/dbl-64/s_llrint.c (__llrint): Likewise.
+ * sysdeps/ieee754/dbl-64/s_llround.c (__llround): Likewise.
+ * sysdeps/ieee754/dbl-64/s_lrint.c (__lrint): Likewise.
+ * sysdeps/ieee754/dbl-64/s_lround.c (__lround): Likewise.
+ * sysdeps/ieee754/dbl-64/s_modf.c (__modf): Likewise.
+ * sysdeps/ieee754/dbl-64/s_nextup.c (__nextup): Likewise.
+ * sysdeps/ieee754/dbl-64/s_remquo.c (__remquo): Likewise.
+ * sysdeps/ieee754/dbl-64/s_round.c (__round): Likewise.
+ * sysdeps/ieee754/dbl-64/s_trunc.c (__trunc): Likewise.
+ * sysdeps/ieee754/dbl-64/wordsize-64/s_issignaling.c
+ (__issignaling): Likewise.
+ * sysdeps/ieee754/flt-32/e_atan2f.c (__ieee754_atan2f): Likewise.
+ * sysdeps/ieee754/flt-32/e_fmodf.c (__ieee754_fmodf): Likewise.
+ * sysdeps/ieee754/flt-32/e_gammaf_r.c (__ieee754_gammaf_r):
+ Likewise.
+ * sysdeps/ieee754/flt-32/e_jnf.c (__ieee754_ynf): Likewise.
+ * sysdeps/ieee754/flt-32/e_log10f.c (__ieee754_log10f): Likewise.
+ * sysdeps/ieee754/flt-32/e_powf.c (__ieee754_powf): Likewise.
+ * sysdeps/ieee754/flt-32/e_rem_pio2f.c (__ieee754_rem_pio2f):
+ Likewise.
+ * sysdeps/ieee754/flt-32/e_remainderf.c (__ieee754_remainderf):
+ Likewise.
+ * sysdeps/ieee754/flt-32/e_sqrtf.c (__ieee754_sqrtf): Likewise.
+ * sysdeps/ieee754/flt-32/s_ceilf.c (__ceilf): Likewise.
+ * sysdeps/ieee754/flt-32/s_copysignf.c (__copysignf): Likewise.
+ * sysdeps/ieee754/flt-32/s_erff.c (__erff): Likewise.
+ (__erfcf): Likewise.
+ * sysdeps/ieee754/flt-32/s_expm1f.c (__expm1f): Likewise.
+ * sysdeps/ieee754/flt-32/s_finitef.c (FINITEF): Likewise.
+ * sysdeps/ieee754/flt-32/s_floorf.c (__floorf): Likewise.
+ * sysdeps/ieee754/flt-32/s_fpclassifyf.c (__fpclassifyf):
+ Likewise.
+ * sysdeps/ieee754/flt-32/s_isnanf.c (__isnanf): Likewise.
+ * sysdeps/ieee754/flt-32/s_issignalingf.c (__issignalingf):
+ Likewise.
+ * sysdeps/ieee754/flt-32/s_llrintf.c (__llrintf): Likewise.
+ * sysdeps/ieee754/flt-32/s_llroundf.c (__llroundf): Likewise.
+ * sysdeps/ieee754/flt-32/s_lrintf.c (__lrintf): Likewise.
+ * sysdeps/ieee754/flt-32/s_lroundf.c (__lroundf): Likewise.
+ * sysdeps/ieee754/flt-32/s_modff.c (__modff): Likewise.
+ * sysdeps/ieee754/flt-32/s_remquof.c (__remquof): Likewise.
+ * sysdeps/ieee754/flt-32/s_roundf.c (__roundf): Likewise.
+ * sysdeps/ieee754/ldbl-128/e_acoshl.c (__ieee754_acoshl):
+ Likewise.
+ * sysdeps/ieee754/ldbl-128/e_atan2l.c (__ieee754_atan2l):
+ Likewise.
+ * sysdeps/ieee754/ldbl-128/e_atanhl.c (__ieee754_atanhl):
+ Likewise.
+ * sysdeps/ieee754/ldbl-128/e_fmodl.c (__ieee754_fmodl): Likewise.
+ * sysdeps/ieee754/ldbl-128/e_gammal_r.c (__ieee754_gammal_r):
+ Likewise.
+ * sysdeps/ieee754/ldbl-128/e_hypotl.c (__ieee754_hypotl):
+ Likewise.
+ * sysdeps/ieee754/ldbl-128/e_jnl.c (__ieee754_jnl): Likewise.
+ (__ieee754_ynl): Likewise.
+ * sysdeps/ieee754/ldbl-128/e_powl.c (__ieee754_powl): Likewise.
+ * sysdeps/ieee754/ldbl-128/e_rem_pio2l.c (__ieee754_rem_pio2l):
+ Likewise.
+ * sysdeps/ieee754/ldbl-128/e_remainderl.c (__ieee754_remainderl):
+ Likewise.
+ * sysdeps/ieee754/ldbl-128/e_sinhl.c (__ieee754_sinhl): Likewise.
+ * sysdeps/ieee754/ldbl-128/k_cosl.c (__kernel_cosl): Likewise.
+ * sysdeps/ieee754/ldbl-128/k_sincosl.c (__kernel_sincosl):
+ Likewise.
+ * sysdeps/ieee754/ldbl-128/k_sinl.c (__kernel_sinl): Likewise.
+ * sysdeps/ieee754/ldbl-128/s_ceill.c (__ceill): Likewise.
+ * sysdeps/ieee754/ldbl-128/s_copysignl.c (__copysignl): Likewise.
+ * sysdeps/ieee754/ldbl-128/s_erfl.c (__erfcl): Likewise.
+ * sysdeps/ieee754/ldbl-128/s_fabsl.c (__fabsl): Likewise.
+ * sysdeps/ieee754/ldbl-128/s_finitel.c (__finitel): Likewise.
+ * sysdeps/ieee754/ldbl-128/s_floorl.c (__floorl): Likewise.
+ * sysdeps/ieee754/ldbl-128/s_fpclassifyl.c (__fpclassifyl):
+ Likewise.
+ * sysdeps/ieee754/ldbl-128/s_frexpl.c (__frexpl): Likewise.
+ * sysdeps/ieee754/ldbl-128/s_isnanl.c (__isnanl): Likewise.
+ * sysdeps/ieee754/ldbl-128/s_issignalingl.c (__issignalingl):
+ Likewise.
+ * sysdeps/ieee754/ldbl-128/s_llrintl.c (__llrintl): Likewise.
+ * sysdeps/ieee754/ldbl-128/s_llroundl.c (__llroundl): Likewise.
+ * sysdeps/ieee754/ldbl-128/s_lrintl.c (__lrintl): Likewise.
+ * sysdeps/ieee754/ldbl-128/s_lroundl.c (__lroundl): Likewise.
+ * sysdeps/ieee754/ldbl-128/s_modfl.c (__modfl): Likewise.
+ * sysdeps/ieee754/ldbl-128/s_nearbyintl.c (__nearbyintl):
+ Likewise.
+ * sysdeps/ieee754/ldbl-128/s_nextafterl.c (__nextafterl):
+ Likewise.
+ * sysdeps/ieee754/ldbl-128/s_nexttoward.c (__nexttoward):
+ Likewise.
+ * sysdeps/ieee754/ldbl-128/s_nexttowardf.c (__nexttowardf):
+ Likewise.
+ * sysdeps/ieee754/ldbl-128/s_nextupl.c (__nextupl): Likewise.
+ * sysdeps/ieee754/ldbl-128/s_remquol.c (__remquol): Likewise.
+ * sysdeps/ieee754/ldbl-128/s_rintl.c (__rintl): Likewise.
+ * sysdeps/ieee754/ldbl-128/s_roundl.c (__roundl): Likewise.
+ * sysdeps/ieee754/ldbl-128/s_tanhl.c (__tanhl): Likewise.
+ * sysdeps/ieee754/ldbl-128/s_truncl.c (__truncl): Likewise.
+ * sysdeps/ieee754/ldbl-128ibm/e_fmodl.c (__ieee754_fmodl):
+ Likewise.
+ * sysdeps/ieee754/ldbl-128ibm/e_gammal_r.c (__ieee754_gammal_r):
+ Likewise.
+ * sysdeps/ieee754/ldbl-128ibm/e_powl.c (__ieee754_powl): Likewise.
+ * sysdeps/ieee754/ldbl-128ibm/e_rem_pio2l.c (__ieee754_rem_pio2l):
+ Likewise.
+ * sysdeps/ieee754/ldbl-128ibm/e_remainderl.c
+ (__ieee754_remainderl): Likewise.
+ * sysdeps/ieee754/ldbl-128ibm/k_cosl.c (__kernel_cosl): Likewise.
+ * sysdeps/ieee754/ldbl-128ibm/k_sinl.c (__kernel_sinl): Likewise.
+ * sysdeps/ieee754/ldbl-128ibm/s_fabsl.c (__fabsl): Likewise.
+ * sysdeps/ieee754/ldbl-128ibm/s_fpclassifyl.c (___fpclassifyl):
+ Likewise.
+ * sysdeps/ieee754/ldbl-128ibm/s_modfl.c (__modfl): Likewise.
+ * sysdeps/ieee754/ldbl-128ibm/s_nexttowardf.c (__nexttowardf):
+ Likewise.
+ * sysdeps/ieee754/ldbl-128ibm/s_remquol.c (__remquol): Likewise.
+ * sysdeps/ieee754/ldbl-96/e_acoshl.c (__ieee754_acoshl): Likewise.
+ * sysdeps/ieee754/ldbl-96/e_asinl.c (__ieee754_asinl): Likewise.
+ * sysdeps/ieee754/ldbl-96/e_atanhl.c (__ieee754_atanhl): Likewise.
+ * sysdeps/ieee754/ldbl-96/e_coshl.c (__ieee754_coshl): Likewise.
+ * sysdeps/ieee754/ldbl-96/e_gammal_r.c (__ieee754_gammal_r):
+ Likewise.
+ * sysdeps/ieee754/ldbl-96/e_hypotl.c (__ieee754_hypotl): Likewise.
+ * sysdeps/ieee754/ldbl-96/e_j0l.c (__ieee754_j0l): Likewise.
+ (__ieee754_y0l): Likewise.
+ (pzero): Likewise.
+ (qzero): Likewise.
+ * sysdeps/ieee754/ldbl-96/e_j1l.c (__ieee754_j1l): Likewise.
+ (__ieee754_y1l): Likewise.
+ (pone): Likewise.
+ (qone): Likewise.
+ * sysdeps/ieee754/ldbl-96/e_jnl.c (__ieee754_jnl): Likewise.
+ (__ieee754_ynl): Likewise.
+ * sysdeps/ieee754/ldbl-96/e_lgammal_r.c (sin_pi): Likewise.
+ (__ieee754_lgammal_r): Likewise.
+ * sysdeps/ieee754/ldbl-96/e_rem_pio2l.c (__ieee754_rem_pio2l):
+ Likewise.
+ * sysdeps/ieee754/ldbl-96/e_sinhl.c (__ieee754_sinhl): Likewise.
+ * sysdeps/ieee754/ldbl-96/s_copysignl.c (__copysignl): Likewise.
+ * sysdeps/ieee754/ldbl-96/s_erfl.c (__erfl): Likewise.
+ (__erfcl): Likewise.
+ * sysdeps/ieee754/ldbl-96/s_frexpl.c (__frexpl): Likewise.
+ * sysdeps/ieee754/ldbl-96/s_issignalingl.c (__issignalingl):
+ Likewise.
+ * sysdeps/ieee754/ldbl-96/s_llrintl.c (__llrintl): Likewise.
+ * sysdeps/ieee754/ldbl-96/s_llroundl.c (__llroundl): Likewise.
+ * sysdeps/ieee754/ldbl-96/s_lrintl.c (__lrintl): Likewise.
+ * sysdeps/ieee754/ldbl-96/s_lroundl.c (__lroundl): Likewise.
+ * sysdeps/ieee754/ldbl-96/s_modfl.c (__modfl): Likewise.
+ * sysdeps/ieee754/ldbl-96/s_nexttoward.c (__nexttoward): Likewise.
+ * sysdeps/ieee754/ldbl-96/s_nexttowardf.c (__nexttowardf):
+ Likewise.
+ * sysdeps/ieee754/ldbl-96/s_nextupl.c (__nextupl): Likewise.
+ * sysdeps/ieee754/ldbl-96/s_remquol.c (__remquol): Likewise.
+ * sysdeps/ieee754/ldbl-96/s_roundl.c (__roundl): Likewise.
+ * sysdeps/ieee754/ldbl-96/s_tanhl.c (__tanhl): Likewise.
+ * sysdeps/ieee754/ldbl-opt/s_nexttowardfd.c (__nldbl_nexttowardf):
+ Likewise.
+ * sysdeps/m68k/m680x0/fpu/e_pow.c (s(__ieee754_pow)): Likewise.
+ * sysdeps/m68k/m680x0/fpu/s_fpclassifyl.c (__fpclassifyl):
+ Likewise.
+ * sysdeps/m68k/m680x0/fpu/s_llrint.c (__llrint): Likewise.
+ * sysdeps/m68k/m680x0/fpu/s_llrintf.c (__llrintf): Likewise.
+ * sysdeps/m68k/m680x0/fpu/s_llrintl.c (__llrintl): Likewise.
+ * sysdeps/m68k/m680x0/fpu/s_nextafterl.c (__nextafterl): Likewise.
+ * sysdeps/x86/fpu/powl_helper.c (__powl_helper): Likewise.
+
2017-08-03 Florian Weimer <fweimer@redhat.com>
[BZ #21885]
double __nextafter(double x, double y)
{
int32_t hx,hy,ix,iy;
- u_int32_t lx,ly;
+ uint32_t lx,ly;
EXTRACT_WORDS(hx,lx,x);
EXTRACT_WORDS(hy,ly,y);
float __nexttowardf(float x, long double y)
{
int32_t hx,hy,ix,iy;
- u_int32_t ly;
+ uint32_t ly;
GET_FLOAT_WORD(hx,x);
EXTRACT_WORDS(hy,ly,y);
if((long double) x==y) return y; /* x=y, return y */
if(ix==0) { /* x == 0 */
float u;
- SET_FLOAT_WORD(x,(u_int32_t)(hy&0x80000000)|1);/* return +-minsub*/
+ SET_FLOAT_WORD(x,(uint32_t)(hy&0x80000000)|1);/* return +-minsub*/
u = math_opt_barrier (x);
u = u * u;
math_force_eval (u); /* raise underflow flag */
double value;
struct
{
- u_int32_t msw;
- u_int32_t lsw;
+ uint32_t msw;
+ uint32_t lsw;
} parts;
uint64_t word;
} ieee_double_shape_type;
double value;
struct
{
- u_int32_t lsw;
- u_int32_t msw;
+ uint32_t lsw;
+ uint32_t msw;
} parts;
uint64_t word;
} ieee_double_shape_type;
typedef union
{
float value;
- u_int32_t word;
+ uint32_t word;
} ieee_float_shape_type;
/* Get a 32 bit int from a float. */
int
__fpclassifyl (long double x)
{
- u_int32_t ex, hx, lx;
+ uint32_t ex, hx, lx;
int retval = FP_NORMAL;
GET_LDOUBLE_WORDS (ex, hx, lx, x);
extended format has the normally implicit 1 explicit
present. Sigh! */
lx |= hx & 0x7fffffff;
- se |= (u_int32_t)(lx|(-lx))>>31;
+ se |= (uint32_t)(lx|(-lx))>>31;
se = 0xfffe - se;
- return (int)((u_int32_t)(se))>>16;
+ return (int)((uint32_t)(se))>>16;
}
hidden_def (__isnanl)
weak_alias (__isnanl, isnanl)
long double __nextafterl(long double x, long double y)
{
- u_int32_t hx,hy,ix,iy;
- u_int32_t lx,ly;
+ uint32_t hx,hy,ix,iy;
+ uint32_t lx,ly;
int32_t esx,esy;
GET_LDOUBLE_WORDS(esx,hx,lx,x);
double __nexttoward(double x, long double y)
{
int32_t hx,ix,iy;
- u_int32_t lx,hy,ly,esy;
+ uint32_t lx,hy,ly,esy;
EXTRACT_WORDS(hx,lx,x);
GET_LDOUBLE_WORDS(esy,hy,ly,y);
float __nexttowardf(float x, long double y)
{
int32_t hx,ix,iy;
- u_int32_t hy,ly,esy;
+ uint32_t hy,ly,esy;
GET_FLOAT_WORD(hx,x);
GET_LDOUBLE_WORDS(esy,hy,ly,y);
{
double t;
int32_t hx;
- u_int32_t lx;
+ uint32_t lx;
EXTRACT_WORDS (hx, lx, x);
if (hx < 0x3ff00000) /* x < 1 */
{
{
double t, w;
int32_t ix;
- u_int32_t lx;
+ uint32_t lx;
/* High word of |x|. */
GET_HIGH_WORD (ix, x);
/* |x| in [log(maxdouble), overflowthresold] */
GET_LOW_WORD (lx, x);
- if (ix < 0x408633ce || ((ix == 0x408633ce) && (lx <= (u_int32_t) 0x8fb9f87d)))
+ if (ix < 0x408633ce || ((ix == 0x408633ce) && (lx <= (uint32_t) 0x8fb9f87d)))
{
w = __ieee754_exp (half * fabs (x));
t = half * w;
__ieee754_fmod (double x, double y)
{
int32_t n, hx, hy, hz, ix, iy, sx, i;
- u_int32_t lx, ly, lz;
+ uint32_t lx, ly, lz;
EXTRACT_WORDS (hx, lx, x);
EXTRACT_WORDS (hy, ly, y);
if ((hx < hy) || (lx < ly))
return x; /* |x|<|y| return x */
if (lx == ly)
- return Zero[(u_int32_t) sx >> 31]; /* |x|=|y| return x*0*/
+ return Zero[(uint32_t) sx >> 31]; /* |x|=|y| return x*0*/
}
/* determine ix = ilogb(x) */
else
{
if ((hz | lz) == 0) /* return sign(x)*0 */
- return Zero[(u_int32_t) sx >> 31];
+ return Zero[(uint32_t) sx >> 31];
hx = hz + hz + (lz >> 31); lx = lz + lz;
}
}
/* convert back to floating value and restore the sign */
if ((hx | lx) == 0) /* return sign(x)*0 */
- return Zero[(u_int32_t) sx >> 31];
+ return Zero[(uint32_t) sx >> 31];
while (hx < 0x00100000) /* normalize x */
{
hx = hx + hx + (lx >> 31); lx = lx + lx;
n = -1022 - iy;
if (n <= 20)
{
- lx = (lx >> n) | ((u_int32_t) hx << (32 - n));
+ lx = (lx >> n) | ((uint32_t) hx << (32 - n));
hx >>= n;
}
else if (n <= 31)
__ieee754_gamma_r (double x, int *signgamp)
{
int32_t hx;
- u_int32_t lx;
+ uint32_t lx;
double ret;
EXTRACT_WORDS (hx, lx, x);
return 1.0 / x;
}
if (__builtin_expect (hx < 0, 0)
- && (u_int32_t) hx < 0xfff00000 && __rint (x) == x)
+ && (uint32_t) hx < 0xfff00000 && __rint (x) == x)
{
/* Return value for integer x < 0 is NaN with invalid exception. */
*signgamp = 0;
{
if (ha >= 0x7ff00000) /* Inf or NaN */
{
- u_int32_t low;
+ uint32_t low;
w = a + b; /* for sNaN */
if (issignaling (a) || issignaling (b))
return w;
{
if (hb <= 0x000fffff) /* subnormal b or 0 */
{
- u_int32_t low;
+ uint32_t low;
GET_LOW_WORD (low, b);
if ((hb | low) == 0)
return a;
}
if (k != 0)
{
- u_int32_t high;
+ uint32_t high;
t1 = 1.0;
GET_HIGH_WORD (high, t1);
SET_HIGH_WORD (t1, high + (k << 20));
EXTRACT_WORDS (hx, lx, x);
ix = 0x7fffffff & hx;
/* if J(n,NaN) is NaN */
- if (__glibc_unlikely ((ix | ((u_int32_t) (lx | -lx)) >> 31) > 0x7ff00000))
+ if (__glibc_unlikely ((ix | ((uint32_t) (lx | -lx)) >> 31) > 0x7ff00000))
return x + x;
if (n < 0)
{
EXTRACT_WORDS (hx, lx, x);
ix = 0x7fffffff & hx;
/* if Y(n,NaN) is NaN */
- if (__glibc_unlikely ((ix | ((u_int32_t) (lx | -lx)) >> 31) > 0x7ff00000))
+ if (__glibc_unlikely ((ix | ((uint32_t) (lx | -lx)) >> 31) > 0x7ff00000))
return x + x;
if (__glibc_unlikely ((ix | lx) == 0))
return -HUGE_VAL + x;
}
else
{
- u_int32_t high;
+ uint32_t high;
a = __ieee754_y0 (x);
b = __ieee754_y1 (x);
/* quit if b is -inf */
{
double y, z;
int32_t i, k, hx;
- u_int32_t lx;
+ uint32_t lx;
EXTRACT_WORDS (hx, lx, x);
if (__glibc_unlikely (hx >= 0x7ff00000))
return x + x;
k += (hx >> 20) - 1023;
- i = ((u_int32_t) k & 0x80000000) >> 31;
+ i = ((uint32_t) k & 0x80000000) >> 31;
hx = (hx & 0x000fffff) | ((0x3ff - i) << 20);
y = (double) (k + i);
if (FIX_INT_FP_CONVERT_ZERO && y == 0.0)
{
double hfsq, f, s, z, R, w, t1, t2, dk;
int32_t k, hx, i, j;
- u_int32_t lx;
+ uint32_t lx;
EXTRACT_WORDS (hx, lx, x);
double z, w, t, r, fn;
double tx[3];
int32_t e0, i, j, nx, n, ix, hx;
- u_int32_t low;
+ uint32_t low;
GET_HIGH_WORD (hx, x); /* high word of x */
ix = hx & 0x7fffffff;
}
else
{
- u_int32_t high;
+ uint32_t high;
j = ix >> 20;
y[0] = r - w;
GET_HIGH_WORD (high, y[0]);
{
double t, w, h;
int32_t ix, jx;
- u_int32_t lx;
+ uint32_t lx;
/* High word of |x|. */
GET_HIGH_WORD (jx, x);
/* |x| in [log(maxdouble), overflowthresold] */
GET_LOW_WORD (lx, x);
- if (ix < 0x408633ce || ((ix == 0x408633ce) && (lx <= (u_int32_t) 0x8fb9f87d)))
+ if (ix < 0x408633ce || ((ix == 0x408633ce) && (lx <= (uint32_t) 0x8fb9f87d)))
{
w = __ieee754_exp (0.5 * fabs (x));
t = h * w;
__ceil (double x)
{
int32_t i0, i1, j0;
- u_int32_t i, j;
+ uint32_t i, j;
EXTRACT_WORDS (i0, i1, x);
j0 = ((i0 >> 20) & 0x7ff) - 0x3ff;
if (j0 < 20)
}
else
{
- i = ((u_int32_t) (0xffffffff)) >> (j0 - 20);
+ i = ((uint32_t) (0xffffffff)) >> (j0 - 20);
if ((i1 & i) == 0)
return x; /* x is integral */
if (i0 > 0)
double
__copysign (double x, double y)
{
- u_int32_t hx, hy;
+ uint32_t hx, hy;
GET_HIGH_WORD (hx, x);
GET_HIGH_WORD (hy, y);
SET_HIGH_WORD (x, (hx & 0x7fffffff) | (hy & 0x80000000));
ix = hx & 0x7fffffff;
if (ix >= 0x7ff00000) /* erf(nan)=nan */
{
- i = ((u_int32_t) hx >> 31) << 1;
+ i = ((uint32_t) hx >> 31) << 1;
return (double) (1 - i) + one / x; /* erf(+-inf)=+-1 */
}
ix = hx & 0x7fffffff;
if (ix >= 0x7ff00000) /* erfc(nan)=nan */
{ /* erfc(+-inf)=0,2 */
- double ret = (double) (((u_int32_t) hx >> 31) << 1) + one / x;
+ double ret = (double) (((uint32_t) hx >> 31) << 1) + one / x;
if (FIX_INT_FP_CONVERT_ZERO && ret == 0.0)
return 0.0;
return ret;
{
double y, hi, lo, c, t, e, hxs, hfx, r1, h2, h4, R1, R2, R3;
int32_t k, xsb;
- u_int32_t hx;
+ uint32_t hx;
GET_HIGH_WORD (hx, x);
xsb = hx & 0x80000000; /* sign bit of x */
{
if (hx >= 0x7ff00000)
{
- u_int32_t low;
+ uint32_t low;
GET_LOW_WORD (low, x);
if (((hx & 0xfffff) | low) != 0)
return x + x; /* NaN */
}
if (k <= -2 || k > 56) /* suffice to return exp(x)-1 */
{
- u_int32_t high;
+ uint32_t high;
y = one - (e - x);
GET_HIGH_WORD (high, y);
SET_HIGH_WORD (y, high + (k << 20)); /* add k to y's exponent */
t = one;
if (k < 20)
{
- u_int32_t high;
+ uint32_t high;
SET_HIGH_WORD (t, 0x3ff00000 - (0x200000 >> k)); /* t=1-2^-k */
y = t - (e - x);
GET_HIGH_WORD (high, y);
}
else
{
- u_int32_t high;
+ uint32_t high;
SET_HIGH_WORD (t, ((0x3ff - k) << 20)); /* 2^-k */
y = x - (e + t);
y += one;
{
int32_t hx;
GET_HIGH_WORD (hx, x);
- return (int) ((u_int32_t) ((hx & 0x7ff00000) - 0x7ff00000) >> 31);
+ return (int) ((uint32_t) ((hx & 0x7ff00000) - 0x7ff00000) >> 31);
}
hidden_def (__finite)
weak_alias (__finite, finite)
__floor (double x)
{
int32_t i0, i1, j0;
- u_int32_t i, j;
+ uint32_t i, j;
EXTRACT_WORDS (i0, i1, x);
j0 = ((i0 >> 20) & 0x7ff) - 0x3ff;
if (j0 < 20)
}
else
{
- i = ((u_int32_t) (0xffffffff)) >> (j0 - 20);
+ i = ((uint32_t) (0xffffffff)) >> (j0 - 20);
if ((i1 & i) == 0)
return x; /* x is integral */
if (i0 < 0)
int
__fpclassify (double x)
{
- u_int32_t hx, lx;
+ uint32_t hx, lx;
int retval = FP_NORMAL;
EXTRACT_WORDS (hx, lx, x);
int32_t hx, lx;
EXTRACT_WORDS (hx, lx, x);
hx &= 0x7fffffff;
- hx |= (u_int32_t) (lx | (-lx)) >> 31;
+ hx |= (uint32_t) (lx | (-lx)) >> 31;
hx = 0x7ff00000 - hx;
- return (int) (((u_int32_t) hx) >> 31);
+ return (int) (((uint32_t) hx) >> 31);
}
hidden_def (__isnan)
weak_alias (__isnan, isnan)
__issignaling (double x)
{
#if HIGH_ORDER_BIT_IS_SET_FOR_SNAN
- u_int32_t hxi;
+ uint32_t hxi;
GET_HIGH_WORD (hxi, x);
/* We only have to care about the high-order bit of x's significand, because
having it set (sNaN) already makes the significand different from that
used to designate infinity. */
return (hxi & 0x7ff80000) == 0x7ff80000;
#else
- u_int32_t hxi, lxi;
+ uint32_t hxi, lxi;
EXTRACT_WORDS (hxi, lxi, x);
/* To keep the following comparison simple, toggle the quiet/signaling bit,
so that it is set for sNaNs. This is inverse to IEEE 754-2008 (as well as
__llrint (double x)
{
int32_t j0;
- u_int32_t i1, i0;
+ uint32_t i1, i0;
long long int result;
double w;
double t;
__llround (double x)
{
int32_t j0;
- u_int32_t i1, i0;
+ uint32_t i1, i0;
long long int result;
int sign;
result = (((long long int) i0 << 32) | i1) << (j0 - 52);
else
{
- u_int32_t j = i1 + (0x80000000 >> (j0 - 20));
+ uint32_t j = i1 + (0x80000000 >> (j0 - 20));
if (j < i1)
++i0;
__lrint (double x)
{
int32_t j0;
- u_int32_t i0, i1;
+ uint32_t i0, i1;
double w;
double t;
long int result;
__lround (double x)
{
int32_t j0;
- u_int32_t i1, i0;
+ uint32_t i1, i0;
long int result;
int sign;
result = ((long int) i0 << (j0 - 20)) | ((long int) i1 << (j0 - 52));
else
{
- u_int32_t j = i1 + (0x80000000 >> (j0 - 20));
+ uint32_t j = i1 + (0x80000000 >> (j0 - 20));
if (j < i1)
++i0;
__modf (double x, double *iptr)
{
int32_t i0, i1, j0;
- u_int32_t i;
+ uint32_t i;
EXTRACT_WORDS (i0, i1, x);
j0 = ((i0 >> 20) & 0x7ff) - 0x3ff; /* exponent of x */
if (j0 < 20) /* integer part in high x */
}
else /* fraction part in low x */
{
- i = ((u_int32_t) (0xffffffff)) >> (j0 - 20);
+ i = ((uint32_t) (0xffffffff)) >> (j0 - 20);
if ((i1 & i) == 0) /* x is integral */
{
*iptr = x;
__nextup (double x)
{
int32_t hx, ix;
- u_int32_t lx;
+ uint32_t lx;
EXTRACT_WORDS (hx, lx, x);
ix = hx & 0x7fffffff;
__remquo (double x, double y, int *quo)
{
int32_t hx, hy;
- u_int32_t sx, lx, ly;
+ uint32_t sx, lx, ly;
int cquo, qs;
EXTRACT_WORDS (hx, lx, x);
__round (double x)
{
int32_t i0, j0;
- u_int32_t i1;
+ uint32_t i1;
EXTRACT_WORDS (i0, i1, x);
j0 = ((i0 >> 20) & 0x7ff) - 0x3ff;
}
else
{
- u_int32_t i = 0x000fffff >> j0;
+ uint32_t i = 0x000fffff >> j0;
if (((i0 & i) | i1) == 0)
/* X is integral. */
return x;
}
else
{
- u_int32_t i = 0xffffffff >> (j0 - 20);
+ uint32_t i = 0xffffffff >> (j0 - 20);
if ((i1 & i) == 0)
/* X is integral. */
return x;
- u_int32_t j = i1 + (1 << (51 - j0));
+ uint32_t j = i1 + (1 << (51 - j0));
if (j < i1)
i0 += 1;
i1 = j;
__trunc (double x)
{
int32_t i0, j0;
- u_int32_t i1;
+ uint32_t i1;
int sx;
EXTRACT_WORDS (i0, i1, x);
int
__issignaling (double x)
{
- u_int64_t xi;
+ uint64_t xi;
EXTRACT_WORDS64 (xi, x);
#if HIGH_ORDER_BIT_IS_SET_FOR_SNAN
/* We only have to care about the high-order bit of x's significand, because
switch (m) {
case 0: return z ; /* atan(+,+) */
case 1: {
- u_int32_t zh;
+ uint32_t zh;
GET_FLOAT_WORD(zh,z);
SET_FLOAT_WORD(z,zh ^ 0x80000000);
}
return (x*y)/(x*y);
if(hx<hy) return x; /* |x|<|y| return x */
if(hx==hy)
- return Zero[(u_int32_t)sx>>31]; /* |x|=|y| return x*0*/
+ return Zero[(uint32_t)sx>>31]; /* |x|=|y| return x*0*/
/* determine ix = ilogb(x) */
if(hx<0x00800000) { /* subnormal x */
if(hz<0){hx = hx+hx;}
else {
if(hz==0) /* return sign(x)*0 */
- return Zero[(u_int32_t)sx>>31];
+ return Zero[(uint32_t)sx>>31];
hx = hz+hz;
}
}
/* convert back to floating value and restore the sign */
if(hx==0) /* return sign(x)*0 */
- return Zero[(u_int32_t)sx>>31];
+ return Zero[(uint32_t)sx>>31];
while(hx<0x00800000) { /* normalize x */
hx = hx+hx;
iy -= 1;
return 1.0 / x;
}
if (__builtin_expect (hx < 0, 0)
- && (u_int32_t) hx < 0xff800000 && __rintf (x) == x)
+ && (uint32_t) hx < 0xff800000 && __rintf (x) == x)
{
/* Return value for integer x < 0 is NaN with invalid exception. */
*signgamp = 0;
float ret;
{
int32_t i,hx,ix;
- u_int32_t ib;
+ uint32_t ib;
int32_t sign;
float a, b, temp;
}
if (__builtin_expect(hx >= 0x7f800000, 0)) return x+x;
k += (hx>>23)-127;
- i = ((u_int32_t)k&0x80000000)>>31;
+ i = ((uint32_t)k&0x80000000)>>31;
hx = (hx&0x007fffff)|((0x7f-i)<<23);
y = (float)(k+i);
if (FIX_INT_FP_CONVERT_ZERO && y == 0.0f)
}
/* (x<0)**(non-int) is NaN */
- if(__builtin_expect(((((u_int32_t)hx>>31)-1)|yisint)==0, 0))
+ if(__builtin_expect(((((uint32_t)hx>>31)-1)|yisint)==0, 0))
return (x-x)/(x-x);
/* |y| is huge */
}
s = one; /* s (sign of result -ve**odd) = -1 else = 1 */
- if(((((u_int32_t)hx>>31)-1)|(yisint-1))==0)
+ if(((((uint32_t)hx>>31)-1)|(yisint-1))==0)
s = -one; /* (-ve)**(odd int) */
/* compute y * d2 */
}
else if (__builtin_expect((j&0x7fffffff)>0x43160000, 0))/* z <= -150 */
return s*tiny*tiny; /* underflow */
- else if (__builtin_expect((u_int32_t) j==0xc3160000, 0)){/* z == -150*/
+ else if (__builtin_expect((uint32_t) j==0xc3160000, 0)){/* z == -150*/
if(0.0<=(z-d1)) return s*tiny*tiny; /* underflow */
}
/*
if(n<32&&(int32_t)(ix&0xffffff00)!=npio2_hw[n-1]) {
y[0] = r-w; /* quick check no cancellation */
} else {
- u_int32_t high;
+ uint32_t high;
j = ix>>23;
y[0] = r-w;
GET_FLOAT_WORD(high,y[0]);
__ieee754_remainderf(float x, float p)
{
int32_t hx,hp;
- u_int32_t sx;
+ uint32_t sx;
float p_half;
GET_FLOAT_WORD(hx,x);
float z;
int32_t sign = (int)0x80000000;
int32_t ix,s,q,m,t,i;
- u_int32_t r;
+ uint32_t r;
GET_FLOAT_WORD(ix,x);
__ceilf(float x)
{
int32_t i0,j0;
- u_int32_t i;
+ uint32_t i;
GET_FLOAT_WORD(i0,x);
j0 = ((i0>>23)&0xff)-0x7f;
float __copysignf(float x, float y)
{
- u_int32_t ix,iy;
+ uint32_t ix,iy;
GET_FLOAT_WORD(ix,x);
GET_FLOAT_WORD(iy,y);
SET_FLOAT_WORD(x,(ix&0x7fffffff)|(iy&0x80000000));
GET_FLOAT_WORD(hx,x);
ix = hx&0x7fffffff;
if(ix>=0x7f800000) { /* erf(nan)=nan */
- i = ((u_int32_t)hx>>31)<<1;
+ i = ((uint32_t)hx>>31)<<1;
return (float)(1-i)+one/x; /* erf(+-inf)=+-1 */
}
ix = hx&0x7fffffff;
if(ix>=0x7f800000) { /* erfc(nan)=nan */
/* erfc(+-inf)=0,2 */
- float ret = (float)(((u_int32_t)hx>>31)<<1)+one/x;
+ float ret = (float)(((uint32_t)hx>>31)<<1)+one/x;
if (FIX_INT_FP_CONVERT_ZERO && ret == 0.0f)
return 0.0f;
return ret;
{
float y,hi,lo,c,t,e,hxs,hfx,r1;
int32_t k,xsb;
- u_int32_t hx;
+ uint32_t hx;
GET_FLOAT_WORD(hx,x);
xsb = hx&0x80000000; /* sign bit of x */
{
int32_t ix;
GET_FLOAT_WORD(ix,x);
- return (int)((u_int32_t)((ix&0x7f800000)-0x7f800000)>>31);
+ return (int)((uint32_t)((ix&0x7f800000)-0x7f800000)>>31);
}
hidden_def (__finitef)
weak_alias (__finitef, finitef)
__floorf(float x)
{
int32_t i0,j0;
- u_int32_t i;
+ uint32_t i;
GET_FLOAT_WORD(i0,x);
j0 = ((i0>>23)&0xff)-0x7f;
if(j0<23) {
int
__fpclassifyf (float x)
{
- u_int32_t wx;
+ uint32_t wx;
int retval = FP_NORMAL;
GET_FLOAT_WORD (wx, x);
GET_FLOAT_WORD(ix,x);
ix &= 0x7fffffff;
ix = 0x7f800000 - ix;
- return (int)(((u_int32_t)(ix))>>31);
+ return (int)(((uint32_t)(ix))>>31);
}
hidden_def (__isnanf)
weak_alias (__isnanf, isnanf)
int
__issignalingf (float x)
{
- u_int32_t xi;
+ uint32_t xi;
GET_FLOAT_WORD (xi, x);
#if HIGH_ORDER_BIT_IS_SET_FOR_SNAN
/* We only have to care about the high-order bit of x's significand, because
__llrintf (float x)
{
int32_t j0;
- u_int32_t i0;
+ uint32_t i0;
float w;
float t;
long long int result;
__llroundf (float x)
{
int32_t j0;
- u_int32_t i;
+ uint32_t i;
long long int result;
int sign;
__lrintf (float x)
{
int32_t j0;
- u_int32_t i0;
+ uint32_t i0;
float w;
float t;
long int result;
__lroundf (float x)
{
int32_t j0;
- u_int32_t i;
+ uint32_t i;
long int result;
int sign;
__modff(float x, float *iptr)
{
int32_t i0,j0;
- u_int32_t i;
+ uint32_t i;
GET_FLOAT_WORD(i0,x);
j0 = ((i0>>23)&0xff)-0x7f; /* exponent of x */
if(__builtin_expect(j0<23, 1)) { /* integer part in x */
} else {
i = (0x007fffff)>>j0;
if((i0&i)==0) { /* x is integral */
- u_int32_t ix;
+ uint32_t ix;
*iptr = x;
GET_FLOAT_WORD(ix,x);
SET_FLOAT_WORD(x,ix&0x80000000); /* return +-0 */
__remquof (float x, float y, int *quo)
{
int32_t hx,hy;
- u_int32_t sx;
+ uint32_t sx;
int cquo, qs;
GET_FLOAT_WORD (hx, x);
}
else
{
- u_int32_t i = 0x007fffff >> j0;
+ uint32_t i = 0x007fffff >> j0;
if ((i0 & i) == 0)
/* X is integral. */
return x;
__ieee754_acoshl(_Float128 x)
{
_Float128 t;
- u_int64_t lx;
+ uint64_t lx;
int64_t hx;
GET_LDOUBLE_WORDS64(hx,lx,x);
if(hx<0x3fff000000000000LL) { /* x < 1 */
{
_Float128 z;
int64_t k,m,hx,hy,ix,iy;
- u_int64_t lx,ly;
+ uint64_t lx,ly;
GET_LDOUBLE_WORDS64(hx,lx,x);
ix = hx&0x7fffffffffffffffLL;
switch (m) {
case 0: return z ; /* atan(+,+) */
case 1: {
- u_int64_t zh;
+ uint64_t zh;
GET_LDOUBLE_MSW64(zh,z);
SET_LDOUBLE_MSW64(z,zh ^ 0x8000000000000000ULL);
}
__ieee754_atanhl(_Float128 x)
{
_Float128 t;
- u_int32_t jx, ix;
+ uint32_t jx, ix;
ieee854_long_double_shape_type u;
u.value = x;
__ieee754_fmodl (_Float128 x, _Float128 y)
{
int64_t n,hx,hy,hz,ix,iy,sx,i;
- u_int64_t lx,ly,lz;
+ uint64_t lx,ly,lz;
GET_LDOUBLE_WORDS64(hx,lx,x);
GET_LDOUBLE_WORDS64(hy,ly,y);
if(hx<=hy) {
if((hx<hy)||(lx<ly)) return x; /* |x|<|y| return x */
if(lx==ly)
- return Zero[(u_int64_t)sx>>63]; /* |x|=|y| return x*0*/
+ return Zero[(uint64_t)sx>>63]; /* |x|=|y| return x*0*/
}
/* determine ix = ilogb(x) */
if(hz<0){hx = hx+hx+(lx>>63); lx = lx+lx;}
else {
if((hz|lz)==0) /* return sign(x)*0 */
- return Zero[(u_int64_t)sx>>63];
+ return Zero[(uint64_t)sx>>63];
hx = hz+hz+(lz>>63); lx = lz+lz;
}
}
/* convert back to floating value and restore the sign */
if((hx|lx)==0) /* return sign(x)*0 */
- return Zero[(u_int64_t)sx>>63];
+ return Zero[(uint64_t)sx>>63];
while(hx<0x0001000000000000LL) { /* normalize x */
hx = hx+hx+(lx>>63); lx = lx+lx;
iy -= 1;
} else { /* subnormal output */
n = -16382 - iy;
if(n<=48) {
- lx = (lx>>n)|((u_int64_t)hx<<(64-n));
+ lx = (lx>>n)|((uint64_t)hx<<(64-n));
hx >>= n;
} else if (n<=63) {
lx = (hx<<(64-n))|(lx>>n); hx = sx;
__ieee754_gammal_r (_Float128 x, int *signgamp)
{
int64_t hx;
- u_int64_t lx;
+ uint64_t lx;
_Float128 ret;
GET_LDOUBLE_WORDS64 (hx, lx, x);
*signgamp = 0;
return 1.0 / x;
}
- if (hx < 0 && (u_int64_t) hx < 0xffff000000000000ULL && __rintl (x) == x)
+ if (hx < 0 && (uint64_t) hx < 0xffff000000000000ULL && __rintl (x) == x)
{
/* Return value for integer x < 0 is NaN with invalid exception. */
*signgamp = 0;
k=0;
if(ha > 0x5f3f000000000000LL) { /* a>2**8000 */
if(ha >= 0x7fff000000000000LL) { /* Inf or NaN */
- u_int64_t low;
+ uint64_t low;
w = a+b; /* for sNaN */
if (issignaling (a) || issignaling (b))
return w;
}
if(hb < 0x20bf000000000000LL) { /* b < 2**-8000 */
if(hb <= 0x0000ffffffffffffLL) { /* subnormal b or 0 */
- u_int64_t low;
+ uint64_t low;
GET_LDOUBLE_LSW64(low,b);
if((hb|low)==0) return a;
t1=0;
w = __ieee754_sqrtl(t1*y1-(w*(-w)-(t1*y2+t2*b)));
}
if(k!=0) {
- u_int64_t high;
+ uint64_t high;
t1 = 1;
GET_LDOUBLE_MSW64(high,t1);
SET_LDOUBLE_MSW64(t1,high+(k<<48));
_Float128
__ieee754_jnl (int n, _Float128 x)
{
- u_int32_t se;
+ uint32_t se;
int32_t i, ix, sgn;
_Float128 a, b, temp, di, ret;
_Float128 z, w;
_Float128
__ieee754_ynl (int n, _Float128 x)
{
- u_int32_t se;
+ uint32_t se;
int32_t i, ix;
int32_t sign;
_Float128 a, b, temp, ret;
_Float128 y1, t1, t2, r, s, sgn, t, u, v, w;
_Float128 s2, s_h, s_l, t_h, t_l, ay;
int32_t i, j, k, yisint, n;
- u_int32_t ix, iy;
+ uint32_t ix, iy;
int32_t hx, hy;
ieee854_long_double_shape_type o, p, q;
}
/* (x<0)**(non-int) is NaN */
- if (((((u_int32_t) hx >> 31) - 1) | yisint) == 0)
+ if (((((uint32_t) hx >> 31) - 1) | yisint) == 0)
return (x - x) / (x - x);
/* sgn (sign of result -ve**odd) = -1 else = 1 */
sgn = one;
- if (((((u_int32_t) hx >> 31) - 1) | (yisint - 1)) == 0)
+ if (((((uint32_t) hx >> 31) - 1) | (yisint - 1)) == 0)
sgn = -one; /* (-ve)**(odd int) */
/* |y| is huge.
_Float128 z, w, t;
double tx[8];
int64_t exp, n, ix, hx;
- u_int64_t lx;
+ uint64_t lx;
GET_LDOUBLE_WORDS64 (hx, lx, x);
ix = hx & 0x7fffffffffffffffLL;
__ieee754_remainderl(_Float128 x, _Float128 p)
{
int64_t hx,hp;
- u_int64_t sx,lx,lp;
+ uint64_t sx,lx,lp;
_Float128 p_half;
GET_LDOUBLE_WORDS64(hx,lx,x);
__ieee754_sinhl (_Float128 x)
{
_Float128 t, w, h;
- u_int32_t jx, ix;
+ uint32_t jx, ix;
ieee854_long_double_shape_type u;
/* Words of |x|. */
{
_Float128 h, l, z, sin_l, cos_l_m1;
int64_t ix;
- u_int32_t tix, hix, index;
+ uint32_t tix, hix, index;
GET_LDOUBLE_MSW64 (ix, x);
- tix = ((u_int64_t)ix) >> 32;
+ tix = ((uint64_t)ix) >> 32;
tix &= ~0x80000000; /* tix = |x|'s high 32 bits */
if (tix < 0x3ffc3000) /* |x| < 0.1484375 */
{
case 2: index = (hix - 0x3ffc3000) >> 10; break;
}
- SET_LDOUBLE_WORDS64(h, ((u_int64_t)hix) << 32, 0);
+ SET_LDOUBLE_WORDS64(h, ((uint64_t)hix) << 32, 0);
l = y - (h - x);
z = l * l;
sin_l = l*(ONE+z*(SSIN1+z*(SSIN2+z*(SSIN3+z*(SSIN4+z*SSIN5)))));
{
_Float128 h, l, z, sin_l, cos_l_m1;
int64_t ix;
- u_int32_t tix, hix, index;
+ uint32_t tix, hix, index;
GET_LDOUBLE_MSW64 (ix, x);
- tix = ((u_int64_t)ix) >> 32;
+ tix = ((uint64_t)ix) >> 32;
tix &= ~0x80000000; /* tix = |x|'s high 32 bits */
if (tix < 0x3ffc3000) /* |x| < 0.1484375 */
{
case 2: index = (hix - 0x3ffc3000) >> 10; break;
}
- SET_LDOUBLE_WORDS64(h, ((u_int64_t)hix) << 32, 0);
+ SET_LDOUBLE_WORDS64(h, ((uint64_t)hix) << 32, 0);
if (iy)
l = y - (h - x);
else
{
_Float128 h, l, z, sin_l, cos_l_m1;
int64_t ix;
- u_int32_t tix, hix, index;
+ uint32_t tix, hix, index;
GET_LDOUBLE_MSW64 (ix, x);
- tix = ((u_int64_t)ix) >> 32;
+ tix = ((uint64_t)ix) >> 32;
tix &= ~0x80000000; /* tix = |x|'s high 32 bits */
if (tix < 0x3ffc3000) /* |x| < 0.1484375 */
{
case 2: index = (hix - 0x3ffc3000) >> 10; break;
}
- SET_LDOUBLE_WORDS64(h, ((u_int64_t)hix) << 32, 0);
+ SET_LDOUBLE_WORDS64(h, ((uint64_t)hix) << 32, 0);
if (iy)
l = (ix < 0 ? -y : y) - (h - x);
else
_Float128 __ceill(_Float128 x)
{
int64_t i0,i1,j0;
- u_int64_t i,j;
+ uint64_t i,j;
GET_LDOUBLE_WORDS64(i0,i1,x);
j0 = ((i0>>48)&0x7fff)-0x3fff;
if(j0<48) {
_Float128 __copysignl(_Float128 x, _Float128 y)
{
- u_int64_t hx,hy;
+ uint64_t hx,hy;
GET_LDOUBLE_MSW64(hx,x);
GET_LDOUBLE_MSW64(hy,y);
SET_LDOUBLE_MSW64(x,(hx&0x7fffffffffffffffULL)
if (ix >= 0x7fff0000)
{ /* erfc(nan)=nan */
/* erfc(+-inf)=0,2 */
- return (_Float128) (((u_int32_t) sign >> 31) << 1) + one / x;
+ return (_Float128) (((uint32_t) sign >> 31) << 1) + one / x;
}
if (ix < 0x3ffd0000) /* |x| <1/4 */
_Float128 __fabsl(_Float128 x)
{
- u_int64_t hx;
+ uint64_t hx;
GET_LDOUBLE_MSW64(hx,x);
SET_LDOUBLE_MSW64(x,hx&0x7fffffffffffffffLL);
return x;
{
int64_t hx;
GET_LDOUBLE_MSW64(hx,x);
- return (int)((u_int64_t)((hx&0x7fff000000000000LL)
- -0x7fff000000000000LL)>>63);
+ return (int)((uint64_t)((hx&0x7fff000000000000LL)
+ -0x7fff000000000000LL)>>63);
}
mathx_hidden_def (__finitel)
weak_alias (__finitel, finitel)
_Float128 __floorl(_Float128 x)
{
int64_t i0,i1,j0;
- u_int64_t i,j;
+ uint64_t i,j;
GET_LDOUBLE_WORDS64(i0,i1,x);
j0 = ((i0>>48)&0x7fff)-0x3fff;
if(j0<48) {
int
__fpclassifyl (_Float128 x)
{
- u_int64_t hx, lx;
+ uint64_t hx, lx;
int retval = FP_NORMAL;
GET_LDOUBLE_WORDS64 (hx, lx, x);
_Float128 __frexpl(_Float128 x, int *eptr)
{
- u_int64_t hx, lx, ix;
+ uint64_t hx, lx, ix;
GET_LDOUBLE_WORDS64(hx,lx,x);
ix = 0x7fffffffffffffffULL&hx;
*eptr = 0;
int64_t hx,lx;
GET_LDOUBLE_WORDS64(hx,lx,x);
hx &= 0x7fffffffffffffffLL;
- hx |= (u_int64_t)(lx|(-lx))>>63;
+ hx |= (uint64_t)(lx|(-lx))>>63;
hx = 0x7fff000000000000LL - hx;
- return (int)((u_int64_t)hx>>63);
+ return (int)((uint64_t)hx>>63);
}
mathx_hidden_def (__isnanl)
weak_alias (__isnanl, isnanl)
int
__issignalingl (_Float128 x)
{
- u_int64_t hxi, lxi __attribute__ ((unused));
+ uint64_t hxi, lxi __attribute__ ((unused));
GET_LDOUBLE_WORDS64 (hxi, lxi, x);
#if HIGH_ORDER_BIT_IS_SET_FOR_SNAN
/* We only have to care about the high-order bit of x's significand, because
__llrintl (_Float128 x)
{
int32_t j0;
- u_int64_t i0,i1;
+ uint64_t i0,i1;
_Float128 w;
_Float128 t;
long long int result;
__llroundl (_Float128 x)
{
int64_t j0;
- u_int64_t i1, i0;
+ uint64_t i1, i0;
long long int result;
int sign;
result = ((long long int) i0 << (j0 - 48)) | (i1 << (j0 - 112));
else
{
- u_int64_t j = i1 + (0x8000000000000000ULL >> (j0 - 48));
+ uint64_t j = i1 + (0x8000000000000000ULL >> (j0 - 48));
if (j < i1)
++i0;
__lrintl (_Float128 x)
{
int32_t j0;
- u_int64_t i0,i1;
+ uint64_t i0,i1;
_Float128 w;
_Float128 t;
long int result;
__lroundl (_Float128 x)
{
int64_t j0;
- u_int64_t i1, i0;
+ uint64_t i1, i0;
long int result;
int sign;
result = ((long int) i0 << (j0 - 48)) | (i1 << (j0 - 112));
else
{
- u_int64_t j = i1 + (0x8000000000000000ULL >> (j0 - 48));
+ uint64_t j = i1 + (0x8000000000000000ULL >> (j0 - 48));
if (j < i1)
++i0;
_Float128 __modfl(_Float128 x, _Float128 *iptr)
{
int64_t i0,i1,j0;
- u_int64_t i;
+ uint64_t i;
GET_LDOUBLE_WORDS64(i0,i1,x);
j0 = ((i0>>48)&0x7fff)-0x3fff; /* exponent of x */
if(j0<48) { /* integer part in high x */
{
fenv_t env;
int64_t i0,j0,sx;
- u_int64_t i1 __attribute__ ((unused));
+ uint64_t i1 __attribute__ ((unused));
_Float128 w,t;
GET_LDOUBLE_WORDS64(i0,i1,x);
- sx = (((u_int64_t)i0)>>63);
+ sx = (((uint64_t)i0)>>63);
j0 = ((i0>>48)&0x7fff)-0x3fff;
if(j0<112) {
if(j0<0) {
_Float128 __nextafterl(_Float128 x, _Float128 y)
{
int64_t hx,hy,ix,iy;
- u_int64_t lx,ly;
+ uint64_t lx,ly;
GET_LDOUBLE_WORDS64(hx,lx,x);
GET_LDOUBLE_WORDS64(hy,ly,y);
{
int32_t hx,ix;
int64_t hy,iy;
- u_int32_t lx;
- u_int64_t ly;
+ uint32_t lx;
+ uint64_t ly;
EXTRACT_WORDS(hx,lx,x);
GET_LDOUBLE_WORDS64(hy,ly,y);
if((long double) x==y) return y; /* x=y, return y */
if((ix|lx)==0) { /* x == 0 */
double u;
- INSERT_WORDS(x,(u_int32_t)((hy>>32)&0x80000000),1);/* return +-minsub */
+ INSERT_WORDS(x,(uint32_t)((hy>>32)&0x80000000),1);/* return +-minsub */
u = math_opt_barrier (x);
u = u * u;
math_force_eval (u); /* raise underflow flag */
{
int32_t hx,ix;
int64_t hy,iy;
- u_int64_t ly;
+ uint64_t ly;
GET_FLOAT_WORD(hx,x);
GET_LDOUBLE_WORDS64(hy,ly,y);
if((long double) x==y) return y; /* x=y, return y */
if(ix==0) { /* x == 0 */
float u;
- SET_FLOAT_WORD(x,(u_int32_t)((hy>>32)&0x80000000)|1);/* return +-minsub*/
+ SET_FLOAT_WORD(x,(uint32_t)((hy>>32)&0x80000000)|1);/* return +-minsub*/
u = math_opt_barrier (x);
u = u * u;
math_force_eval (u); /* raise underflow flag */
__nextupl (_Float128 x)
{
int64_t hx, ix;
- u_int64_t lx;
+ uint64_t lx;
GET_LDOUBLE_WORDS64 (hx, lx, x);
ix = hx & 0x7fffffffffffffffLL;
__remquol (_Float128 x, _Float128 y, int *quo)
{
int64_t hx,hy;
- u_int64_t sx,lx,ly,qs;
+ uint64_t sx,lx,ly,qs;
int cquo;
GET_LDOUBLE_WORDS64 (hx, lx, x);
_Float128 __rintl(_Float128 x)
{
int64_t i0,j0,sx;
- u_int64_t i1 __attribute__ ((unused));
+ uint64_t i1 __attribute__ ((unused));
_Float128 w,t;
GET_LDOUBLE_WORDS64(i0,i1,x);
- sx = (((u_int64_t)i0)>>63);
+ sx = (((uint64_t)i0)>>63);
j0 = ((i0>>48)&0x7fff)-0x3fff;
if(j0<112) {
if(j0<0) {
__roundl (_Float128 x)
{
int32_t j0;
- u_int64_t i1, i0;
+ uint64_t i1, i0;
GET_LDOUBLE_WORDS64 (i0, i1, x);
j0 = ((i0 >> 48) & 0x7fff) - 0x3fff;
}
else
{
- u_int64_t i = 0x0000ffffffffffffLL >> j0;
+ uint64_t i = 0x0000ffffffffffffLL >> j0;
if (((i0 & i) | i1) == 0)
/* X is integral. */
return x;
}
else
{
- u_int64_t i = -1ULL >> (j0 - 48);
+ uint64_t i = -1ULL >> (j0 - 48);
if ((i1 & i) == 0)
/* X is integral. */
return x;
- u_int64_t j = i1 + (1LL << (111 - j0));
+ uint64_t j = i1 + (1LL << (111 - j0));
if (j < i1)
i0 += 1;
i1 = j;
__tanhl (_Float128 x)
{
_Float128 t, z;
- u_int32_t jx, ix;
+ uint32_t jx, ix;
ieee854_long_double_shape_type u;
/* Words of |x|. */
__truncl (_Float128 x)
{
int32_t j0;
- u_int64_t i0, i1, sx;
+ uint64_t i0, i1, sx;
GET_LDOUBLE_WORDS64 (i0, i1, x);
sx = i0 & 0x8000000000000000ULL;
if(hz<0){hx = hx+hx+(lx>>63); lx = lx+lx;}
else {
if((hz|lz)==0) /* return sign(x)*0 */
- return Zero[(u_int64_t)sx>>63];
+ return Zero[(uint64_t)sx>>63];
hx = hz+hz+(lz>>63); lx = lz+lz;
}
}
/* convert back to floating value and restore the sign */
if((hx|lx)==0) /* return sign(x)*0 */
- return Zero[(u_int64_t)sx>>63];
+ return Zero[(uint64_t)sx>>63];
while(hx<0x0001000000000000LL) { /* normalize x */
hx = hx+hx+(lx>>63); lx = lx+lx;
iy -= 1;
n = -1022 - iy;
/* We know 1 <= N <= 52, and that there are no nonzero
bits in places below 2^-1074. */
- lx = (lx >> n) | ((u_int64_t) hx << (64 - n));
+ lx = (lx >> n) | ((uint64_t) hx << (64 - n));
hx >>= n;
x = ldbl_insert_mantissa((sx>>63), -1023, hx, lx);
x *= one; /* create necessary signal */
*signgamp = 0;
return 1.0 / x;
}
- if (hx < 0 && (u_int64_t) hx < 0xfff0000000000000ULL && __rintl (x) == x)
+ if (hx < 0 && (uint64_t) hx < 0xfff0000000000000ULL && __rintl (x) == x)
{
/* Return value for integer x < 0 is NaN with invalid exception. */
*signgamp = 0;
}
/* (x<0)**(non-int) is NaN */
- if (((((u_int32_t) hx >> 31) - 1) | yisint) == 0)
+ if (((((uint32_t) hx >> 31) - 1) | yisint) == 0)
return (x - x) / (x - x);
/* sgn (sign of result -ve**odd) = -1 else = 1 */
sgn = one;
- if (((((u_int32_t) hx >> 31) - 1) | (yisint - 1)) == 0)
+ if (((((uint32_t) hx >> 31) - 1) | (yisint - 1)) == 0)
sgn = -one; /* (-ve)**(odd int) */
/* |y| is huge.
double tx[8];
int exp;
int64_t n, ix, hx, ixd;
- u_int64_t lxd;
+ uint64_t lxd;
double xhi;
xhi = ldbl_high (x);
__ieee754_remainderl(long double x, long double p)
{
int64_t hx,hp;
- u_int64_t sx,lx,lp;
+ uint64_t sx,lx,lp;
long double p_half;
double xhi, xlo, phi, plo;
xhi = ldbl_high (x);
EXTRACT_WORDS64 (ix, xhi);
- tix = ((u_int64_t)ix) >> 32;
+ tix = ((uint64_t)ix) >> 32;
tix &= ~0x80000000; /* tix = |x|'s high 32 bits */
if (tix < 0x3fc30000) /* |x| < 0.1484375 */
{
{
long double h, l, z, sin_l, cos_l_m1;
int64_t ix;
- u_int32_t tix, hix, index;
+ uint32_t tix, hix, index;
double xhi, hhi;
xhi = ldbl_high (x);
EXTRACT_WORDS64 (ix, xhi);
- tix = ((u_int64_t)ix) >> 32;
+ tix = ((uint64_t)ix) >> 32;
tix &= ~0x80000000; /* tix = |x|'s high 32 bits */
if (tix < 0x3fc30000) /* |x| < 0.1484375 */
{
long double __fabsl(long double x)
{
- u_int64_t hx, lx;
+ uint64_t hx, lx;
double xhi, xlo;
ldbl_unpack (x, &xhi, &xlo);
int
___fpclassifyl (long double x)
{
- u_int64_t hx, lx;
+ uint64_t hx, lx;
int retval = FP_NORMAL;
double xhi, xlo;
long double __modfl(long double x, long double *iptr)
{
int64_t i0,i1,j0;
- u_int64_t i;
+ uint64_t i;
double xhi, xlo;
ldbl_unpack (x, &xhi, &xlo);
if((long double) x==y) return y; /* x=y, return y */
if(ix==0) { /* x == 0 */
float u;
- SET_FLOAT_WORD(x,(u_int32_t)((hy>>32)&0x80000000)|1);/* return +-minsub*/
+ SET_FLOAT_WORD(x,(uint32_t)((hy>>32)&0x80000000)|1);/* return +-minsub*/
u = math_opt_barrier (x);
u = u * u;
math_force_eval (u); /* raise underflow flag */
__remquol (long double x, long double y, int *quo)
{
int64_t hx,hy;
- u_int64_t sx,lx,ly,qs;
+ uint64_t sx,lx,ly,qs;
int cquo;
double xhi, xlo, yhi, ylo;
__ieee754_acoshl(long double x)
{
long double t;
- u_int32_t se,i0,i1;
+ uint32_t se,i0,i1;
GET_LDOUBLE_WORDS(se,i0,i1,x);
if(se<0x3fff || se & 0x8000) { /* x < 1 */
return (x-x)/(x-x);
{
long double t, w, p, q, c, r, s;
int32_t ix;
- u_int32_t se, i0, i1, k;
+ uint32_t se, i0, i1, k;
GET_LDOUBLE_WORDS (se, i0, i1, x);
ix = se & 0x7fff;
{
long double t;
int32_t ix;
- u_int32_t se,i0,i1;
+ uint32_t se,i0,i1;
GET_LDOUBLE_WORDS(se,i0,i1,x);
ix = se&0x7fff;
if ((ix+((((i0&0x7fffffff)|i1)|(-((i0&0x7fffffff)|i1)))>>31))>0x3fff)
{
long double t,w;
int32_t ex;
- u_int32_t mx,lx;
+ uint32_t mx,lx;
/* High word of |x|. */
GET_LDOUBLE_WORDS(ex,mx,lx,x);
long double
__ieee754_gammal_r (long double x, int *signgamp)
{
- u_int32_t es, hx, lx;
+ uint32_t es, hx, lx;
long double ret;
GET_LDOUBLE_WORDS (es, hx, lx, x);
long double __ieee754_hypotl(long double x, long double y)
{
long double a,b,t1,t2,y1,y2,w;
- u_int32_t j,k,ea,eb;
+ uint32_t j,k,ea,eb;
GET_LDOUBLE_EXP(ea,x);
ea &= 0x7fff;
k=0;
if(__builtin_expect(ea > 0x5f3f,0)) { /* a>2**8000 */
if(ea == 0x7fff) { /* Inf or NaN */
- u_int32_t exp __attribute__ ((unused));
- u_int32_t high,low;
+ uint32_t exp __attribute__ ((unused));
+ uint32_t high,low;
w = a+b; /* for sNaN */
if (issignaling (a) || issignaling (b))
return w;
}
if(__builtin_expect(eb < 0x20bf, 0)) { /* b < 2**-8000 */
if(eb == 0) { /* subnormal b or 0 */
- u_int32_t exp __attribute__ ((unused));
- u_int32_t high,low;
+ uint32_t exp __attribute__ ((unused));
+ uint32_t high,low;
GET_LDOUBLE_WORDS(exp,high,low,b);
if((high|low)==0) return a;
SET_LDOUBLE_WORDS(t1, 0x7ffd, 0x80000000, 0); /* t1=2^16382 */
/* medium size a and b */
w = a-b;
if (w>b) {
- u_int32_t high;
+ uint32_t high;
GET_LDOUBLE_MSW(high,a);
SET_LDOUBLE_WORDS(t1,ea,high,0);
t2 = a-t1;
w = __ieee754_sqrtl(t1*t1-(b*(-b)-t2*(a+t1)));
} else {
- u_int32_t high;
+ uint32_t high;
GET_LDOUBLE_MSW(high,b);
a = a+a;
SET_LDOUBLE_WORDS(y1,eb,high,0);
w = __ieee754_sqrtl(t1*y1-(w*(-w)-(t1*y2+t2*b)));
}
if(k!=0) {
- u_int32_t exp;
+ uint32_t exp;
t1 = 1.0;
GET_LDOUBLE_EXP(exp,t1);
SET_LDOUBLE_EXP(t1,exp+k);
{
long double z, s, c, ss, cc, r, u, v;
int32_t ix;
- u_int32_t se;
+ uint32_t se;
GET_LDOUBLE_EXP (se, x);
ix = se & 0x7fff;
{
long double z, s, c, ss, cc, u, v;
int32_t ix;
- u_int32_t se, i0, i1;
+ uint32_t se, i0, i1;
GET_LDOUBLE_WORDS (se, i0, i1, x);
ix = se & 0x7fff;
const long double *p, *q;
long double z, r, s;
int32_t ix;
- u_int32_t se, i0, i1;
+ uint32_t se, i0, i1;
GET_LDOUBLE_WORDS (se, i0, i1, x);
ix = se & 0x7fff;
const long double *p, *q;
long double s, r, z;
int32_t ix;
- u_int32_t se, i0, i1;
+ uint32_t se, i0, i1;
GET_LDOUBLE_WORDS (se, i0, i1, x);
ix = se & 0x7fff;
{
long double z, c, r, s, ss, cc, u, v, y;
int32_t ix;
- u_int32_t se;
+ uint32_t se;
GET_LDOUBLE_EXP (se, x);
ix = se & 0x7fff;
{
long double z, s, c, ss, cc, u, v;
int32_t ix;
- u_int32_t se, i0, i1;
+ uint32_t se, i0, i1;
GET_LDOUBLE_WORDS (se, i0, i1, x);
ix = se & 0x7fff;
const long double *p, *q;
long double z, r, s;
int32_t ix;
- u_int32_t se, i0, i1;
+ uint32_t se, i0, i1;
GET_LDOUBLE_WORDS (se, i0, i1, x);
ix = se & 0x7fff;
const long double *p, *q;
static long double s, r, z;
int32_t ix;
- u_int32_t se, i0, i1;
+ uint32_t se, i0, i1;
GET_LDOUBLE_WORDS (se, i0, i1, x);
ix = se & 0x7fff;
long double
__ieee754_jnl (int n, long double x)
{
- u_int32_t se, i0, i1;
+ uint32_t se, i0, i1;
int32_t i, ix, sgn;
long double a, b, temp, di, ret;
long double z, w;
long double
__ieee754_ynl (int n, long double x)
{
- u_int32_t se, i0, i1;
+ uint32_t se, i0, i1;
int32_t i, ix;
int32_t sign;
long double a, b, temp, ret;
{
long double y, z;
int n, ix;
- u_int32_t se, i0, i1;
+ uint32_t se, i0, i1;
GET_LDOUBLE_WORDS (se, i0, i1, x);
ix = se & 0x7fff;
{
long double t, y, z, nadj, p, p1, p2, q, r, w;
int i, ix;
- u_int32_t se, i0, i1;
+ uint32_t se, i0, i1;
*signgamp = 1;
GET_LDOUBLE_WORDS (se, i0, i1, x);
{
double tx[3], ty[3];
int32_t se, j0;
- u_int32_t i0, i1;
+ uint32_t i0, i1;
int sx;
int n, exp;
__ieee754_sinhl(long double x)
{
long double t,w,h;
- u_int32_t jx,ix,i0,i1;
+ uint32_t jx,ix,i0,i1;
/* Words of |x|. */
GET_LDOUBLE_WORDS(jx,i0,i1,x);
long double __copysignl(long double x, long double y)
{
- u_int32_t es1,es2;
+ uint32_t es1,es2;
GET_LDOUBLE_EXP(es1,x);
GET_LDOUBLE_EXP(es2,y);
SET_LDOUBLE_EXP(x,(es1&0x7fff)|(es2&0x8000));
{
long double R, S, P, Q, s, y, z, r;
int32_t ix, i;
- u_int32_t se, i0, i1;
+ uint32_t se, i0, i1;
GET_LDOUBLE_WORDS (se, i0, i1, x);
ix = se & 0x7fff;
{
int32_t hx, ix;
long double R, S, P, Q, s, y, z, r;
- u_int32_t se, i0, i1;
+ uint32_t se, i0, i1;
GET_LDOUBLE_WORDS (se, i0, i1, x);
ix = se & 0x7fff;
long double __frexpl(long double x, int *eptr)
{
- u_int32_t se, hx, ix, lx;
+ uint32_t se, hx, ix, lx;
GET_LDOUBLE_WORDS(se,hx,lx,x);
ix = 0x7fff&se;
*eptr = 0;
int
__issignalingl (long double x)
{
- u_int32_t exi, hxi, lxi;
+ uint32_t exi, hxi, lxi;
GET_LDOUBLE_WORDS (exi, hxi, lxi, x);
#if HIGH_ORDER_BIT_IS_SET_FOR_SNAN
# error not implemented
__llrintl (long double x)
{
int32_t se,j0;
- u_int32_t i0, i1;
+ uint32_t i0, i1;
long long int result;
long double w;
long double t;
__llroundl (long double x)
{
int32_t j0;
- u_int32_t se, i1, i0;
+ uint32_t se, i1, i0;
long long int result;
int sign;
return j0 < -1 ? 0 : sign;
else
{
- u_int32_t j = i0 + (0x40000000 >> j0);
+ uint32_t j = i0 + (0x40000000 >> j0);
if (j < i0)
{
j >>= 1;
result = (((long long int) i0 << 32) | i1) << (j0 - 63);
else
{
- u_int32_t j = i1 + (0x80000000 >> (j0 - 31));
+ uint32_t j = i1 + (0x80000000 >> (j0 - 31));
result = (long long int) i0;
if (j < i1)
__lrintl (long double x)
{
int32_t se,j0;
- u_int32_t i0, i1;
+ uint32_t i0, i1;
long int result;
long double w;
long double t;
__lroundl (long double x)
{
int32_t j0;
- u_int32_t se, i1, i0;
+ uint32_t se, i1, i0;
long int result;
int sign;
return j0 < -1 ? 0 : sign;
else
{
- u_int32_t j = i0 + (0x40000000 >> j0);
+ uint32_t j = i0 + (0x40000000 >> j0);
if (j < i0)
{
j >>= 1;
result = ((long int) i0 << (j0 - 31)) | (i1 << (j0 - 63));
else
{
- u_int32_t j = i1 + (0x80000000 >> (j0 - 31));
+ uint32_t j = i1 + (0x80000000 >> (j0 - 31));
unsigned long int ures = i0;
if (j < i1)
__modfl(long double x, long double *iptr)
{
int32_t i0,i1,j0;
- u_int32_t i,se;
+ uint32_t i,se;
GET_LDOUBLE_WORDS(se,i0,i1,x);
j0 = (se&0x7fff)-0x3fff; /* exponent of x */
if(j0<32) { /* integer part in high x */
SET_LDOUBLE_WORDS(x,se&0x8000,0,0); /* return +-0 */
return x;
} else { /* fraction part in low x */
- i = ((u_int32_t)(0x7fffffff))>>(j0-32);
+ i = ((uint32_t)(0x7fffffff))>>(j0-32);
if((i1&i)==0) { /* x is integral */
*iptr = x;
SET_LDOUBLE_WORDS(x,se&0x8000,0,0); /* return +-0 */
double __nexttoward(double x, long double y)
{
int32_t hx,ix,iy;
- u_int32_t lx,hy,ly,esy;
+ uint32_t lx,hy,ly,esy;
EXTRACT_WORDS(hx,lx,x);
GET_LDOUBLE_WORDS(esy,hy,ly,y);
float __nexttowardf(float x, long double y)
{
int32_t hx,ix,iy;
- u_int32_t hy,ly,esy;
+ uint32_t hy,ly,esy;
GET_FLOAT_WORD(hx,x);
GET_LDOUBLE_WORDS(esy,hy,ly,y);
long double
__nextupl (long double x)
{
- u_int32_t hx, ix;
- u_int32_t lx;
+ uint32_t hx, ix;
+ uint32_t lx;
int32_t esx;
GET_LDOUBLE_WORDS (esx, hx, lx, x);
__remquol (long double x, long double p, int *quo)
{
int32_t ex,ep,hx,hp;
- u_int32_t sx,lx,lp;
+ uint32_t sx,lx,lp;
int cquo,qs;
GET_LDOUBLE_WORDS (ex, hx, lx, x);
__roundl (long double x)
{
int32_t j0;
- u_int32_t se, i1, i0;
+ uint32_t se, i1, i0;
GET_LDOUBLE_WORDS (se, i0, i1, x);
j0 = (se & 0x7fff) - 0x3fff;
}
else
{
- u_int32_t i = 0x7fffffff >> j0;
+ uint32_t i = 0x7fffffff >> j0;
if (((i0 & i) | i1) == 0)
/* X is integral. */
return x;
- u_int32_t j = i0 + (0x40000000 >> j0);
+ uint32_t j = i0 + (0x40000000 >> j0);
if (j < i0)
se += 1;
i0 = (j & ~i) | 0x80000000;
}
else
{
- u_int32_t i = 0xffffffff >> (j0 - 31);
+ uint32_t i = 0xffffffff >> (j0 - 31);
if ((i1 & i) == 0)
/* X is integral. */
return x;
- u_int32_t j = i1 + (1 << (62 - j0));
+ uint32_t j = i1 + (1 << (62 - j0));
if (j < i1)
{
- u_int32_t k = i0 + 1;
+ uint32_t k = i0 + 1;
if (k < i0)
{
se += 1;
{
long double t,z;
int32_t se;
- u_int32_t j0,j1,ix;
+ uint32_t j0,j1,ix;
/* High word of |x|. */
GET_LDOUBLE_WORDS(se,j0,j1,x);
float __nldbl_nexttowardf(float x, double y)
{
int32_t hx,hy,ix,iy;
- u_int32_t ly;
+ uint32_t ly;
GET_FLOAT_WORD(hx,x);
EXTRACT_WORDS(hy,ly,y);
if((double) x==y) return y; /* x=y, return y */
if(ix==0) { /* x == 0 */
float u;
- SET_FLOAT_WORD(x,(u_int32_t)(hy&0x80000000)|1);/* return +-minsub*/
+ SET_FLOAT_WORD(x,(uint32_t)(hy&0x80000000)|1);/* return +-minsub*/
u = math_opt_barrier (x);
u = u * u;
math_force_eval (u); /* raise underflow flag */
this format and rounding won't change the result. */
{
int32_t exponent;
- u_int32_t i0, i1;
+ uint32_t i0, i1;
GET_LDOUBLE_WORDS (exponent, i0, i1, y);
exponent = (exponent & 0x7fff) - 0x3fff;
if (exponent <= 31
int
__fpclassifyl (long double x)
{
- u_int32_t ex, hx, lx;
+ uint32_t ex, hx, lx;
int retval = FP_NORMAL;
GET_LDOUBLE_WORDS (ex, hx, lx, x);
__llrint (double x)
{
int32_t e;
- u_int32_t h, l, s;
+ uint32_t h, l, s;
long long int result;
x = __m81_u(__rint) (x);
__llrintf (float x)
{
int32_t e;
- u_int32_t i, s;
+ uint32_t i, s;
long long int result;
x = __m81_u(__rintf) (x);
__llrintl (long double x)
{
int32_t e, s;
- u_int32_t h, l;
+ uint32_t h, l;
long long int result;
x = __m81_u(__rintl) (x);
long double __nextafterl(long double x, long double y)
{
int32_t ix,iy,esx,esy;
- u_int32_t hx,hy,lx,ly;
+ uint32_t hx,hy,lx,ly;
GET_LDOUBLE_WORDS(esx,hx,lx,x);
GET_LDOUBLE_WORDS(esy,hy,ly,y);
corrected for by adding log2 (e) * X_FRAC_LOW to the final
result. */
int32_t se;
- u_int32_t i0, i1;
+ uint32_t i0, i1;
GET_LDOUBLE_WORDS (se, i0, i1, x_frac);
x_frac_low = x_frac;
i1 &= 0xffffffe0;
long double w = x_frac - 1;
long double w_hi, w_lo;
int32_t se;
- u_int32_t i0, i1;
+ uint32_t i0, i1;
GET_LDOUBLE_WORDS (se, i0, i1, w);
i0 &= 0xffff0000;
i1 = 0;
projects / platform / upstream / glibc.git / commitdiff