67 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#include <clc/utils.h>10 11#if __CLC_FPSIZE == 3212_CLC_OVERLOAD _CLC_DEF __CLC_GENTYPE13__clc_frexp(__CLC_GENTYPE x, __CLC_ADDRESS_SPACE __CLC_INTN *ep) {14 __CLC_INTN i = __CLC_AS_INTN(x);15 __CLC_INTN ai = i & EXSIGNBIT_SP32;16 17 // Scale subnormal by 2^26 without multiplying18 __CLC_INTN is_subnormal = ai > 0 && ai < 0x00800000;19 __CLC_GENTYPE s = __CLC_AS_GENTYPE(ai | 0x0d800000) - 0x1.0p-100f;20 ai = __clc_select(ai, __CLC_AS_INTN(s), is_subnormal);21 __CLC_INTN e = (ai >> EXPSHIFTBITS_SP32) - 126 -22 __clc_select((__CLC_INTN)0, (__CLC_INTN)26, is_subnormal);23 24 i = (i & (__CLC_INTN)SIGNBIT_SP32) | (__CLC_INTN)HALFEXPBITS_SP32 |25 (ai & (__CLC_INTN)MANTBITS_SP32);26 27 __CLC_INTN is_inf_nan_or_zero =28 ai == (__CLC_INTN)0 || __clc_isinf(x) || __clc_isnan(x);29 *ep = __clc_select(e, (__CLC_INTN)0, is_inf_nan_or_zero);30 return __clc_select(__CLC_AS_GENTYPE(i), x, is_inf_nan_or_zero);31}32#endif33 34#if __CLC_FPSIZE == 1635 36_CLC_OVERLOAD _CLC_DEF __CLC_GENTYPE37__clc_frexp(__CLC_GENTYPE x, __CLC_ADDRESS_SPACE __CLC_INTN *ep) {38 return __CLC_CONVERT_HALFN(__clc_frexp(__CLC_CONVERT_FLOATN(x), ep));39}40 41#endif42 43#if __CLC_FPSIZE == 6444 45_CLC_OVERLOAD _CLC_DEF __CLC_GENTYPE46__clc_frexp(__CLC_GENTYPE x, __CLC_ADDRESS_SPACE __CLC_INTN *ep) {47 __CLC_LONGN i = __CLC_AS_LONGN(x);48 __CLC_LONGN ai = i & EXSIGNBIT_DP64;49 50 // Scale subnormal by 2^54 without multiplying51 __CLC_LONGN is_subnormal = ai > 0 && ai < 0x0010000000000000L;52 __CLC_GENTYPE s = __CLC_AS_GENTYPE(ai | 0x0370000000000000L) - 0x1.0p-968;53 ai = __clc_select(ai, __CLC_AS_LONGN(s), is_subnormal);54 __CLC_LONGN e = (ai >> EXPSHIFTBITS_DP64) - (__CLC_LONGN)1022 -55 __clc_select((__CLC_LONGN)0, (__CLC_LONGN)54, is_subnormal);56 57 i = (i & (__CLC_LONGN)SIGNBIT_DP64) | (__CLC_LONGN)HALFEXPBITS_DP64 |58 (ai & (__CLC_LONGN)MANTBITS_DP64);59 60 __CLC_LONGN is_inf_nan_or_zero =61 x == __CLC_FP_LIT(0.0) || __clc_isinf(x) || __clc_isnan(x);62 *ep = __CLC_CONVERT_INTN(__clc_select(e, 0L, is_inf_nan_or_zero));63 return __clc_select(__CLC_AS_GENTYPE(i), x, is_inf_nan_or_zero);64}65 66#endif67