Import sources from 2011-03-30 iso image

1 files changed, 31 insertions, 0 deletions

diff --git a/sys/src/libc/port/fmod.c b/sys/src/libc/port/fmod.c
new file mode 100755
index 000000000..3b97a5b43
--- /dev/null
+++ b/sys/src/libc/port/fmod.c
@@ -0,0 +1,31 @@
+#include <u.h>
+#include <libc.h>
+
+/*
+ * floating-point mod function without infinity or NaN checking
+ */
+double
+fmod (double x, double y)
+{
+	int sign, yexp, rexp;
+	double r, yfr, rfr;
+
+	if (y == 0)
+		return x;
+	if (y < 0)
+		y = -y;
+	yfr = frexp(y, &yexp);
+	sign = 0;
+	if(x < 0) {
+		r = -x;
+		sign++;
+	} else
+		r = x;
+	while(r >= y) {
+		rfr = frexp(r, &rexp);
+		r -= ldexp(y, rexp - yexp - (rfr < yfr));
+	}
+	if(sign)
+		r = -r;
+	return r;
+}