From fcbf63e62c627deae76c1b8cb8c0876c536ed811 Mon Sep 17 00:00:00 2001
From: Jari Vetoniemi <jari.vetoniemi@indooratlas.com>
Date: Mon, 16 Mar 2020 18:49:26 +0900
Subject: Fresh start

---
 jni/ruby/missing/nextafter.c | 77 ++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 77 insertions(+)
 create mode 100644 jni/ruby/missing/nextafter.c

(limited to 'jni/ruby/missing/nextafter.c')

diff --git a/jni/ruby/missing/nextafter.c b/jni/ruby/missing/nextafter.c
new file mode 100644
index 0000000..dd1f1f2
--- /dev/null
+++ b/jni/ruby/missing/nextafter.c
@@ -0,0 +1,77 @@
+#include "ruby/missing.h"
+
+#include <math.h>
+#include <float.h>
+
+/* This function doesn't set errno.  It should on POSIX, though. */
+
+double
+nextafter(double x, double y)
+{
+    double x1, x2, d;
+    int e;
+
+    if (isnan(x))
+        return x;
+    if (isnan(y))
+        return y;
+
+    if (x == y)
+        return y;
+
+    if (x == 0) {
+        /* the minimum "subnormal" float */
+        x1 = ldexp(0.5, DBL_MIN_EXP - DBL_MANT_DIG + 1);
+        if (x1 == 0)
+            x1 = DBL_MIN; /* the minimum "normal" float */
+        if (0 < y)
+            return x1;
+        else
+            return -x1;
+    }
+
+    if (x < 0) {
+        if (isinf(x))
+            return -DBL_MAX;
+        if (x == -DBL_MAX && y < 0 && isinf(y))
+            return y;
+    }
+    else {
+        if (isinf(x))
+            return DBL_MAX;
+        if (x == DBL_MAX && 0 < y && isinf(y))
+            return y;
+    }
+
+    x1 = frexp(x, &e);
+
+    if (x < y) {
+        d = DBL_EPSILON/2;
+        if (x1 == -0.5) {
+            x1 *= 2;
+            e--;
+        }
+    }
+    else {
+        d = -DBL_EPSILON/2;
+        if (x1 == 0.5) {
+            x1 *= 2;
+            e--;
+        }
+    }
+
+    if (e < DBL_MIN_EXP) {
+        d = ldexp(d, DBL_MIN_EXP-e);
+    }
+
+    x2 = x1 + d;
+
+    if (x2 == 0.0) {
+        if (x1 < 0)
+            return -0.0;
+        else
+            return +0.0;
+    }
+
+    return ldexp(x2, e);
+}
-- 
cgit v1.2.3