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
|
/*========================== begin_copyright_notice ============================
Copyright (C) 2021 Intel Corporation
SPDX-License-Identifier: MIT
============================= end_copyright_notice ===========================*/
#include "../imf.h"
#pragma OPENCL FP_CONTRACT OFF
typedef struct
{
unsigned int _AbsMask;
unsigned int _Half;
unsigned int _SgnMask;
unsigned int _zero;
unsigned int _c4;
unsigned int _c3;
unsigned int _c2;
unsigned int _c1;
unsigned int _c0;
} __internal_ssincospi_data_t;
static __constant __internal_ssincospi_data_t __internal_ssincospi_data = {
0x7FFFFFFFu,
0x3F000000u,
0x80000000u,
0x00000000u,
0x3d9f0000u,
0xbf1929adu,
0x40233479u,
0xc0a55de2u,
0x40490fdbu,
};
void __ocl_svml_sincospif (float a, __private float *b, __private float *c)
{
float va1;
float vr1;
float vr2;
va1 = a;
{
float AbsMask;
float fN;
unsigned int iN;
float Rs;
float Rs2;
float Rc;
float Rc2;
float aRs;
float Half;
float c4;
float c3;
float c2;
float c1;
float c0;
float sgn_N;
float spoly;
float cpoly;
float zero;
float sgn_x;
float sgn_c;
float sgn_s;
float SgnMask;
float Rs_msk;
float Rc_msk;
AbsMask = as_float (__internal_ssincospi_data._AbsMask);
fN = SPIRV_OCL_BUILTIN(rint, _f32, ) (va1);
Rs = (va1 - fN);
iN = ((int) ((-SPIRV_OCL_BUILTIN(fabs, _f64, ) (fN))));
Half = as_float (__internal_ssincospi_data._Half);
aRs = as_float ((as_uint (Rs) & as_uint (AbsMask)));
Rc = (Half - aRs);
Rs2 = (Rs * Rs);
Rc2 = (Rc * Rc);
iN = ((unsigned int) (iN) << (31));
sgn_N = as_float (iN);
c4 = as_float (__internal_ssincospi_data._c4);
c3 = as_float (__internal_ssincospi_data._c3);
c2 = as_float (__internal_ssincospi_data._c2);
spoly = SPIRV_OCL_BUILTIN(fma, _f32_f32_f32, ) (c4, Rs2, c3);
cpoly = SPIRV_OCL_BUILTIN(fma, _f32_f32_f32, ) (c4, Rc2, c3);
spoly = SPIRV_OCL_BUILTIN(fma, _f32_f32_f32, ) (spoly, Rs2, c2);
cpoly = SPIRV_OCL_BUILTIN(fma, _f32_f32_f32, ) (cpoly, Rc2, c2);
zero = as_float (__internal_ssincospi_data._zero);
SgnMask = as_float (__internal_ssincospi_data._SgnMask);
Rs_msk = as_float (((unsigned int) (-(signed int) (!(Rs == zero)))));
sgn_x = as_float ((as_uint (va1) & as_uint (SgnMask)));
Rc_msk = as_float (((unsigned int) (-(signed int) (!(Rc == zero)))));
sgn_s = as_float ((((~as_uint (Rs_msk)) & as_uint (sgn_x)) | (as_uint (Rs_msk) & as_uint (sgn_N))));
sgn_c = as_float ((as_uint (sgn_N) & as_uint (Rc_msk)));
c1 = as_float (__internal_ssincospi_data._c1);
c0 = as_float (__internal_ssincospi_data._c0);
spoly = SPIRV_OCL_BUILTIN(fma, _f32_f32_f32, ) (spoly, Rs2, c1);
cpoly = SPIRV_OCL_BUILTIN(fma, _f32_f32_f32, ) (cpoly, Rc2, c1);
spoly = SPIRV_OCL_BUILTIN(fma, _f32_f32_f32, ) (spoly, Rs2, c0);
cpoly = SPIRV_OCL_BUILTIN(fma, _f32_f32_f32, ) (cpoly, Rc2, c0);
Rs = as_float ((as_uint (Rs) ^ as_uint (sgn_s)));
Rc = as_float ((as_uint (Rc) ^ as_uint (sgn_c)));
vr1 = (spoly * Rs);
vr2 = (cpoly * Rc);
}
((float *) b)[0] = vr1;
((float *) c)[0] = vr2;
return;
}
|
