/sysdeps/ieee754/ldbl-128ibm/ |
A D | s_nextafterl.c | 57 u = u * u; in __nextafterl() 62 long double u; in __nextafterl() local 78 return u; in __nextafterl() 86 u = u * u; in __nextafterl() 107 u = yhi * 0x1p-105; in __nextafterl() 110 u = yhi; in __nextafterl() 112 return x - u; in __nextafterl() 122 return u; in __nextafterl() 130 u = u * u; in __nextafterl() 153 u = yhi; in __nextafterl() [all …]
|
A D | s_atanl.c | 179 long double t, u, p, q; in __atanl() local 232 u = 0.125 * k; in __atanl() 234 t = (x - u) / (1.0 + x * u); in __atanl() 238 u = t * t; in __atanl() 239 p = ((((p4 * u) + p3) * u + p2) * u + p1) * u + p0; in __atanl() 240 q = ((((u + q4) * u + q3) * u + q2) * u + q1) * u + q0; in __atanl() 241 u = t * u * p / q + t; in __atanl() 244 u = atantbl[k] + u; in __atanl() 246 return (-u); in __atanl() 248 return u; in __atanl()
|
A D | math_ldbl.h | 39 u.ld = x; in ldbl_extract_mantissa() 42 lo = ((uint64_t) u.d[1].ieee.mantissa0 << 32) | u.d[1].ieee.mantissa1; in ldbl_extract_mantissa() 73 ediff = u.d[0].ieee.exponent - u.d[1].ieee.exponent - 53; in ldbl_extract_mantissa() 84 if (u.d[0].ieee.negative != u.d[1].ieee.negative in ldbl_extract_mantissa() 181 return u.ld; in ldbl_insert_mantissa() 190 u.d[0].d = a; in default_ldbl_pack() 191 u.d[1].d = aa; in default_ldbl_pack() 192 return u.ld; in default_ldbl_pack() 199 u.ld = l; in default_ldbl_unpack() 200 *a = u.d[0].d; in default_ldbl_unpack() [all …]
|
A D | math-nan-payload-ldouble.h | 22 union ibm_extended_long_double u; \ 23 u.ld = (flt); \ 24 u.d[0].ieee_nan.mantissa0 = (mant) >> 32; \ 25 u.d[0].ieee_nan.mantissa1 = (mant); \ 26 if ((u.d[0].ieee.mantissa0 | u.d[0].ieee.mantissa1) != 0) \ 27 (flt) = u.ld; \
|
A D | mpn2ldbl.c | 33 union ibm_extended_long_double u; in __mpn_construct_long_double() local 38 u.d[0].ieee.negative = sign; in __mpn_construct_long_double() 39 u.d[1].ieee.negative = sign; in __mpn_construct_long_double() 41 u.d[1].ieee.exponent = 0; in __mpn_construct_long_double() 81 u.d[0].ieee.exponent = 0; in __mpn_construct_long_double() 115 u.d[0].ieee.exponent++; in __mpn_construct_long_double() 125 u.d[1].ieee.negative = !sign; in __mpn_construct_long_double() 150 u.d[1].ieee.negative = 0; in __mpn_construct_long_double() 152 u.d[1].ieee.mantissa1 = lo; in __mpn_construct_long_double() 154 u.d[0].ieee.mantissa1 = hi; in __mpn_construct_long_double() [all …]
|
A D | s_nexttowardf.c | 46 float u; in __nexttowardf() local 48 u = math_opt_barrier (x); in __nexttowardf() 49 u = u * u; in __nexttowardf() 50 math_force_eval (u); /* raise underflow flag */ in __nexttowardf() 68 float u = x+x; /* overflow */ in __nexttowardf() local 69 math_force_eval (u); in __nexttowardf() 73 float u = x*x; in __nexttowardf() local 74 math_force_eval (u); /* raise underflow flag */ in __nexttowardf()
|
/sysdeps/ieee754/ldbl-96/ |
A D | s_cbrtl.c | 42 long double xm, u; in __cbrtl() local 55 u = (((-1.34661104733595206551E-1 * xm in __cbrtl() 61 u *= factor[2 + xe % 3]; in __cbrtl() 62 u = __ldexpl (x > 0.0 ? u : -u, xe / 3); in __cbrtl() 64 u -= (u - (x / (u * u))) * third; in __cbrtl() 65 u -= (u - (x / (u * u))) * third; in __cbrtl() 66 return u; in __cbrtl()
|
A D | math-nan-payload-ldouble.h | 22 union ieee854_long_double u; \ 23 u.d = (flt); \ 24 u.ieee_nan.mantissa0 = (mant) >> 32; \ 25 u.ieee_nan.mantissa1 = (mant); \ 26 if ((u.ieee.mantissa0 | u.ieee.mantissa1) != 0) \ 27 (flt) = u.d; \
|
A D | mpn2ldbl.c | 31 union ieee854_long_double u; in __mpn_construct_long_double() local 33 u.ieee.negative = sign; in __mpn_construct_long_double() 34 u.ieee.exponent = expt + IEEE854_LONG_DOUBLE_BIAS; in __mpn_construct_long_double() 36 u.ieee.mantissa1 = frac_ptr[0]; in __mpn_construct_long_double() 37 u.ieee.mantissa0 = frac_ptr[1]; in __mpn_construct_long_double() 39 u.ieee.mantissa1 = frac_ptr[0] & (((mp_limb_t) 1 << 32) - 1); in __mpn_construct_long_double() 40 u.ieee.mantissa0 = frac_ptr[0] >> 32; in __mpn_construct_long_double() 45 return u.d; in __mpn_construct_long_double()
|
A D | s_fmal.c | 37 u.d = x; in __fmal() 151 if (u.ieee.exponent) in __fmal() 154 u.d *= 0x1p64L; in __fmal() 174 x = u.d; in __fmal() 224 u.d = a2 + m2; in __fmal() 228 if ((u.ieee.mantissa1 & 1) == 0 && u.ieee.exponent != 0x7fff) in __fmal() 232 return a1 + u.d; in __fmal() 236 if ((u.ieee.mantissa1 & 1) == 0 && u.ieee.exponent != 0x7fff) in __fmal() 246 v.d = a1 + u.d; in __fmal() 253 asm volatile ("" : "=m" (u) : "m" (u)); in __fmal() [all …]
|
A D | s_nexttowardf.c | 40 float u; in __nexttowardf() local 42 u = math_opt_barrier (x); in __nexttowardf() 43 u = u * u; in __nexttowardf() 44 math_force_eval (u); /* raise underflow flag */ in __nexttowardf() 62 float u = x+x; /* overflow */ in __nexttowardf() local 63 math_force_eval (u); in __nexttowardf() 67 float u = x*x; /* underflow */ in __nexttowardf() local 68 math_force_eval (u); /* raise underflow flag */ in __nexttowardf()
|
/sysdeps/ieee754/float128/ |
A D | math-nan-payload-float128.h | 24 union ieee854_float128 u; \ 25 u.d = (flt); \ 26 u.ieee_nan.mantissa0 = 0; \ 27 u.ieee_nan.mantissa1 = 0; \ 28 u.ieee_nan.mantissa2 = (mant) >> 32; \ 29 u.ieee_nan.mantissa3 = (mant); \ 30 if ((u.ieee.mantissa0 | u.ieee.mantissa1 \ 31 | u.ieee.mantissa2 | u.ieee.mantissa3) != 0) \ 32 (flt) = u.d; \
|
A D | mpn2float128.c | 33 union ieee854_float128 u; in __mpn_construct_float128() local 35 u.ieee.negative = sign; in __mpn_construct_float128() 36 u.ieee.exponent = expt + IEEE854_FLOAT128_BIAS; in __mpn_construct_float128() 38 u.ieee.mantissa3 = frac_ptr[0]; in __mpn_construct_float128() 39 u.ieee.mantissa2 = frac_ptr[1]; in __mpn_construct_float128() 40 u.ieee.mantissa1 = frac_ptr[2]; in __mpn_construct_float128() 41 u.ieee.mantissa0 = frac_ptr[3] & (((mp_limb_t) 1 in __mpn_construct_float128() 44 u.ieee.mantissa3 = frac_ptr[0] & (((mp_limb_t) 1 << 32) - 1); in __mpn_construct_float128() 45 u.ieee.mantissa2 = frac_ptr[0] >> 32; in __mpn_construct_float128() 47 u.ieee.mantissa0 = (frac_ptr[1] >> 32) & (((mp_limb_t) 1 in __mpn_construct_float128() [all …]
|
/sysdeps/ieee754/ldbl-128/ |
A D | math-nan-payload-ldouble.h | 22 union ieee854_long_double u; \ 23 u.d = (flt); \ 24 u.ieee_nan.mantissa0 = 0; \ 25 u.ieee_nan.mantissa1 = 0; \ 26 u.ieee_nan.mantissa2 = (mant) >> 32; \ 27 u.ieee_nan.mantissa3 = (mant); \ 28 if ((u.ieee.mantissa0 | u.ieee.mantissa1 \ 29 | u.ieee.mantissa2 | u.ieee.mantissa3) != 0) \ 30 (flt) = u.d; \
|
A D | mpn2ldbl.c | 31 union ieee854_long_double u; in __mpn_construct_long_double() local 33 u.ieee.negative = sign; in __mpn_construct_long_double() 34 u.ieee.exponent = expt + IEEE854_LONG_DOUBLE_BIAS; in __mpn_construct_long_double() 36 u.ieee.mantissa3 = frac_ptr[0]; in __mpn_construct_long_double() 37 u.ieee.mantissa2 = frac_ptr[1]; in __mpn_construct_long_double() 38 u.ieee.mantissa1 = frac_ptr[2]; in __mpn_construct_long_double() 39 u.ieee.mantissa0 = frac_ptr[3] & (((mp_limb_t) 1 in __mpn_construct_long_double() 42 u.ieee.mantissa3 = frac_ptr[0] & (((mp_limb_t) 1 << 32) - 1); in __mpn_construct_long_double() 43 u.ieee.mantissa2 = frac_ptr[0] >> 32; in __mpn_construct_long_double() 45 u.ieee.mantissa0 = (frac_ptr[1] >> 32) & (((mp_limb_t) 1 in __mpn_construct_long_double() [all …]
|
A D | e_atanhl.c | 48 ieee854_long_double_shape_type u; in __ieee754_atanhl() local 50 u.value = x; in __ieee754_atanhl() 51 jx = u.parts32.w0; in __ieee754_atanhl() 53 u.parts32.w0 = ix; in __ieee754_atanhl() 56 if (u.value == one) in __ieee754_atanhl() 68 t = u.value+u.value; in __ieee754_atanhl() 69 t = 0.5*__log1pl(t+t*u.value/(one-u.value)); in __ieee754_atanhl() 71 t = 0.5*__log1pl((u.value+u.value)/(one-u.value)); in __ieee754_atanhl()
|
A D | s_atanl.c | 180 _Float128 t, u, p, q; in __atanl() local 236 u = L(0.125) * k; in __atanl() 238 t = (x - u) / (1.0 + x * u); in __atanl() 242 u = t * t; in __atanl() 243 p = ((((p4 * u) + p3) * u + p2) * u + p1) * u + p0; in __atanl() 244 q = ((((u + q4) * u + q3) * u + q2) * u + q1) * u + q0; in __atanl() 245 u = t * u * p / q + t; in __atanl() 248 u = atantbl[k] + u; in __atanl() 250 return (-u); in __atanl() 252 return u; in __atanl()
|
A D | s_fmal.c | 45 u.d = x; in __fmal() 161 if (u.ieee.exponent) in __fmal() 164 u.d *= L(0x1p113); in __fmal() 184 x = u.d; in __fmal() 234 u.d = a2 + m2; in __fmal() 238 if ((u.ieee.mantissa3 & 1) == 0 && u.ieee.exponent != 0x7fff) in __fmal() 242 return a1 + u.d; in __fmal() 246 if ((u.ieee.mantissa3 & 1) == 0 && u.ieee.exponent != 0x7fff) in __fmal() 256 v.d = a1 + u.d; in __fmal() 263 asm volatile ("" : "=m" (u) : "m" (u)); in __fmal() [all …]
|
A D | s_nexttowardf.c | 43 float u; in __nexttowardf() local 45 u = math_opt_barrier (x); in __nexttowardf() 46 u = u * u; in __nexttowardf() 47 math_force_eval (u); /* raise underflow flag */ in __nexttowardf() 65 float u = x+x; /* overflow */ in __nexttowardf() local 66 math_force_eval (u); in __nexttowardf() 70 float u = x*x; /* underflow */ in __nexttowardf() local 71 math_force_eval (u); /* raise underflow flag */ in __nexttowardf()
|
/sysdeps/generic/ |
A D | math-nan-payload-double.h | 22 union ieee754_double u; \ 23 u.d = (flt); \ 24 u.ieee_nan.mantissa0 = (mant) >> 32; \ 25 u.ieee_nan.mantissa1 = (mant); \ 26 if ((u.ieee.mantissa0 | u.ieee.mantissa1) != 0) \ 27 (flt) = u.d; \
|
/sysdeps/ieee754/dbl-64/ |
A D | mpn2dbl.c | 30 union ieee754_double u; in __mpn_construct_double() local 32 u.ieee.negative = negative; in __mpn_construct_double() 33 u.ieee.exponent = expt + IEEE754_DOUBLE_BIAS; in __mpn_construct_double() 35 u.ieee.mantissa1 = frac_ptr[0]; in __mpn_construct_double() 36 u.ieee.mantissa0 = frac_ptr[1] & (((mp_limb_t) 1 in __mpn_construct_double() 39 u.ieee.mantissa1 = frac_ptr[0] & (((mp_limb_t) 1 << 32) - 1); in __mpn_construct_double() 40 u.ieee.mantissa0 = (frac_ptr[0] >> 32) & (((mp_limb_t) 1 in __mpn_construct_double() 46 return u.d; in __mpn_construct_double()
|
A D | s_fma.c | 48 u.d = x; in __fma() 160 if (u.ieee.exponent) in __fma() 163 u.d *= 0x1p53; in __fma() 183 x = u.d; in __fma() 233 u.d = a2 + m2; in __fma() 239 v.d = a1 + u.d; in __fma() 249 if ((u.ieee.mantissa1 & 1) == 0 && u.ieee.exponent != 0x7ff) in __fma() 250 u.ieee.mantissa1 |= j; in __fma() 252 return a1 + u.d; in __fma() 256 if ((u.ieee.mantissa1 & 1) == 0 && u.ieee.exponent != 0x7ff) in __fma() [all …]
|
/sysdeps/ieee754/flt-32/ |
A D | s_nextafterf.c | 40 float u; in __nextafterf() local 42 u = math_opt_barrier (x); in __nextafterf() 43 u = u*u; in __nextafterf() 44 math_force_eval (u); /* raise underflow flag */ in __nextafterf() 62 float u = x+x; /* overflow */ in __nextafterf() local 63 math_force_eval (u); in __nextafterf() 67 float u = x*x; /* underflow */ in __nextafterf() local 68 math_force_eval (u); /* raise underflow flag */ in __nextafterf()
|
/sysdeps/i386/fpu/ |
A D | s_nexttowardf.c | 43 float u; in __nexttowardf() local 45 u = math_opt_barrier (x); in __nexttowardf() 46 u = u * u; in __nexttowardf() 47 math_force_eval (u); /* raise underflow flag */ in __nexttowardf() 65 float u = x+x; /* overflow */ in __nexttowardf() local 66 math_force_eval (u); in __nexttowardf() 70 float u = x*x; /* underflow */ in __nexttowardf() local 71 math_force_eval (u); /* raise underflow flag */ in __nexttowardf()
|
/sysdeps/ieee754/ldbl-opt/ |
A D | s_nexttowardfd.c | 47 float u; in __nldbl_nexttowardf() local 49 u = math_opt_barrier (x); in __nldbl_nexttowardf() 50 u = u * u; in __nldbl_nexttowardf() 51 math_force_eval (u); /* raise underflow flag */ in __nldbl_nexttowardf() 67 float u = x+x; /* overflow */ in __nldbl_nexttowardf() local 68 math_force_eval (u); in __nldbl_nexttowardf() 72 float u = x*x; /* underflow */ in __nldbl_nexttowardf() local 73 math_force_eval (u); /* raise underflow flag */ in __nldbl_nexttowardf()
|