math: add fma TODO comments about the underflow issue
[oweals/musl.git] / src / math / modff.c
1 #include "libm.h"
2
3 float modff(float x, float *iptr)
4 {
5         union {float x; uint32_t n;} u = {x};
6         uint32_t mask;
7         int e;
8
9         e = (int)(u.n>>23 & 0xff) - 0x7f;
10
11         /* no fractional part */
12         if (e >= 23) {
13                 *iptr = x;
14                 if (e == 0x80 && u.n<<9 != 0) { /* nan */
15                         return x;
16                 }
17                 u.n &= 0x80000000;
18                 return u.x;
19         }
20         /* no integral part */
21         if (e < 0) {
22                 u.n &= 0x80000000;
23                 *iptr = u.x;
24                 return x;
25         }
26
27         mask = 0x007fffff>>e;
28         if ((u.n & mask) == 0) {
29                 *iptr = x;
30                 u.n &= 0x80000000;
31                 return u.x;
32         }
33         u.n &= ~mask;
34         *iptr = u.x;
35         STRICT_ASSIGN(float, x, x - *iptr);
36         return x;
37 }