69 lines · plain
1//===----------------------------------------------------------------------===//2//3// Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions.4// See https://llvm.org/LICENSE.txt for license information.5// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception6//7//===----------------------------------------------------------------------===//8 9#if __CLC_FPSIZE == 3210 11_CLC_OVERLOAD _CLC_DEF __CLC_FLOATN __clc_sin(__CLC_FLOATN x) {12 __CLC_FLOATN absx = __clc_fabs(x);13 14 __CLC_FLOATN r0, r1;15 __CLC_INTN regn = __clc_argReductionS(&r0, &r1, absx);16 17 __CLC_FLOATN ss = __clc_sinf_piby4(r0, r1);18 __CLC_FLOATN cc = __clc_cosf_piby4(r0, r1);19 20 __CLC_FLOATN s = (regn & 1) != 0 ? cc : ss;21 s = __CLC_AS_FLOATN(__CLC_AS_INTN(s) ^ ((regn > 1) << 31) ^22 (__CLC_AS_INTN(x) ^ __CLC_AS_INTN(absx)));23 24 s = __clc_select(s, __CLC_GENTYPE_NAN, __clc_isnan(x) || __clc_isinf(x));25 26 // Subnormals27 s = x == 0.0f ? x : s;28 29 return s;30}31 32#elif __CLC_FPSIZE == 1633 34_CLC_OVERLOAD _CLC_DEF __CLC_GENTYPE __clc_sin(__CLC_GENTYPE x) {35 return __CLC_CONVERT_GENTYPE(__clc_sin(__CLC_CONVERT_FLOATN(x)));36}37 38#elif __CLC_FPSIZE == 6439 40_CLC_OVERLOAD _CLC_DEF __CLC_GENTYPE __clc_sin(__CLC_GENTYPE x) {41 __CLC_GENTYPE absx = __clc_fabs(x);42 43 __CLC_BIT_INTN is_medium = absx < 0x1.0p+47;44 45 __CLC_INTN regn_m, regn_l;46 __CLC_GENTYPE r_m, r_l, rr_m, rr_l;47 48 __clc_remainder_piby2_medium(absx, &r_m, &rr_m, ®n_m);49 __clc_remainder_piby2_large(absx, &r_l, &rr_l, ®n_l);50 51 __CLC_GENTYPE r = is_medium ? r_m : r_l;52 __CLC_GENTYPE rr = is_medium ? rr_m : rr_l;53 __CLC_INTN regn = __CLC_CONVERT_INTN(is_medium) ? regn_m : regn_l;54 55 __CLC_GENTYPE sinval, cosval;56 __clc_sincos_piby4(r, rr, &sinval, &cosval);57 58 __CLC_LONGN s =59 __CLC_AS_LONGN(__CLC_CONVERT_BIT_INTN((regn & 1) != 0) ? cosval : sinval);60 61 s ^= (__CLC_CONVERT_BIT_INTN(regn > 1) << 63) ^62 (__CLC_CONVERT_BIT_INTN(x < 0.0) << 63);63 64 return __clc_isinf(x) | __clc_isnan(x) ? __CLC_GENTYPE_NAN65 : __CLC_AS_GENTYPE(s);66}67 68#endif69