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
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
|
/* Configuration for libmvec routines.
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/>. */
#ifndef _VECMATH_CONFIG_H
#define _VECMATH_CONFIG_H
#include <math_private.h>
/* Deprecated config option from Arm Optimized Routines which ensures
fp exceptions are correctly triggered. This is not intended to be
supported in GLIBC, however we keep it for ease of development. */
#define WANT_SIMD_EXCEPT 0
/* Return ptr but hide its value from the compiler so accesses through it
cannot be optimized based on the contents. */
#define ptr_barrier(ptr) \
({ \
__typeof (ptr) __ptr = (ptr); \
__asm("" : "+r"(__ptr)); \
__ptr; \
})
#define V_LOG_POLY_ORDER 6
#define V_LOG_TABLE_BITS 7
extern const struct v_log_data
{
/* Shared data for vector log and log-derived routines (e.g. asinh). */
double poly[V_LOG_POLY_ORDER - 1];
double ln2;
struct
{
double invc, logc;
} table[1 << V_LOG_TABLE_BITS];
} __v_log_data attribute_hidden;
#define V_EXP_TAIL_TABLE_BITS 8
extern const uint64_t __v_exp_tail_data[1 << V_EXP_TAIL_TABLE_BITS] attribute_hidden;
#define V_EXP_TABLE_BITS 7
extern const uint64_t __v_exp_data[1 << V_EXP_TABLE_BITS] attribute_hidden;
#define V_LOG2_TABLE_BITS 7
extern const struct v_log2_data
{
double poly[5];
double invln2;
struct
{
double invc, log2c;
} table[1 << V_LOG2_TABLE_BITS];
} __v_log2_data attribute_hidden;
#define V_LOG10_TABLE_BITS 7
extern const struct v_log10_data
{
double poly[5];
double invln10, log10_2;
struct
{
double invc, log10c;
} table[1 << V_LOG10_TABLE_BITS];
} __v_log10_data attribute_hidden;
extern const struct erff_data
{
struct
{
float erf, scale;
} tab[513];
} __erff_data attribute_hidden;
extern const struct sv_erff_data
{
float erf[513];
float scale[513];
} __sv_erff_data attribute_hidden;
extern const struct erf_data
{
struct
{
double erf, scale;
} tab[769];
} __erf_data attribute_hidden;
extern const struct sv_erf_data
{
double erf[769];
double scale[769];
} __sv_erf_data attribute_hidden;
extern const struct erfc_data
{
struct
{
double erfc, scale;
} tab[3488];
} __erfc_data attribute_hidden;
extern const struct erfcf_data
{
struct
{
float erfc, scale;
} tab[645];
} __erfcf_data attribute_hidden;
/* Some data for AdvSIMD and SVE pow's internal exp and log. */
#define V_POW_EXP_TABLE_BITS 8
extern const struct v_pow_exp_data
{
double poly[3];
double n_over_ln2, ln2_over_n_hi, ln2_over_n_lo, shift;
uint64_t sbits[1 << V_POW_EXP_TABLE_BITS];
} __v_pow_exp_data attribute_hidden;
#define V_POW_LOG_TABLE_BITS 7
extern const struct v_pow_log_data
{
double poly[7]; /* First coefficient is 1. */
double ln2_hi, ln2_lo;
double invc[1 << V_POW_LOG_TABLE_BITS];
double logc[1 << V_POW_LOG_TABLE_BITS];
double logctail[1 << V_POW_LOG_TABLE_BITS];
} __v_pow_log_data attribute_hidden;
/* Some data for SVE powf's internal exp and log. */
#define V_POWF_EXP2_TABLE_BITS 5
#define V_POWF_EXP2_N (1 << V_POWF_EXP2_TABLE_BITS)
#define V_POWF_LOG2_TABLE_BITS 5
#define V_POWF_LOG2_N (1 << V_POWF_LOG2_TABLE_BITS)
extern const struct v_powf_data
{
double invc[V_POWF_LOG2_N];
double logc[V_POWF_LOG2_N];
uint64_t scale[V_POWF_EXP2_N];
} __v_powf_data attribute_hidden;
#endif
|
