aboutsummaryrefslogtreecommitdiff
path: root/circuitpython/lib/libm/nearbyintf.c
diff options
context:
space:
mode:
authorRaghuram Subramani <raghus2247@gmail.com>2022-06-19 19:47:51 +0530
committerRaghuram Subramani <raghus2247@gmail.com>2022-06-19 19:47:51 +0530
commit4fd287655a72b9aea14cdac715ad5b90ed082ed2 (patch)
tree65d393bc0e699dd12d05b29ba568e04cea666207 /circuitpython/lib/libm/nearbyintf.c
parent0150f70ce9c39e9e6dd878766c0620c85e47bed0 (diff)
add circuitpython code
Diffstat (limited to 'circuitpython/lib/libm/nearbyintf.c')
-rw-r--r--circuitpython/lib/libm/nearbyintf.c24
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;
+}