1 files changed, 33 insertions, 0 deletions
diff --git a/contrib/libf2c/libF77/d_mod.c b/contrib/libf2c/libF77/d_mod.c
new file mode 100644
index 0000000..15cedef
--- /dev/null
+++ b/contrib/libf2c/libF77/d_mod.c
@@ -0,0 +1,33 @@
+#include "f2c.h"
+
+#ifdef IEEE_drem
+double drem (double, double);
+#else
+#undef abs
+#include <math.h>
+#endif
+double
+d_mod (doublereal * x, doublereal * y)
+{
+#ifdef IEEE_drem
+  double xa, ya, z;
+  if ((ya = *y) < 0.)
+    ya = -ya;
+  z = drem (xa = *x, ya);
+  if (xa > 0)
+    {
+      if (z < 0)
+	z += ya;
+    }
+  else if (z > 0)
+    z -= ya;
+  return z;
+#else
+  double quotient;
+  if ((quotient = *x / *y) >= 0)
+    quotient = floor (quotient);
+  else
+    quotient = -floor (-quotient);
+  return (*x - (*y) * quotient);
+#endif
+}