1 files changed, 33 insertions, 0 deletions
diff --git a/src/math/modf.c b/src/math/modf.c
new file mode 100644
index 00000000..73c5d3f3
--- /dev/null
+++ b/src/math/modf.c
@@ -0,0 +1,33 @@
+# define TGSOURCE "modf.c"
+#include <math.h>
+#include "nonstd/tgmath.h"
+#include "errno.h"
+
+/** decompose floating-point numbers **/
+TYPE TGFN(modf)(TYPE value, TYPE *iptr)
+{
+	if (0) {
+		errno = ERANGE; /* The result cannot be represented */
+		/* RETURN_FAILURE(CONSTANT(HUGE_VAL), A range error occurred); */
+		return TGHUGE;
+	}
+
+	/* RETURN_SUCCESS(the signed fractional part of ARGUMENT(value)); */
+	(void)iptr;
+	return value;
+}
+
+/***
+break ARGUMENT(value) into integral and fractional parts.
+Each part has the same sign as ARGUMENT(value). The integral part is converted to
+the same floating-point type as ARGUMENT(value) and stored in the object pointed to
+by ARGUMENT(iptr).
+***/
+
+/*
+IMPLEMENTATION(The value returned on a domain error, CONSTANT(HUGE_VAL))
+LINK(m)
+*/
+/*
+STDC(1)
+*/