aboutsummaryrefslogtreecommitdiff
path: root/compiler-rt/lib/builtins/floattitf.c
diff options
context:
space:
mode:
Diffstat (limited to 'compiler-rt/lib/builtins/floattitf.c')
-rw-r--r--compiler-rt/lib/builtins/floattitf.c55
1 files changed, 6 insertions, 49 deletions
diff --git a/compiler-rt/lib/builtins/floattitf.c b/compiler-rt/lib/builtins/floattitf.c
index fff0755c3bb4..5dffe22fdb4e 100644
--- a/compiler-rt/lib/builtins/floattitf.c
+++ b/compiler-rt/lib/builtins/floattitf.c
@@ -16,6 +16,11 @@
#include "fp_lib.h"
#include "int_lib.h"
+#if defined(CRT_HAS_TF_MODE)
+#define SRC_I128
+#define DST_QUAD
+#include "int_to_fp_impl.inc"
+
// Returns: convert a ti_int to a fp_t, rounding toward even.
// Assumption: fp_t is a IEEE 128 bit floating point type
@@ -25,54 +30,6 @@
// mmmm | mmmm mmmm mmmm mmmm mmmm mmmm mmmm mmmm | mmmm mmmm mmmm mmmm mmmm
// mmmm mmmm mmmm
-#if defined(CRT_HAS_TF_MODE)
-COMPILER_RT_ABI fp_t __floattitf(ti_int a) {
- if (a == 0)
- return 0.0;
- const unsigned N = sizeof(ti_int) * CHAR_BIT;
- const ti_int s = a >> (N - 1);
- a = (a ^ s) - s;
- int sd = N - __clzti2(a); // number of significant digits
- int e = sd - 1; // exponent
- if (sd > TF_MANT_DIG) {
- // start: 0000000000000000000001xxxxxxxxxxxxxxxxxxxxxxPQxxxxxxxxxxxxxxxxxx
- // finish: 000000000000000000000000000000000000001xxxxxxxxxxxxxxxxxxxxxxPQR
- // 12345678901234567890123456
- // 1 = msb 1 bit
- // P = bit LDBL_MANT_DIG-1 bits to the right of 1
- // Q = bit LDBL_MANT_DIG bits to the right of 1
- // R = "or" of all bits to the right of Q
- switch (sd) {
- case TF_MANT_DIG + 1:
- a <<= 1;
- break;
- case TF_MANT_DIG + 2:
- break;
- default:
- a = ((tu_int)a >> (sd - (TF_MANT_DIG + 2))) |
- ((a & ((tu_int)(-1) >> ((N + TF_MANT_DIG + 2) - sd))) != 0);
- };
- // finish:
- a |= (a & 4) != 0; // Or P into R
- ++a; // round - this step may add a significant bit
- a >>= 2; // dump Q and R
- // a is now rounded to LDBL_MANT_DIG or LDBL_MANT_DIG+1 bits
- if (a & ((tu_int)1 << TF_MANT_DIG)) {
- a >>= 1;
- ++e;
- }
- // a is now rounded to LDBL_MANT_DIG bits
- } else {
- a <<= (TF_MANT_DIG - sd);
- // a is now rounded to LDBL_MANT_DIG bits
- }
-
- long_double_bits fb;
- fb.u.high.all = (s & 0x8000000000000000LL) // sign
- | (du_int)(e + 16383) << 48 // exponent
- | ((a >> 64) & 0x0000ffffffffffffLL); // significand
- fb.u.low.all = (du_int)(a);
- return fb.f;
-}
+COMPILER_RT_ABI fp_t __floattitf(ti_int a) { return __floatXiYf__(a); }
#endif