/*========================== begin_copyright_notice ============================

Copyright (C) 2020-2021 Intel Corporation

SPDX-License-Identifier: MIT

============================= end_copyright_notice ===========================*/

#include "../imf.h"
#pragma OPENCL FP_CONTRACT OFF
typedef struct
{
    unsigned int _dAbsMask;
    unsigned int _dHalf;
    unsigned int _dSignBit;
    unsigned int _dC1_h;
    unsigned int _dC2;
    unsigned int _dC3;
    unsigned int _dC4;
    unsigned int _dC5;
} __internal_ssinpi_la_data_avx512_t;
static __constant __internal_ssinpi_la_data_avx512_t __internal_ssinpi_la_data_avx512 = {
    0x7FFFFFFFu,
    0x3f000000u,
    0x80000000u,
    0x40490fdbu,
    0xc0a55de2u,
    0x40233479u,
    0xbf1929adu,
    0x3d9f0be5u,
};

typedef struct
{
    unsigned int _sAbsMask;
    unsigned int _sSignBit;
    unsigned int _sReductionRangeVal;
    unsigned int _sRangeVal;
    unsigned int _sPiToRad;

    unsigned int _sA3;
    unsigned int _sA5;
    unsigned int _sA7;
    unsigned int _sA9;
    unsigned int _sA5_FMA;
    unsigned int _sA7_FMA;
    unsigned int _sA9_FMA;

    unsigned int _sRShifter;

} __internal_ssinpi_la_data_t;
static __constant __internal_ssinpi_la_data_t __internal_ssinpi_la_data = {
    0x7fffffffu,
    0x80000000u,
    0x4A800000u,
    0x7f800000u,
    0x40490FDBu,

    0xBE2AAAA6u,
    0x3c088769u,
    0xb94fb7ebu,
    0x362ede96u,
    0x3c088767u,
    0xb94fb6c5u,
    0x362ec33bu,

    0x4B400000u,
};
static __constant int_float __ssinpi_la_c4 = { 0x3d9f0be5u };
static __constant int_float __ssinpi_la_c3 = { 0xbf1929adu };
static __constant int_float __ssinpi_la_c2 = { 0x40233479u };
static __constant int_float __ssinpi_la_c1 = { 0xc0a55de2u };
static __constant int_float __ssinpi_la_c0 = { 0x40490fdau };
static __constant int_float __ssinpi_la_max_norm = { 0x7f7fffffu };

__attribute__((always_inline))
inline int __internal_ssinpi_la_cout (float *a, float *pres)
{
    int nRet = 0;
    float x = *a;
    {
        float fN, fNi, R, R2, poly;
        int iN, sgn;
        int_float xin, res;

        fN = __spirv_ocl_rint(x);

        fNi = -__spirv_ocl_fabs(fN);
        iN = (int) fNi;

        xin.f = x;

        sgn = iN << 31;

        R = x - fN;
        R2 = R * R;

        poly = __spirv_ocl_fma(R2,__ssinpi_la_c4.f,__ssinpi_la_c3.f);
        poly = __spirv_ocl_fma(poly,R2,__ssinpi_la_c2.f);
        poly = __spirv_ocl_fma(poly,R2,__ssinpi_la_c1.f);
        poly = __spirv_ocl_fma(poly,R2,__ssinpi_la_c0.f);

        res.f = R * poly;

        res.w = (res.w == 0) ? (xin.w & 0x80000000) : (res.w ^ sgn);

        *pres = res.f;
    }
    nRet = (__spirv_ocl_fabs(x) > __ssinpi_la_max_norm.f) ? 1 : 0;
    return nRet;

}

float __ocl_svml_sinpif (float a)
{

    float va1;
    float vr1;
    unsigned int vm;

    float r;

    va1 = a;;

    __internal_ssinpi_la_cout (&va1, &vr1);
    r = vr1;;

    return r;

}
