crt: arm: Add implementations of remainder and remquo for msvcrt Signed-off-by: Martin Storsjö <martin@martin.st>

diff --git a/mingw-w64-crt/Makefile.am b/mingw-w64-crt/Makefile.am
index 3fea74e..64d627e 100644
--- a/mingw-w64-crt/Makefile.am
+++ b/mingw-w64-crt/Makefile.am

@@ -259,6 +259,12 @@
   math/arm/truncf.S \
   math/arm-common/copysignl.c \
   math/arm-common/log2.c \
+  math/arm-common/remainder.c \
+  math/arm-common/remainderf.c \
+  math/arm-common/remainderl.c \
+  math/arm-common/remquo.c \
+  math/arm-common/remquof.c \
+  math/arm-common/remquol.c \
   math/arm-common/scalbn.c
 endif
 
@@ -266,6 +272,12 @@
   $(src_msvcrt) \
   math/arm-common/copysignl.c \
   math/arm-common/log2.c \
+  math/arm-common/remainder.c \
+  math/arm-common/remainderf.c \
+  math/arm-common/remainderl.c \
+  math/arm-common/remquo.c \
+  math/arm-common/remquof.c \
+  math/arm-common/remquol.c \
   math/arm-common/scalbn.c \
   math/arm64/exp2.S \
   math/arm64/exp2f.S \

diff --git a/mingw-w64-crt/math/arm-common/remainder.c b/mingw-w64-crt/math/arm-common/remainder.c
new file mode 100644
index 0000000..b97ccf1
--- /dev/null
+++ b/mingw-w64-crt/math/arm-common/remainder.c

@@ -0,0 +1,14 @@
+/**
+ * This file has no copyright assigned and is placed in the Public Domain.
+ * This file is part of the mingw-w64 runtime package.
+ * No warranty is given; refer to the file DISCLAIMER.PD within this package.
+ */
+
+#include <math.h>
+#include <errno.h>
+
+double remainder(double x, double y)
+{
+  int iret;
+  return remquo(x, y, &iret);
+}

diff --git a/mingw-w64-crt/math/arm-common/remainderf.c b/mingw-w64-crt/math/arm-common/remainderf.c
new file mode 100644
index 0000000..1e2946c
--- /dev/null
+++ b/mingw-w64-crt/math/arm-common/remainderf.c

@@ -0,0 +1,14 @@
+/**
+ * This file has no copyright assigned and is placed in the Public Domain.
+ * This file is part of the mingw-w64 runtime package.
+ * No warranty is given; refer to the file DISCLAIMER.PD within this package.
+ */
+
+#include <math.h>
+#include <errno.h>
+
+float remainderf(float x, float y)
+{
+  int iret;
+  return remquof(x, y, &iret);
+}

diff --git a/mingw-w64-crt/math/arm-common/remainderl.c b/mingw-w64-crt/math/arm-common/remainderl.c
new file mode 100644
index 0000000..16fef94
--- /dev/null
+++ b/mingw-w64-crt/math/arm-common/remainderl.c

@@ -0,0 +1,17 @@
+/**
+ * This file has no copyright assigned and is placed in the Public Domain.
+ * This file is part of the mingw-w64 runtime package.
+ * No warranty is given; refer to the file DISCLAIMER.PD within this package.
+ */
+
+#include <math.h>
+#include <errno.h>
+
+long double remainderl(long double x, long double y)
+{
+#if defined(__arm__) || defined(_ARM_) || defined(__aarch64__) || defined(_ARM64_)
+  return remainder(x, y);
+#else
+#error Not supported on your platform yet
+#endif
+}

diff --git a/mingw-w64-crt/math/arm-common/remquo.c b/mingw-w64-crt/math/arm-common/remquo.c
new file mode 100644
index 0000000..6a97ade
--- /dev/null
+++ b/mingw-w64-crt/math/arm-common/remquo.c

@@ -0,0 +1,31 @@
+/**
+ * This file has no copyright assigned and is placed in the Public Domain.
+ * This file is part of the mingw-w64 runtime package.
+ * No warranty is given; refer to the file DISCLAIMER.PD within this package.
+ */
+
+#include <math.h>
+#include <errno.h>
+
+double remquo(double x, double y, int *quo)
+{
+  if (isinf(x) || y == 0)
+    return NAN;
+  double div = x/y;
+  double integral;
+  double frac = modf(div, &integral);
+  int iintegral = (int)integral;
+  if (frac == 0.5) {
+    if (iintegral & 1)
+      *quo = iintegral + 1;
+    else
+      *quo = iintegral;
+  } else if (frac == -0.5) {
+    if (iintegral & 1)
+      *quo = iintegral - 1;
+    else
+      *quo = iintegral;
+  } else
+    *quo = round(div);
+  return x - *quo * y;
+}

diff --git a/mingw-w64-crt/math/arm-common/remquof.c b/mingw-w64-crt/math/arm-common/remquof.c
new file mode 100644
index 0000000..c395b2b
--- /dev/null
+++ b/mingw-w64-crt/math/arm-common/remquof.c

@@ -0,0 +1,31 @@
+/**
+ * This file has no copyright assigned and is placed in the Public Domain.
+ * This file is part of the mingw-w64 runtime package.
+ * No warranty is given; refer to the file DISCLAIMER.PD within this package.
+ */
+
+#include <math.h>
+#include <errno.h>
+
+float remquof(float x, float y, int *quo)
+{
+  if (isinf(x) || y == 0)
+    return NAN;
+  float div = x/y;
+  float integral;
+  float frac = modff(div, &integral);
+  int iintegral = (int)integral;
+  if (frac == 0.5) {
+    if (iintegral & 1)
+      *quo = iintegral + 1;
+    else
+      *quo = iintegral;
+  } else if (frac == -0.5) {
+    if (iintegral & 1)
+      *quo = iintegral - 1;
+    else
+      *quo = iintegral;
+  } else
+    *quo = roundf(div);
+  return x - *quo * y;
+}

diff --git a/mingw-w64-crt/math/arm-common/remquol.c b/mingw-w64-crt/math/arm-common/remquol.c
new file mode 100644
index 0000000..c9fb1cb
--- /dev/null
+++ b/mingw-w64-crt/math/arm-common/remquol.c

@@ -0,0 +1,17 @@
+/**
+ * This file has no copyright assigned and is placed in the Public Domain.
+ * This file is part of the mingw-w64 runtime package.
+ * No warranty is given; refer to the file DISCLAIMER.PD within this package.
+ */
+
+#include <math.h>
+#include <errno.h>
+
+long double remquol(long double x, long double y, int *quo)
+{
+#if defined(__arm__) || defined(_ARM_) || defined(__aarch64__) || defined(_ARM64_)
+    return remquo(x, y, quo);
+#else
+#error Not supported on your platform yet
+#endif
+}