| /src/contrib/arm-optimized-routines/math/aarch64/ |
| H A D | sinpi_3u5.c | 41 double r = asdouble (asuint64 (x) & ~0x8000000000000000); in arm_math_sinpi() 47 return asdouble (sign); in arm_math_sinpi() 52 return asdouble (sign); in arm_math_sinpi() 64 return asdouble (sign ^ iy); in arm_math_sinpi() 84 return asdouble (asuint64 (y) ^ sign); in arm_math_sinpi()
|
| H A D | sincospi_4u.c | 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() 138 *out_sin = asdouble (sign); in arm_math_sincospi()
|
| H A D | cospi_3u5.c | 37 double ax = asdouble (asuint64 (x) & ~0x8000000000000000); in arm_math_cospi() 81 return asdouble (asuint64 (y) ^ sign); in arm_math_cospi()
|
| H A D | tanpi_2u5.c | 71 return asdouble (x_sign); in arm_math_tanpi() 118 r_x = asdouble (asuint64 (r_x) ^ sign); in arm_math_tanpi()
|
| /src/contrib/arm-optimized-routines/math/aarch64/experimental/ |
| H A D | asinh_2u5.c | 47 double ax = asdouble (ia); 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 | atan_2u5.c | 40 return asdouble (asuint64 (PiOver2) ^ sign); in atan() 57 az = asdouble (ia); in atan() 63 return asdouble (asuint64 (y) ^ sign); in atan()
|
| H A D | erfc_1u8.c | 68 double a = asdouble (ia); in erfc() 73 return asdouble (sign >> 1) + 1.0 / x; /* Special cases. */ in erfc() 126 double off = asdouble (sign >> 1); in erfc() 127 double fac = asdouble (asuint64 (0x1p-128) | sign); in erfc()
|
| H A D | atanh_3u.c | 39 double f = asdouble (u_red) - 1; in log1p_inline() 74 double halfsign = asdouble (Half | sign); in atanh() 75 double ax = asdouble (ia); in atanh()
|
| H A D | exp_inline.h | 44 scale = asdouble (sbits); in exp_inline_special_case() 50 scale = asdouble (sbits); in exp_inline_special_case() 136 tail = asdouble (T[idx]); in exp_inline() 153 scale = asdouble (sbits); in exp_inline()
|
| H A D | erf_2u5.c | 59 double a = asdouble (ia); 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 | atan2_2u5.c | 122 double ax = asdouble (iax); in atan2() 123 double ay = asdouble (iay); in atan2() 149 return asdouble (asuint64 (ret) ^ sign_xy); in atan2()
|
| H A D | cosh_2u.c | 26 double t = exp_inline (asdouble (iax) / 2, 0); in specialcase() 50 double ax = asdouble (iax); in cosh()
|
| H A D | sinh_3u.c | 28 double ax = asdouble (iax); in sinh() 30 double halfsign = asdouble (Half | sign); in sinh()
|
| H A D | cbrt_2u.c | 40 double m = frexp (asdouble (iax), &e); in cbrt() 65 return asdouble (asuint64 (ldexp (a * T (2 + e % 3), e / 3)) | sign); in cbrt()
|
| H A D | asin_3u.c | 65 double ax = asdouble (ia); in asin() 96 return asdouble (asuint64 (y) | sign); in asin()
|
| H A D | tanh_3u.c | 45 double t = asdouble ((uint64_t) (i + 1023) << 52); in expm1_inline() 65 return asdouble (One | sign); in tanh()
|
| H A D | acos_2u.c | 58 double ax = asdouble (ia); in acos() 88 return PiOver2 - asdouble (asuint64 (p) | sign); in acos()
|
| H A D | log1p_2u.c | 66 return __math_invalid (asdouble (ia)); in log1p() 112 double f = asdouble (u_red) - 1; in log1p()
|
| /src/contrib/arm-optimized-routines/math/ |
| H A D | pow.c | 52 z = asdouble (iz); in log_inline() 66 double_t zhi = asdouble ((iz + (1ULL << 31)) & (-1ULL << 32)); in log_inline() 137 scale = asdouble (sbits); in specialcase() 144 scale = asdouble (sbits); in specialcase() 161 y = asdouble (sbits & 0x8000000000000000); in specialcase() 226 tail = asdouble (T[idx]); in exp_inline() 243 scale = asdouble (sbits); in exp_inline() 366 double_t yhi = asdouble (iy & -1ULL << 27); in pow() 368 double_t lhi = asdouble (asuint64 (hi) & -1ULL << 27); in pow()
|
| H A D | exp10.c | 31 scale = asdouble (sbits); in special_case() 38 scale = asdouble (sbits); in special_case() 113 double_t tail = asdouble (__exp_data.tab[i]); in exp10() 129 double_t s = asdouble (sbits); in exp10()
|
| H A D | exp2.c | 41 scale = asdouble (sbits); in specialcase() 47 scale = asdouble (sbits); in specialcase() 117 tail = asdouble (T[idx]); in exp2() 134 scale = asdouble (sbits); in exp2()
|
| H A D | exp.c | 43 scale = asdouble (sbits); in specialcase() 49 scale = asdouble (sbits); in specialcase() 135 tail = asdouble (T[idx]); in exp_inline() 152 scale = asdouble (sbits); in exp_inline()
|
| H A D | log2.c | 59 rhi = asdouble (asuint64 (r) & -1ULL << 32); in log2() 100 z = asdouble (iz); in log2() 114 rhi = asdouble (asuint64 (r) & -1ULL << 32); in log2()
|
| /src/contrib/arm-optimized-routines/math/aarch64/advsimd/ |
| H A D | finite_pow.h | 57 double z = asdouble (iz); in log_inline() 109 scale = asdouble (sbits); in special_case() 116 scale = asdouble (sbits); in special_case() 134 y = asdouble (sbits & 0x8000000000000000); in special_case() 192 double scale = asdouble (sbits); in exp_inline() 239 double scale = asdouble (sbits); in exp_nosignbias()
|
| /src/contrib/arm-optimized-routines/math/test/ |
| H A D | ulp.c | 161 return (struct args_d1){asdouble (next1 (g))}; in next_d1() 169 return (struct args_d2){asdouble (x), asdouble (x2)}; in next_d2() 457 printf ("%s in [%a;%a]", f->name, asdouble (gen->start), in printgen_d1() 458 asdouble (gen->start + gen->len)); in printgen_d1() 463 printf ("%s in [%a;%a] x [%a;%a]", f->name, asdouble (gen->start), in printgen_d2() 464 asdouble (gen->start + gen->len), asdouble (gen->start2), in printgen_d2() 465 asdouble (gen->start2 + gen->len2)); in printgen_d2() 571 #define asfloat_d asdouble
|