62 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_DEF _CLC_OVERLOAD __CLC_GENTYPE __clc_tan(__CLC_GENTYPE x) {12 __CLC_GENTYPE absx = __clc_fabs(x);13 __CLC_UINTN x_signbit = __CLC_AS_UINTN(x) & SIGNBIT_SP32;14 15 __CLC_GENTYPE r0, r1;16 __CLC_INTN regn = __clc_argReductionS(&r0, &r1, absx);17 18 __CLC_GENTYPE t = __clc_tanf_piby4(r0 + r1, regn);19 t = __CLC_AS_GENTYPE(__CLC_AS_UINTN(t) ^ x_signbit);20 21 t = __clc_select(t, __CLC_GENTYPE_NAN, __clc_isnan(x) || __clc_isinf(x));22 // Take care of subnormals23 t = (x == 0.0f) ? x : t;24 return t;25}26 27#elif __CLC_FPSIZE == 6428 29_CLC_DEF _CLC_OVERLOAD __CLC_GENTYPE __clc_tan(__CLC_GENTYPE x) {30 __CLC_GENTYPE y = __clc_fabs(x);31 32 __CLC_BIT_INTN is_medium = y < 0x1.0p+30;33 34 __CLC_INTN regn_m, regn_l;35 __CLC_GENTYPE r_m, r_l, rr_m, rr_l;36 37 __clc_remainder_piby2_medium(y, &r_m, &rr_m, ®n_m);38 __clc_remainder_piby2_large(y, &r_l, &rr_l, ®n_l);39 40 __CLC_GENTYPE r = is_medium ? r_m : r_l;41 __CLC_GENTYPE rr = is_medium ? rr_m : rr_l;42 __CLC_INTN regn = __CLC_CONVERT_INTN(is_medium) ? regn_m : regn_l;43 44 __CLC_GENTYPE lead, tail;45 __clc_tan_piby4(r, rr, &lead, &tail);46 47 __CLC_LONGN t =48 __CLC_AS_LONGN(__CLC_CONVERT_BIT_INTN((regn & 1) != 0) ? tail : lead);49 t ^= __CLC_CONVERT_BIT_INTN(x < 0.0) << 63;50 51 return __clc_isnan(x) || __clc_isinf(x) ? __CLC_GENTYPE_NAN52 : __CLC_AS_GENTYPE(t);53}54 55#elif __CLC_FPSIZE == 1656 57_CLC_OVERLOAD _CLC_DEF __CLC_GENTYPE __clc_tan(__CLC_GENTYPE x) {58 return __CLC_CONVERT_GENTYPE(__clc_tan(__CLC_CONVERT_FLOATN(x)));59}60 61#endif62