/freebsd/contrib/arm-optimized-routines/math/aarch64/advsimd/ |
H A D | finite_pow.h | 41 return asuint64 (x) >> 52; in top12() 162 return asuint64 (x) >> 63 ? __math_uflow (sign_bias) in exp_inline() 165 double res_uoflow = asuint64 (x) >> 63 ? 0.0 : INFINITY; in exp_inline() 213 return asuint64 (x) >> 63 ? __math_uflow (0) : __math_oflow (0); in exp_nosignbias() 215 return asuint64 (x) >> 63 ? 0.0 : INFINITY; in exp_nosignbias() 266 return 2 * i - 1 >= 2 * asuint64 (INFINITY) - 1; in zeroinfnan() 276 ix = asuint64 (x); in pow_scalar_special_case() 277 iy = asuint64 (y); in pow_scalar_special_case() 291 if (ix == asuint64 (1.0)) in pow_scalar_special_case() 293 if (2 * ix > 2 * asuint64 (INFINITY) in pow_scalar_special_case() [all …]
|
/freebsd/contrib/arm-optimized-routines/math/ |
H A D | pow.c | 31 return asuint64 (x) >> 52; in top12() 194 if (asuint64 (x) >> 63) in exp_inline() 212 ki = asuint64 (kd) >> 16; in exp_inline() 217 ki = asuint64 (kd); in exp_inline() 270 return 2 * i - 1 >= 2 * asuint64 (INFINITY) - 1; in zeroinfnan() 280 ix = asuint64 (x); in pow() 281 iy = asuint64 (y); in pow() 295 if (ix == asuint64 (1.0)) in pow() 297 if (2 * ix > 2 * asuint64 (INFINITY) in pow() 298 || 2 * iy > 2 * asuint64 (INFINITY)) in pow() [all …]
|
H A D | exp2.c | 74 return asuint64 (x) >> 52; in top12() 94 if (asuint64 (x) == asuint64 (-INFINITY)) in exp2() 98 if (!(asuint64 (x) >> 63)) in exp2() 100 else if (asuint64 (x) >= asuint64 (-1075.0)) in exp2() 103 if (2 * asuint64 (x) > 2 * asuint64 (928.0)) in exp2() 111 ki = asuint64 (kd); /* k. */ in exp2()
|
H A D | log2.c | 28 return asuint64 (x) >> 48; in top16() 40 ix = asuint64 (x); in log2() 44 # define LO asuint64 (1.0 - 0x1.5b51p-5) in log2() 45 # define HI asuint64 (1.0 + 0x1.6ab2p-5) in log2() 51 if (WANT_ROUNDING && unlikely (ix == asuint64 (1.0))) in log2() 59 rhi = asdouble (asuint64 (r) & -1ULL << 32); in log2() 82 if (ix == asuint64 (INFINITY)) /* log(inf) == inf. */ in log2() 87 ix = asuint64 (x * 0x1p52); in log2() 114 rhi = asdouble (asuint64 (r) & -1ULL << 32); in log2()
|
H A D | log.c | 28 return asuint64 (x) >> 48; in top16() 40 ix = asuint64 (x); in log() 44 # define LO asuint64 (1.0 - 0x1p-5) in log() 45 # define HI asuint64 (1.0 + 0x1.1p-5) in log() 47 # define LO asuint64 (1.0 - 0x1p-4) in log() 48 # define HI asuint64 (1.0 + 0x1.09p-4) in log() 54 if (WANT_ROUNDING && unlikely (ix == asuint64 (1.0))) in log() 106 if (ix == asuint64 (INFINITY)) /* log(inf) == inf. */ in log() 111 ix = asuint64 (x * 0x1p52); in log()
|
H A D | exp.c | 76 return asuint64 (x) >> 52; in top12() 98 if (asuint64 (x) == asuint64 (-INFINITY)) in exp_inline() 102 if (asuint64 (x) >> 63) in exp_inline() 120 ki = asuint64 (kd) >> 16; in exp_inline() 125 ki = asuint64 (kd); in exp_inline()
|
H A D | exp10.c | 66 uint64_t ix = asuint64 (x); in exp10() 76 return ix == asuint64 (-INFINITY) ? 0.0 : x + 1.0; in exp10() 95 ki = asuint64 (kd); in exp10()
|
H A D | powf.c | 91 ki = asuint64 (kd); in exp2_inline() 194 if (unlikely ((asuint64 (ylogx) >> 47 & 0xffff) in powf() 195 >= asuint64 (126.0 * POWF_SCALE) >> 47)) in powf()
|
/freebsd/contrib/arm-optimized-routines/math/aarch64/ |
H A D | sinpi_3u5.c | 41 double r = asdouble (asuint64 (x) & ~0x8000000000000000); in arm_math_sinpi() 42 uint64_t sign = asuint64 (x) & 0x8000000000000000; in arm_math_sinpi() 63 uint64_t iy = ((m & 1) << 63) ^ asuint64 (1.0); in arm_math_sinpi() 69 sign ^= (asuint64 (n) << 63); in arm_math_sinpi() 84 return asdouble (asuint64 (y) ^ sign); in arm_math_sinpi()
|
H A D | sincospi_4u.c | 34 return (asuint64 (x) >> 52) & 0x7ff; in abstop12() 75 uint64_t sign = asuint64 (x) & 0x8000000000000000; in arm_math_sincospi() 106 = asdouble (asuint64 (fma (-4 * ar2_s, ar_s, ss * ar_s)) ^ ss_sign); in arm_math_sincospi() 108 = asdouble (asuint64 (fma (-4 * ar2_c, ar_c, cc * ar_c)) ^ cc_sign); in arm_math_sincospi() 126 uint64_t ax = asuint64 (x) & 0x7fffffffffffffff; in arm_math_sincospi()
|
H A D | tanpi_2u5.c | 59 uint64_t xabs_12 = asuint64 (x) >> 52 & 0x7ff; in arm_math_tanpi() 70 uint64_t x_sign = asuint64 (x) & SIGN_MASK; in arm_math_tanpi() 117 uint64_t sign = asuint64 (x_reduced) & SIGN_MASK; in arm_math_tanpi() 118 r_x = asdouble (asuint64 (r_x) ^ sign); in arm_math_tanpi() 121 uint64_t is_finite = asuint64 (abs_x_reduced); in arm_math_tanpi()
|
H A D | cospi_3u5.c | 37 double ax = asdouble (asuint64 (x) & ~0x8000000000000000); in arm_math_cospi() 62 uint64_t sign = asuint64 (n) << 63; in arm_math_cospi() 81 return asdouble (asuint64 (y) ^ sign); in arm_math_cospi()
|
/freebsd/contrib/arm-optimized-routines/math/aarch64/experimental/ |
H A D | log10_2u.c | 22 #define LO asuint64 (1.0 - 0x1p-4) 23 #define HI asuint64 (1.0 + 0x1.09p-4) 29 return asuint64 (x) >> 48; in top16() 49 ix = asuint64 (x); in log10() 56 if (WANT_ROUNDING && unlikely (ix == asuint64 (1.0))) in log10() 86 if (ix == asuint64 (INFINITY)) /* log10(inf) == inf. */ in log10() 91 ix = asuint64 (x * 0x1p52); in log10()
|
H A D | erf_2u5.c | 48 uint64_t ix = asuint64 (x); in arm_math_erf() 61 uint64_t i = asuint64 (z) - asuint64 (Shift); in arm_math_erf() 87 return asdouble (asuint64 (y) | sign); in arm_math_erf() 95 return asdouble (sign | asuint64 (1.0)); in arm_math_erf()
|
H A D | exp_inline.h | 77 return asuint64 (x) >> 52; in top12() 99 if (asuint64 (x) == asuint64 (-INFINITY)) in exp_inline() 103 if (asuint64 (x) >> 63) in exp_inline() 121 ki = asuint64 (kd) >> 16; in exp_inline() 126 ki = asuint64 (kd); in exp_inline()
|
H A D | asinh_2u5.c | 45 uint64_t ix = asuint64 (x); in asinh() 63 return asdouble (asuint64 (y) | sign); in asinh() 68 return asdouble (asuint64 (log (ax) + Ln2) | sign); in asinh() 71 return asdouble (asuint64 (log (ax + sqrt (ax * ax + 1))) | sign); in asinh()
|
H A D | erfc_1u8.c | 66 uint64_t ix = asuint64 (x); in erfc() 83 if (unlikely (ia < asuint64 (0x1p-511))) in erfc() 89 uint64_t i = asuint64 (z) - asuint64 (Shift); in erfc() 127 double fac = asdouble (asuint64 (0x1p-128) | sign); in erfc()
|
H A D | atan2_2u5.c | 31 uint64_t fi = asuint64 (f); in biased_exponent() 42 uint64_t ix = asuint64 (x); in atan2() 43 uint64_t iy = asuint64 (y); in atan2() 149 return asdouble (asuint64 (ret) ^ sign_xy); in atan2()
|
H A D | atan_2u5.c | 26 uint64_t ix = asuint64 (x); in atan() 40 return asdouble (asuint64 (PiOver2) ^ sign); in atan() 63 return asdouble (asuint64 (y) ^ sign); in atan()
|
H A D | cbrt_2u.c | 30 uint64_t ix = asuint64 (x); in cbrt() 65 return asdouble (asuint64 (ldexp (a * T (2 + e % 3), e / 3)) | sign); in cbrt()
|
H A D | asin_3u.c | 62 uint64_t ix = asuint64 (x); in asin() 96 return asdouble (asuint64 (y) | sign); in asin()
|
H A D | acos_2u.c | 55 uint64_t ix = asuint64 (x); in acos() 88 return PiOver2 - asdouble (asuint64 (p) | sign); in acos()
|
H A D | atanh_3u.c | 31 uint64_t mi = asuint64 (m); in log1p_inline() 64 uint64_t ix = asuint64 (x); in atanh()
|
H A D | log1p_2u.c | 43 uint64_t ix = asuint64 (x); in log1p() 102 uint64_t mi = asuint64 (m); in log1p()
|
/freebsd/contrib/arm-optimized-routines/math/aarch64/sve/ |
H A D | pow.c | 130 return 2 * i - 1 >= 2 * asuint64 (INFINITY) - 1; in zeroinfnan() 138 2 * asuint64 (INFINITY) - 1); in sv_zeroinfnan() 347 uint64_t ix = asuint64 (x); in pow_sc() 348 uint64_t iy = asuint64 (y); in pow_sc() 354 if (ix == asuint64 (1.0)) in pow_sc() 356 if (2 * ix > 2 * asuint64 (INFINITY) || 2 * iy > 2 * asuint64 (INFINITY)) in pow_sc() 358 if (2 * ix == 2 * asuint64 (1.0)) in pow_sc() 360 if ((2 * ix < 2 * asuint64 (1.0)) == !(iy >> 63)) in pow_sc()
|