1 files changed, 38 insertions, 0 deletions
diff --git a/libm/s_remquo.c b/libm/s_remquo.c
new file mode 100644
index 000000000..18853e211
--- /dev/null
+++ b/libm/s_remquo.c
@@ -0,0 +1,38 @@
+/* Copyright (C) 2002 by  Red Hat, Incorporated. All rights reserved.
+ *
+ * Permission to use, copy, modify, and distribute this software
+ * is freely granted, provided that this notice is preserved.
+ */
+
+#include "math.h"
+#include "math_private.h"
+
+libm_hidden_proto(remquo)
+#ifdef __STDC__
+	double remquo(double x, double y, int *quo)	/* wrapper remquo */
+#else
+	double remquo(x,y,quo)			/* wrapper remquo */
+	double x,y;
+        int *quo;
+#endif
+{
+        int signx, signy, signres;
+        int mswx;
+        int mswy;
+        double x_over_y;
+
+        GET_HIGH_WORD(mswx, x);
+        GET_HIGH_WORD(mswy, y);
+
+        signx = (mswx & 0x80000000) >> 31;
+        signy = (mswy & 0x80000000) >> 31;
+
+        signres = (signx ^ signy) ? -1 : 1;
+
+        x_over_y = fabs(x / y);
+
+        *quo = signres * (lrint(x_over_y) & 0x7f);
+
+        return remainder(x,y);
+}
+libm_hidden_def(remquo)