diff options
| author | Raghuram Subramani <raghus2247@gmail.com> | 2022-06-19 19:47:51 +0530 |
|---|---|---|
| committer | Raghuram Subramani <raghus2247@gmail.com> | 2022-06-19 19:47:51 +0530 |
| commit | 4fd287655a72b9aea14cdac715ad5b90ed082ed2 (patch) | |
| tree | 65d393bc0e699dd12d05b29ba568e04cea666207 /circuitpython/lib/libm/nearbyintf.c | |
| parent | 0150f70ce9c39e9e6dd878766c0620c85e47bed0 (diff) | |
add circuitpython code
Diffstat (limited to 'circuitpython/lib/libm/nearbyintf.c')
| -rw-r--r-- | circuitpython/lib/libm/nearbyintf.c | 24 |
1 files changed, 24 insertions, 0 deletions
diff --git a/circuitpython/lib/libm/nearbyintf.c b/circuitpython/lib/libm/nearbyintf.c new file mode 100644 index 0000000..5ac3178 --- /dev/null +++ b/circuitpython/lib/libm/nearbyintf.c @@ -0,0 +1,24 @@ +// adapted from the rintf() function from musl-1.1.16 + +#include "libm.h" + +float nearbyintf(float x) +{ + union {float f; uint32_t i;} u = {x}; + int e = u.i>>23 & 0xff; + int s = u.i>>31; + float_t y; + + if (e >= 0x7f+23) + return x; + if (s) + y = x - 0x1p23f + 0x1p23f; + else + y = x + 0x1p23f - 0x1p23f; +#pragma GCC diagnostic push +#pragma GCC diagnostic ignored "-Wfloat-equal" + if (y == 0) + return s ? -0.0f : 0.0f; +#pragma GCC diagnostic pop + return y; +} |
