Index: lib/builtins/fixtfsi.c =================================================================== --- /dev/null +++ lib/builtins/fixtfsi.c @@ -0,0 +1,62 @@ +//===-- lib/fixtfsi.c - Quad-precision -> integer conversion ------*- C -*-===// +// +// The LLVM Compiler Infrastructure +// +// This file is dual licensed under the MIT and the University of Illinois Open +// Source Licenses. See LICENSE.TXT for details. +// +//===----------------------------------------------------------------------===// +// +// This file implements quad-precision to integer conversion for the +// compiler-rt library. No range checking is performed; the behavior of this +// conversion is undefined for out of range values in the C standard. +// +//===----------------------------------------------------------------------===// + +// __uint128_t is undefined in 32-bit machine toolchain +#if defined(__SIZEOF_INT128__) && __SIZEOF_INT128__ == 16 + +#define QUAD_PRECISION +#include "fp_lib.h" + +#include "int_lib.h" + +int __fixtfsi(fp_t a) { + + // Break a into sign, exponent, significand + const rep_t aRep = toRep(a); + const rep_t aAbs = aRep & absMask; + const int sign = aRep & signBit ? -1 : 1; + const int exponent = (aAbs >> significandBits) - exponentBias; + const rep_t significand = (aAbs & significandMask) | implicitBit; + + // If 0 < exponent < significandBits, right shift to get the result. + if ((unsigned int)exponent < significandBits) { + // value exceed INT_MAX, reuturn INT_MAX + if (exponent >= sizeof(int) * CHAR_BIT){ + return sign * 0x7fffffff; + } + else { + return sign * (significand >> (significandBits - exponent)); + } + } + // If exponent is negative, the result is zero. + else if (exponent < 0) { + return 0; + } + // If significandBits < exponent, left shift to get the result. This shift + // may end up being larger than the type width, which incurs undefined + // behavior, but the conversion itself is undefined in that case, so + // whatever the compiler decides to do is fine. + else { + // value exceed INT_MAX, reuturn INT_MAX + if (exponent >= sizeof(int) * CHAR_BIT){ + return sign * 0x7fffffff; + } + else { + return sign * (significand << (exponent - significandBits)); + } + } +} + +#endif Index: lib/builtins/fixunstfsi.c =================================================================== --- /dev/null +++ lib/builtins/fixunstfsi.c @@ -0,0 +1,52 @@ +/* ===-- fixunstfsi.c - Implement __fixunstfsi -----------------------------=== + * + * The LLVM Compiler Infrastructure + * + * This file is dual licensed under the MIT and the University of Illinois Open + * Source Licenses. See LICENSE.TXT for details. + * + * ===----------------------------------------------------------------------=== + * + * This file implements __fixunstfsi for the compiler_rt library. + * + * ===----------------------------------------------------------------------=== + */ + +// __uint128_t is undefined in 32-bit machine toolchain +#if defined(__SIZEOF_INT128__) && __SIZEOF_INT128__ == 16 + +#include "int_lib.h" + +/* Returns: convert a to a unsigned int, rounding toward zero. + * Negative values all become zero. + */ + +/* Assumption: long double is a IEEE 128 bit floating point type + * su_int is a 32 bit integral type + * value in long double is representable in su_int or is negative + * (no range checking performed) + */ + +/* seee eeee eeee eeee mmmm mmmm mmmm mmmm | mmmm mmmm mmmm mmmm mmmm mmmm mmmm mmmm + mmmm mmmm mmmm mmmm mmmm mmmm mmmm mmmm | mmmm mmmm mmmm mmmm mmmm mmmm mmmm mmmm */ + +COMPILER_RT_ABI su_int +__fixunstfsi(long double a) +{ + long_double_bits fb; + fb.f = a; + int e = ((fb.u.high.s.high & 0x7FFF0000) >> 16) - 16383; + if (e < 0 || (fb.u.high.s.high & 0x80000000)) + return 0; + // out ou range, return UINT_MAX + if (e >= sizeof(su_int) * CHAR_BIT){ + return 0xffffffff; + } + return ( + 0x80000000u | + ((fb.u.high.s.high & 0x0000FFFF) << 15) | + (fb.u.high.s.low >> 17) + ) >> (31 - e); +} + +#endif Index: test/builtins/Unit/fixtfsi_test.c =================================================================== --- /dev/null +++ test/builtins/Unit/fixtfsi_test.c @@ -0,0 +1,61 @@ +//===--------------- fixtfsi_test.c - Test __fixtfsi ----------------------===// +// +// The LLVM Compiler Infrastructure +// +// This file is dual licensed under the MIT and the University of Illinois Open +// Source Licenses. See LICENSE.TXT for details. +// +//===----------------------------------------------------------------------===// +// +// This file tests __fixtfsi for the compiler_rt library. +// +//===----------------------------------------------------------------------===// + +#if defined(__SIZEOF_INT128__) && __SIZEOF_INT128__ == 16 + +#include +#include +#include + +int __fixtfsi(long double a); + +int test__fixtfsi(long double a) +{ + int x = __fixtfsi(a); + int expected = (int)a; + + if (x != expected) + { + printf("error in test__fixtfsi(%.20Lf) = %d, " + "expected %d\n", a, x, expected); + } + return x != expected; +} + +char assumption_1[sizeof(long double) * CHAR_BIT == 128] = {0}; + +#endif + +int main() +{ +#if defined(__SIZEOF_INT128__) && __SIZEOF_INT128__ == 16 + //if (test__fixtfsi(0x1.0p+16383)) + // return 1; + if (test__fixtfsi(0)) + return 1; + if (test__fixtfsi(0x1.23456789abcdefp+5)) + return 1; + if (test__fixtfsi(0x1.23456789abcdefp-3)) + return 1; + if (test__fixtfsi(0x1.23456789abcdefp+20)) + return 1; + if (test__fixtfsi(0x1.23456789abcdefp+40)) + return 1; + if (test__fixtfsi(0x1.23456789abcdefp+256)) + return 1; +#else + printf("skipped\n"); + +#endif + return 0; +} Index: test/builtins/Unit/fixunstfsi_test.c =================================================================== --- /dev/null +++ test/builtins/Unit/fixunstfsi_test.c @@ -0,0 +1,63 @@ +//===--------------- fixunstfsi_test.c - Test __fixunstfsi ----------------===// +// +// The LLVM Compiler Infrastructure +// +// This file is dual licensed under the MIT and the University of Illinois Open +// Source Licenses. See LICENSE.TXT for details. +// +//===----------------------------------------------------------------------===// +// +// This file tests __fixunstfsi for the compiler_rt library. +// +//===----------------------------------------------------------------------===// + +#if defined(__SIZEOF_INT128__) && __SIZEOF_INT128__ == 16 + +#include +#include +#include + +unsigned int __fixunstfsi(long double a); + +int test__fixunstfsi(long double a) +{ + unsigned int x = __fixunstfsi(a); + unsigned int expected = (unsigned int)a; + + if (x != expected) + { + printf("error in test__fixunstfsi(%.20Lf) = %u, " + "expected %u\n", a, x, expected); + } + return x != expected; +} + +char assumption_1[sizeof(long double) * CHAR_BIT == 128] = {0}; + +#endif + +int main() +{ +#if defined(__SIZEOF_INT128__) && __SIZEOF_INT128__ == 16 + //if (test__fixunstfsi(0x1.0p+16383)) + // return 1; + if (test__fixunstfsi(0)) + return 1; + if (test__fixunstfsi(0x1.23456789abcdefp+5)) + return 1; + if (test__fixunstfsi(0x1.23456789abcdefp-3)) + return 1; + if (test__fixunstfsi(0x1.23456789abcdefp+20)) + return 1; + if (test__fixunstfsi(0x1.23456789abcdefp+40)) + return 1; + if (test__fixunstfsi(0x1.23456789abcdefp+256)) + return 1; + if (test__fixunstfsi(-0x1.23456789abcdefp+3)) + return 1; +#else + printf("skipped\n"); + +#endif + return 0; +}