pl/math/sv_asinh_3u0.c

31 __sv_log_inline (svfloat64_t x, const svbool_t pg)
37   svuint64_t tmp = svsub_x (pg, ix, OFF);
39       = svand_x (pg, svlsr_x (pg, tmp, (51 - V_LOG_TABLE_BITS)), (N - 1) << 1);
40   svint64_t k = svasr_x (pg, svreinterpret_s64 (tmp), 52);
41   svuint64_t iz = svsub_x (pg, ix, svand_x (pg, tmp, 0xfffULL << 52));
43   svfloat64_t invc = svld1_gather_index (pg, &__v_log_data.table[0].invc, i);
44   svfloat64_t logc = svld1_gather_index (pg, &__v_log_data.table[0].logc, i);
45   svfloat64_t r = svmla_x (pg, sv_f64 (-1.0), invc, z);
46   svfloat64_t kd = svcvt_f64_x (pg, k);
47   svfloat64_t hi = svmla_x (pg, svadd_x (pg, logc, r), kd, Ln2);
48   svfloat64_t r2 = svmul_x (pg, r, r);
49   svfloat64_t y = svmla_x (pg, sv_f64 (A (2)), r, A (3));
50   svfloat64_t p = svmla_x (pg, sv_f64 (A (0)), r, A (1));
51   y = svmla_x (pg, y, r2, A (4));
52   y = svmla_x (pg, p, r2, y);
53   y = svmla_x (pg, hi, r2, y);
68 svfloat64_t SV_NAME_D1 (asinh) (svfloat64_t x, const svbool_t pg)
71   svuint64_t iax = svbic_x (pg, ix, SignMask);
72   svuint64_t sign = svand_x (pg, ix, SignMask);
74   svuint64_t top12 = svlsr_x (pg, iax, 52);
76   svbool_t ge1 = svcmpge (pg, top12, OneTop);
77   svbool_t special = svcmpge (pg, top12, HugeBound);
82   if (likely (svptest_any (pg, ge1)))
84       svfloat64_t axax = svmul_x (pg, ax, ax);
86 	  svadd_x (pg, ax, svsqrt_x (pg, svadd_x (pg, axax, 1))), pg);
95   if (likely (svptest_any (pg, svnot_z (pg, ge1))))
97       svfloat64_t x2 = svmul_x (pg, ax, ax);
98       svfloat64_t z2 = svmul_x (pg, x2, x2);
99       svfloat64_t z4 = svmul_x (pg, z2, z2);
100       svfloat64_t z8 = svmul_x (pg, z4, z4);
101       svfloat64_t z16 = svmul_x (pg, z8, z8);
103 	  = sv_estrin_17_f64_x (pg, x2, z2, z4, z8, z16, __asinh_data.poly);
104       option_2 = svmla_x (pg, ax, p, svmul_x (pg, x2, ax));
111   y = svreinterpret_f64 (sveor_x (pg, svreinterpret_u64 (y), sign));
113   if (unlikely (svptest_any (pg, special)))