1*f3087befSAndrew Turner /*
2*f3087befSAndrew Turner * Double-precision vector sincos function.
3*f3087befSAndrew Turner *
4*f3087befSAndrew Turner * Copyright (c) 2023-2024, Arm Limited.
5*f3087befSAndrew Turner * SPDX-License-Identifier: MIT OR Apache-2.0 WITH LLVM-exception
6*f3087befSAndrew Turner */
7*f3087befSAndrew Turner
8*f3087befSAndrew Turner /* Define _GNU_SOURCE in order to include sincos declaration. If building
9*f3087befSAndrew Turner pre-GLIBC 2.1, or on a non-GNU conforming system, this routine will need to
10*f3087befSAndrew Turner be linked against the scalar sincosf from math/. */
11*f3087befSAndrew Turner #define _GNU_SOURCE
12*f3087befSAndrew Turner
13*f3087befSAndrew Turner #include "sv_math.h"
14*f3087befSAndrew Turner #include "sv_sincos_common.h"
15*f3087befSAndrew Turner #include "test_defs.h"
16*f3087befSAndrew Turner
17*f3087befSAndrew Turner #include <math.h>
18*f3087befSAndrew Turner
19*f3087befSAndrew Turner /* sincos not available for all scalar libm implementations. */
20*f3087befSAndrew Turner #ifndef __GLIBC__
21*f3087befSAndrew Turner static void
sincos(double x,double * out_sin,double * out_cos)22*f3087befSAndrew Turner sincos (double x, double *out_sin, double *out_cos)
23*f3087befSAndrew Turner {
24*f3087befSAndrew Turner *out_sin = sin (x);
25*f3087befSAndrew Turner *out_cos = cos (x);
26*f3087befSAndrew Turner }
27*f3087befSAndrew Turner #endif
28*f3087befSAndrew Turner
29*f3087befSAndrew Turner static void NOINLINE
special_case(svfloat64_t x,svbool_t special,double * out_sin,double * out_cos)30*f3087befSAndrew Turner special_case (svfloat64_t x, svbool_t special, double *out_sin,
31*f3087befSAndrew Turner double *out_cos)
32*f3087befSAndrew Turner {
33*f3087befSAndrew Turner svbool_t p = svptrue_pat_b64 (SV_VL1);
34*f3087befSAndrew Turner for (int i = 0; i < svcntd (); i++)
35*f3087befSAndrew Turner {
36*f3087befSAndrew Turner if (svptest_any (special, p))
37*f3087befSAndrew Turner sincos (svlastb (p, x), out_sin + i, out_cos + i);
38*f3087befSAndrew Turner p = svpnext_b64 (svptrue_b64 (), p);
39*f3087befSAndrew Turner }
40*f3087befSAndrew Turner }
41*f3087befSAndrew Turner
42*f3087befSAndrew Turner /* Double-precision vector function allowing calculation of both sin and cos in
43*f3087befSAndrew Turner one function call, using shared argument reduction and separate polynomials.
44*f3087befSAndrew Turner Largest observed error is for sin, 3.22 ULP:
45*f3087befSAndrew Turner sv_sincos_sin (0x1.d70eef40f39b1p+12) got -0x1.ffe9537d5dbb7p-3
46*f3087befSAndrew Turner want -0x1.ffe9537d5dbb4p-3. */
47*f3087befSAndrew Turner void
_ZGVsMxvl8l8_sincos(svfloat64_t x,double * out_sin,double * out_cos,svbool_t pg)48*f3087befSAndrew Turner _ZGVsMxvl8l8_sincos (svfloat64_t x, double *out_sin, double *out_cos,
49*f3087befSAndrew Turner svbool_t pg)
50*f3087befSAndrew Turner {
51*f3087befSAndrew Turner const struct sv_sincos_data *d = ptr_barrier (&sv_sincos_data);
52*f3087befSAndrew Turner svbool_t special = check_ge_rangeval (pg, x, d);
53*f3087befSAndrew Turner
54*f3087befSAndrew Turner svfloat64x2_t sc = sv_sincos_inline (pg, x, d);
55*f3087befSAndrew Turner
56*f3087befSAndrew Turner svst1 (pg, out_sin, svget2 (sc, 0));
57*f3087befSAndrew Turner svst1 (pg, out_cos, svget2 (sc, 1));
58*f3087befSAndrew Turner
59*f3087befSAndrew Turner if (unlikely (svptest_any (pg, special)))
60*f3087befSAndrew Turner special_case (x, special, out_sin, out_cos);
61*f3087befSAndrew Turner }
62*f3087befSAndrew Turner
63*f3087befSAndrew Turner TEST_DISABLE_FENV (_ZGVsMxv_sincos_sin)
64*f3087befSAndrew Turner TEST_DISABLE_FENV (_ZGVsMxv_sincos_cos)
65*f3087befSAndrew Turner TEST_ULP (_ZGVsMxv_sincos_sin, 2.73)
66*f3087befSAndrew Turner TEST_ULP (_ZGVsMxv_sincos_cos, 2.73)
67*f3087befSAndrew Turner #define SV_SINCOS_INTERVAL(lo, hi, n) \
68*f3087befSAndrew Turner TEST_SYM_INTERVAL (_ZGVsMxv_sincos_sin, lo, hi, n) \
69*f3087befSAndrew Turner TEST_SYM_INTERVAL (_ZGVsMxv_sincos_cos, lo, hi, n)
70*f3087befSAndrew Turner SV_SINCOS_INTERVAL (0, 0x1p-63, 50000)
71*f3087befSAndrew Turner SV_SINCOS_INTERVAL (0x1p-63, 0x1p23, 500000)
72*f3087befSAndrew Turner SV_SINCOS_INTERVAL (0x1p23, inf, 10000)
73*f3087befSAndrew Turner CLOSE_SVE_ATTR
74