1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
|
/* Single-precision vector (SVE) exp function.
Copyright (C) 2023-2024 Free Software Foundation, Inc.
This file is part of the GNU C Library.
The GNU C Library is free software; you can redistribute it and/or
modify it under the terms of the GNU Lesser General Public
License as published by the Free Software Foundation; either
version 2.1 of the License, or (at your option) any later version.
The GNU C Library is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
Lesser General Public License for more details.
You should have received a copy of the GNU Lesser General Public
License along with the GNU C Library; if not, see
<https://www.gnu.org/licenses/>. */
#include "sv_math.h"
#include "sv_expf_inline.h"
/* Roughly 87.3. For x < -Thres, the result is subnormal and not handled
correctly by FEXPA. */
#define Thres 0x1.5d5e2ap+6f
static const struct data
{
struct sv_expf_data d;
float thres;
} data = {
.d = SV_EXPF_DATA,
.thres = Thres,
};
static svfloat32_t NOINLINE
special_case (svfloat32_t x, svbool_t special, const struct sv_expf_data *d)
{
return sv_call_f32 (expf, x, expf_inline (x, svptrue_b32 (), d), special);
}
/* Optimised single-precision SVE exp function.
Worst-case error is 1.04 ulp:
SV_NAME_F1 (exp)(0x1.a8eda4p+1) got 0x1.ba74bcp+4
want 0x1.ba74bap+4. */
svfloat32_t SV_NAME_F1 (exp) (svfloat32_t x, const svbool_t pg)
{
const struct data *d = ptr_barrier (&data);
svbool_t is_special_case = svacgt (pg, x, d->thres);
if (__glibc_unlikely (svptest_any (pg, is_special_case)))
return special_case (x, is_special_case, &d->d);
return expf_inline (x, pg, &d->d);
}
|