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
|
/*========================== begin_copyright_notice ============================
Copyright (C) 2017-2021 Intel Corporation
SPDX-License-Identifier: MIT
============================= end_copyright_notice ===========================*/
// SVML code
static float __ocl_svml_sinf_noLUT( float a )
{
float result = 0.0f;
float sN = __spirv_ocl_fma( a, M_1_PI_F, as_float( 0x4B400000 ) );
uint usN = as_uint( sN ) << 31;
sN = sN - as_float( 0x4B400000 );
float sR = __spirv_ocl_fma( sN, as_float( 0xC0490000 ), a );
sR = __spirv_ocl_fma( sN, as_float( 0xBA7DA000 ), sR );
sR = __spirv_ocl_fma( sN, as_float( 0xB4222000 ), sR );
sR = __spirv_ocl_fma( sN, as_float( 0xACB4611A ), sR );
float sR2 = sR * sR;
sR = as_float( as_uint( sR ) ^ usN );
float sP = as_float( 0x362EDEF8 );
sP = __spirv_ocl_fma( sP, sR2, as_float( 0xB94FB7FF ) );
sP = __spirv_ocl_fma( sP, sR2, as_float( 0x3C08876A ) );
sP = __spirv_ocl_fma( sP, sR2, as_float( 0xBE2AAAA6 ) );
sP = sP * sR2;
result = __spirv_ocl_fma( sP, sR, sR );
return result;
}
|
