1*f3087befSAndrew Turner /*
2*f3087befSAndrew Turner * Single-precision SVE hypot(x) 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 #include "sv_math.h"
9*f3087befSAndrew Turner #include "test_sig.h"
10*f3087befSAndrew Turner #include "test_defs.h"
11*f3087befSAndrew Turner
12*f3087befSAndrew Turner #define TinyBound 0x0c800000 /* asuint (0x1p-102). */
13*f3087befSAndrew Turner #define Thres 0x73000000 /* 0x70000000 - TinyBound. */
14*f3087befSAndrew Turner
15*f3087befSAndrew Turner static svfloat32_t NOINLINE
special_case(svfloat32_t sqsum,svfloat32_t x,svfloat32_t y,svbool_t pg,svbool_t special)16*f3087befSAndrew Turner special_case (svfloat32_t sqsum, svfloat32_t x, svfloat32_t y, svbool_t pg,
17*f3087befSAndrew Turner svbool_t special)
18*f3087befSAndrew Turner {
19*f3087befSAndrew Turner return sv_call2_f32 (hypotf, x, y, svsqrt_x (pg, sqsum), special);
20*f3087befSAndrew Turner }
21*f3087befSAndrew Turner
22*f3087befSAndrew Turner /* SVE implementation of single-precision hypot.
23*f3087befSAndrew Turner Maximum error observed is 1.21 ULP:
24*f3087befSAndrew Turner _ZGVsMxvv_hypotf (0x1.6a213cp-19, -0x1.32b982p-26) got 0x1.6a2346p-19
25*f3087befSAndrew Turner want 0x1.6a2344p-19. */
SV_NAME_F2(hypot)26*f3087befSAndrew Turner svfloat32_t SV_NAME_F2 (hypot) (svfloat32_t x, svfloat32_t y,
27*f3087befSAndrew Turner const svbool_t pg)
28*f3087befSAndrew Turner {
29*f3087befSAndrew Turner svfloat32_t sqsum = svmla_x (pg, svmul_x (pg, x, x), y, y);
30*f3087befSAndrew Turner
31*f3087befSAndrew Turner svbool_t special = svcmpge (
32*f3087befSAndrew Turner pg, svsub_x (pg, svreinterpret_u32 (sqsum), TinyBound), Thres);
33*f3087befSAndrew Turner
34*f3087befSAndrew Turner if (unlikely (svptest_any (pg, special)))
35*f3087befSAndrew Turner return special_case (sqsum, x, y, pg, special);
36*f3087befSAndrew Turner
37*f3087befSAndrew Turner return svsqrt_x (pg, sqsum);
38*f3087befSAndrew Turner }
39*f3087befSAndrew Turner
40*f3087befSAndrew Turner TEST_SIG (SV, F, 2, hypot, -10.0, 10.0)
41*f3087befSAndrew Turner TEST_ULP (SV_NAME_F2 (hypot), 0.71)
42*f3087befSAndrew Turner TEST_DISABLE_FENV (SV_NAME_F2 (hypot))
43*f3087befSAndrew Turner TEST_INTERVAL2 (SV_NAME_F2 (hypot), 0, inf, 0, inf, 10000)
44*f3087befSAndrew Turner TEST_INTERVAL2 (SV_NAME_F2 (hypot), 0, inf, -0, -inf, 10000)
45*f3087befSAndrew Turner TEST_INTERVAL2 (SV_NAME_F2 (hypot), -0, -inf, 0, inf, 10000)
46*f3087befSAndrew Turner TEST_INTERVAL2 (SV_NAME_F2 (hypot), -0, -inf, -0, -inf, 10000)
47*f3087befSAndrew Turner CLOSE_SVE_ATTR
48