/freebsd-current/contrib/arm-optimized-routines/pl/math/ |
H A D | atanf_common.h | 26 float z2 = z * z; local 27 float z4 = z2 * z2; 31 z4, z4 * pairwise_poly_3_f32 (z2, z4, __atanf_poly_data.poly + 4), 32 pairwise_poly_3_f32 (z2, z4, __atanf_poly_data.poly)); 35 return fmaf (y, z2 * az, az) + shift;
|
H A D | atan_common.h | 19 double z2 = z * z; local 20 double x2 = z2 * z2; 23 double y = fma (estrin_11_f64 (z2, x2, x4, x8, __atan_poly_data.poly + 8), 24 x8, estrin_7_f64 (z2, x2, x4, __atan_poly_data.poly)); 27 y = fma (y, z2 * az, az);
|
H A D | acos_2u.c | 69 /* Evaluate polynomial Q(|x|) = z + z * z2 * P(z2) with 70 z2 = x ^ 2 and z = |x| , if |x| < 0.5 71 z2 = (1 - |x|) / 2 and z = sqrt(z2), if |x| >= 0.5. */ 72 double z2 = ax < 0.5 ? x * x : fma (-0.5, ax, 0.5); local 73 double z = ax < 0.5 ? ax : sqrt (z2); 76 double z4 = z2 * z2; 79 double p = estrin_11_f64 (z2, z [all...] |
H A D | asin_3u.c | 77 z2 = x ^ 2 and z = |x| , if |x| < 0.5 78 z2 = (1 - |x|) / 2 and z = sqrt(z2), if |x| >= 0.5. */ 79 double z2 = ax < 0.5 ? x * x : fma (-0.5, ax, 0.5); local 80 double z = ax < 0.5 ? ax : sqrt (z2); 83 double z4 = z2 * z2; 86 double p = estrin_11_f64 (z2, z4, z8, z16, __asin_poly); 88 /* Finalize polynomial: z + z * z2 * P(z2) [all...] |
H A D | sv_asin_3u.c | 58 svfloat64_t z2 = svsel (a_ge_half, svmls_x (pg, sv_f64 (0.5), ax, 0.5), local 60 svfloat64_t z = svsqrt_m (ax, a_ge_half, z2); 63 svfloat64_t z4 = svmul_x (pg, z2, z2); 66 svfloat64_t p = sv_estrin_11_f64_x (pg, z2, z4, z8, z16, d->poly); 67 /* Finalize polynomial: z + z * z2 * P(z2). */ 68 p = svmla_x (pg, z, svmul_x (pg, z, z2), p);
|
H A D | sv_asinf_2u5.c | 53 svfloat32_t z2 = svsel (a_ge_half, svmls_x (pg, sv_f32 (0.5), ax, 0.5), local 55 svfloat32_t z = svsqrt_m (ax, a_ge_half, z2); 58 svfloat32_t p = sv_horner_4_f32_x (pg, z2, d->poly); 59 /* Finalize polynomial: z + z * z2 * P(z2). */ 60 p = svmla_x (pg, z, svmul_x (pg, z, z2), p);
|
H A D | asinh_2u5.c | 60 double z2 = x2 * x2; local 61 double z4 = z2 * z2; 63 double p = estrin_17_f64 (x2, z2, z4, z8, z8 * z8, __asinh_data.poly);
|
H A D | acosf_1u4.c | 72 /* Evaluate polynomial Q(|x|) = z + z * z2 * P(z2) with 73 z2 = x ^ 2 and z = |x| , if |x| < 0.5 74 z2 = (1 - |x|) / 2 and z = sqrt(z2), if |x| >= 0.5. */ 75 float z2 = ax < 0.5 ? x * x : fmaf (-0.5f, ax, 0.5f); local 76 float z = ax < 0.5 ? ax : sqrtf (z2); 79 float p = horner_4_f32 (z2, __asinf_poly); 80 /* Finalize polynomial: z + z * z2 * P(z2) [all...] |
H A D | asinf_2u5.c | 75 z2 = x ^ 2 and z = |x| , if |x| < 0.5 76 z2 = (1 - |x|) / 2 and z = sqrt(z2), if |x| >= 0.5. */ 77 float z2 = ax < 0.5 ? x * x : fmaf (-0.5f, ax, 0.5f); local 78 float z = ax < 0.5 ? ax : sqrtf (z2); 81 float p = horner_4_f32 (z2, __asinf_poly); 82 /* Finalize polynomial: z + z * z2 * P(z2). */ 83 p = fmaf (z * z2, p, z);
|
H A D | sv_atanf_2u9.c | 54 svfloat32_t z2 = svmul_x (pg, z, z); local 55 svfloat32_t z4 = svmul_x (pg, z2, z2); 58 svfloat32_t y = sv_estrin_7_f32_x (pg, z2, z4, z8, d->poly); 61 svfloat32_t z3 = svmul_x (pg, z2, az);
|
H A D | v_tanf_3u5.c | 49 float32x4_t z2 = vmulq_f32 (z, z); local 57 z2 = vbslq_f32 (will_uflow, v_f32 (0), z2); 59 float32x4_t z4 = vmulq_f32 (z2, z2); 60 return v_estrin_5_f32 (z, z2, z4, d->poly); 110 float32x4_t z2 = vmulq_f32 (r, r); local 111 float32x4_t p = eval_poly (z2, d); 112 float32x4_t y = vfmaq_f32 (z, vmulq_f32 (z, z2), p);
|
H A D | sv_atan_2u5.c | 60 svfloat64_t z2 = svmul_x (pg, z, z); local 61 svfloat64_t x2 = svmul_x (pg, z2, z2); 66 = svmla_x (pg, sv_estrin_7_f64_x (pg, z2, x2, x4, d->poly), 67 sv_estrin_11_f64_x (pg, z2, x2, x4, x8, d->poly + 8), x8); 70 svfloat64_t z3 = svmul_x (pg, z2, az);
|
H A D | sv_acos_2u.c | 55 /* Evaluate polynomial Q(x) = z + z * z2 * P(z2) with 56 z2 = x ^ 2 and z = |x| , if |x| < 0.5 57 z2 = (1 - |x|) / 2 and z = sqrt(z2), if |x| >= 0.5. */ 58 svfloat64_t z2 = svsel (a_gt_half, svmls_x (pg, sv_f64 (0.5), ax, 0.5), local 60 svfloat64_t z = svsqrt_m (ax, a_gt_half, z2); 63 svfloat64_t z4 = svmul_x (pg, z2, z2); 66 svfloat64_t p = sv_estrin_11_f64_x (pg, z2, z [all...] |
H A D | v_atanf_3u.c | 85 float32x4_t z2 = vmulq_f32 (z, z); local 86 float32x4_t z4 = vmulq_f32 (z2, z2); 89 v_pairwise_poly_3_f32 (z2, z4, d->poly), z4, 90 vmulq_f32 (z4, v_pairwise_poly_3_f32 (z2, z4, d->poly + 4))); 93 y = vaddq_f32 (vfmaq_f32 (az, y, vmulq_f32 (z2, az)), shift);
|
H A D | v_asin_3u.c | 84 float64x2_t z2 = vbslq_f64 (a_lt_half, vmulq_f64 (x, x), local 86 float64x2_t z = vbslq_f64 (a_lt_half, ax, vsqrtq_f64 (z2)); 89 float64x2_t z4 = vmulq_f64 (z2, z2); 92 float64x2_t p = v_estrin_11_f64 (z2, z4, z8, z16, d->poly); 94 /* Finalize polynomial: z + z * z2 * P(z2). */ 95 p = vfmaq_f64 (z, vmulq_f64 (z, z2), p);
|
H A D | sv_acosf_1u4.c | 52 /* Evaluate polynomial Q(x) = z + z * z2 * P(z2) with 53 z2 = x ^ 2 and z = |x| , if |x| < 0.5 54 z2 = (1 - |x|) / 2 and z = sqrt(z2), if |x| >= 0.5. */ 55 svfloat32_t z2 = svsel (a_gt_half, svmls_x (pg, sv_f32 (0.5), ax, 0.5), local 57 svfloat32_t z = svsqrt_m (ax, a_gt_half, z2); 60 svfloat32_t p = sv_horner_4_f32_x (pg, z2, d->poly); 61 /* Finalize polynomial: z + z * z2 * P(z2) [all...] |
H A D | sv_tanf_3u5.c | 86 svfloat32_t z2 = svmul_x (pg, z, z); local 87 svfloat32_t p01 = svmla_lane (sv_f32 (d->c0), z2, odd_coeffs, 0); 88 svfloat32_t p23 = svmla_lane (sv_f32 (d->c2), z2, odd_coeffs, 1); 89 svfloat32_t p45 = svmla_lane (sv_f32 (d->c4), z2, odd_coeffs, 2); 91 svfloat32_t z4 = svmul_x (pg, z2, z2); 97 svfloat32_t y = svmla_x (pg, z, p, svmul_x (pg, z, z2));
|
H A D | v_atan2_3u.c | 91 float64x2_t z2 = vmulq_f64 (z, z); local 92 float64x2_t x2 = vmulq_f64 (z2, z2); 96 = vfmaq_f64 (v_estrin_7_f64 (z2, x2, x4, data_ptr->poly), 97 v_estrin_11_f64 (z2, x2, x4, x8, data_ptr->poly + 8), x8); 100 ret = vfmaq_f64 (z, ret, vmulq_f64 (z2, z));
|
H A D | v_atan2f_3u.c | 86 float32x4_t z2 = vmulq_f32 (z, z); local 87 float32x4_t z4 = vmulq_f32 (z2, z2); 90 v_pairwise_poly_3_f32 (z2, z4, data_ptr->poly), z4, 91 vmulq_f32 (z4, v_pairwise_poly_3_f32 (z2, z4, data_ptr->poly + 4))); 94 ret = vaddq_f32 (vfmaq_f32 (z, ret, vmulq_f32 (z2, z)), shift);
|
/freebsd-current/sys/contrib/libsodium/src/libsodium/crypto_scalarmult/curve25519/sandy2x/ |
H A D | curve25519_sandy2x.c | 22 #define z2 var[2] macro 45 z_51.v[0] = (z2[1] << 26) + z2[0]; 46 z_51.v[1] = (z2[3] << 26) + z2[2]; 47 z_51.v[2] = (z2[5] << 26) + z2[4]; 48 z_51.v[3] = (z2[7] << 26) + z2[6]; 49 z_51.v[4] = (z2[ 65 #undef z2 macro 68 #define z2 macro [all...] |
H A D | fe51_invert.c | 15 fe51 z2; local 25 /* 2 */ fe51_square(&z2,x); 26 /* 4 */ fe51_square(&t,&z2); 29 /* 11 */ fe51_mul(&z11,&z9,&z2);
|
/freebsd-current/sys/contrib/libsodium/src/libsodium/crypto_scalarmult/curve25519/ref10/ |
H A D | x25519_ref10.c | 60 fe25519 z2; local 80 fe25519_0(z2); 90 fe25519_cswap(z2, z3, swap); 93 fe25519_sub(tmp1, x2, z2); 94 fe25519_add(x2, x2, z2); 95 fe25519_add(z2, x3, z3); 97 fe25519_mul(z2, z2, tmp1); 100 fe25519_add(x3, z3, z2); 101 fe25519_sub(z2, z [all...] |
/freebsd-current/contrib/bearssl/src/ec/ |
H A D | ec_c25519_i31.c | 179 uint32_t x1[10], x2[10], x3[10], z2[10], z3[10]; local 223 * Initialise variables x1, x2, z2, x3 and z3. We set all of them 228 br_i31_zero(z2, C255_P[0]); 229 memcpy(x2, z2, sizeof z2); 253 cswap(z2, z3, swap); 258 print_int_mont("z2", z2); 263 c255_add(a, x2, z2); 265 c255_sub(b, x2, z2); [all...] |
H A D | ec_c25519_i15.c | 188 uint16_t x1[18], x2[18], x3[18], z2[18], z3[18]; local 232 * Initialise variables x1, x2, z2, x3 and z3. We set all of them 237 br_i15_zero(z2, C255_P[0]); 238 memcpy(x2, z2, ILEN); 259 cswap(z2, z3, swap); 264 print_int_mont("z2", z2); 269 c255_add(a, x2, z2); 271 c255_sub(b, x2, z2); 297 c255_mul(z2, C255_A2 [all...] |
/freebsd-current/contrib/bearssl/src/hash/ |
H A D | ghash_ctmul64.c | 37 uint64_t z0, z1, z2, z3; local 49 z2 = (x0 * y2) ^ (x1 * y1) ^ (x2 * y0) ^ (x3 * y3); 53 z2 &= (uint64_t)0x4444444444444444; 55 return z0 | z1 | z2 | z3; 100 uint64_t z0, z1, z2, z0h, z1h, z2h; local 123 z2 = bmul64(y2, h2); 127 z2 ^= z0 ^ z1; 134 v1 = z0h ^ z2;
|