/*========================== 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;
}
