xref: /netbsd-src/sys/external/bsd/compiler_rt/dist/lib/builtins/floatdisf.c (revision d31434597b255f3fd7646698e04170d148f1de49)
1156cd587Sjoerg /*===-- floatdisf.c - Implement __floatdisf -------------------------------===
2156cd587Sjoerg  *
3156cd587Sjoerg  *                     The LLVM Compiler Infrastructure
4156cd587Sjoerg  *
5156cd587Sjoerg  * This file is dual licensed under the MIT and the University of Illinois Open
6156cd587Sjoerg  * Source Licenses. See LICENSE.TXT for details.
7156cd587Sjoerg  *
8156cd587Sjoerg  *===----------------------------------------------------------------------===
9156cd587Sjoerg  *
10156cd587Sjoerg  * This file implements __floatdisf for the compiler_rt library.
11156cd587Sjoerg  *
12156cd587Sjoerg  *===----------------------------------------------------------------------===
13156cd587Sjoerg  */
14156cd587Sjoerg 
15156cd587Sjoerg /* Returns: convert a to a float, rounding toward even.*/
16156cd587Sjoerg 
17156cd587Sjoerg /* Assumption: float is a IEEE 32 bit floating point type
18156cd587Sjoerg  *             di_int is a 64 bit integral type
19156cd587Sjoerg  */
20156cd587Sjoerg 
21156cd587Sjoerg /* seee eeee emmm mmmm mmmm mmmm mmmm mmmm */
22156cd587Sjoerg 
23156cd587Sjoerg #include "int_lib.h"
24156cd587Sjoerg 
25156cd587Sjoerg COMPILER_RT_ABI float
__floatdisf(di_int a)26156cd587Sjoerg __floatdisf(di_int a)
27156cd587Sjoerg {
28156cd587Sjoerg     if (a == 0)
29156cd587Sjoerg         return 0.0F;
30156cd587Sjoerg     const unsigned N = sizeof(di_int) * CHAR_BIT;
31156cd587Sjoerg     const di_int s = a >> (N-1);
32156cd587Sjoerg     a = (a ^ s) - s;
33156cd587Sjoerg     int sd = N - __builtin_clzll(a);  /* number of significant digits */
34156cd587Sjoerg     int e = sd - 1;             /* exponent */
35156cd587Sjoerg     if (sd > FLT_MANT_DIG)
36156cd587Sjoerg     {
37156cd587Sjoerg         /*  start:  0000000000000000000001xxxxxxxxxxxxxxxxxxxxxxPQxxxxxxxxxxxxxxxxxx
38156cd587Sjoerg          *  finish: 000000000000000000000000000000000000001xxxxxxxxxxxxxxxxxxxxxxPQR
39156cd587Sjoerg          *                                                12345678901234567890123456
40156cd587Sjoerg          *  1 = msb 1 bit
41156cd587Sjoerg          *  P = bit FLT_MANT_DIG-1 bits to the right of 1
42156cd587Sjoerg          *  Q = bit FLT_MANT_DIG bits to the right of 1
43156cd587Sjoerg          *  R = "or" of all bits to the right of Q
44156cd587Sjoerg          */
45156cd587Sjoerg         switch (sd)
46156cd587Sjoerg         {
47156cd587Sjoerg         case FLT_MANT_DIG + 1:
48156cd587Sjoerg             a <<= 1;
49156cd587Sjoerg             break;
50156cd587Sjoerg         case FLT_MANT_DIG + 2:
51156cd587Sjoerg             break;
52156cd587Sjoerg         default:
53156cd587Sjoerg             a = ((du_int)a >> (sd - (FLT_MANT_DIG+2))) |
54156cd587Sjoerg                 ((a & ((du_int)(-1) >> ((N + FLT_MANT_DIG+2) - sd))) != 0);
55156cd587Sjoerg         };
56156cd587Sjoerg         /* finish: */
57156cd587Sjoerg         a |= (a & 4) != 0;  /* Or P into R */
58156cd587Sjoerg         ++a;  /* round - this step may add a significant bit */
59156cd587Sjoerg         a >>= 2;  /* dump Q and R */
60156cd587Sjoerg         /* a is now rounded to FLT_MANT_DIG or FLT_MANT_DIG+1 bits */
61156cd587Sjoerg         if (a & ((du_int)1 << FLT_MANT_DIG))
62156cd587Sjoerg         {
63156cd587Sjoerg             a >>= 1;
64156cd587Sjoerg             ++e;
65156cd587Sjoerg         }
66156cd587Sjoerg         /* a is now rounded to FLT_MANT_DIG bits */
67156cd587Sjoerg     }
68156cd587Sjoerg     else
69156cd587Sjoerg     {
70156cd587Sjoerg         a <<= (FLT_MANT_DIG - sd);
71156cd587Sjoerg         /* a is now rounded to FLT_MANT_DIG bits */
72156cd587Sjoerg     }
73156cd587Sjoerg     float_bits fb;
74156cd587Sjoerg     fb.u = ((su_int)s & 0x80000000) |  /* sign */
75156cd587Sjoerg            ((e + 127) << 23)       |  /* exponent */
76156cd587Sjoerg            ((su_int)a & 0x007FFFFF);   /* mantissa */
77156cd587Sjoerg     return fb.f;
78156cd587Sjoerg }
793044ee7eSrin 
803044ee7eSrin #if defined(__ARM_EABI__)
81*d3143459Srin #if defined(COMPILER_RT_ARMHF_TARGET)
__aeabi_l2f(di_int a)823044ee7eSrin AEABI_RTABI float __aeabi_l2f(di_int a) {
833044ee7eSrin   return __floatdisf(a);
843044ee7eSrin }
85*d3143459Srin #else
86*d3143459Srin AEABI_RTABI float __aeabi_l2f(di_int a) COMPILER_RT_ALIAS(__floatdisf);
873044ee7eSrin #endif
88*d3143459Srin #endif
89