From 06696f40afe58a231e2531c8bf6d9f0dadb92e51 Mon Sep 17 00:00:00 2001
From: Jakob Kaivo <jkk@ung.org>
Date: Sun, 3 Mar 2019 21:25:50 -0500
Subject: outline details from C18 annex F

---
 src/math/log1p.c | 17 +++++++++++++++++
 1 file changed, 17 insertions(+)

(limited to 'src/math/log1p.c')

diff --git a/src/math/log1p.c b/src/math/log1p.c
index 30e4a8d1..5f3e2684 100644
--- a/src/math/log1p.c
+++ b/src/math/log1p.c
@@ -1,9 +1,26 @@
 # define TGSOURCE "log1p.c"
 #include "_tgmath.h"
 #include <math.h>
+#include "fenv.h"
 
 TYPE TGFN(log1p)(TYPE x)
 {
+	switch (fpclassify(x)) {
+	case FP_ZERO:		return x;
+	case FP_INFINITE:	if (!signbit(x)) { return x; } break;
+	default:		break;
+	}
+
+	if (x == -1) {
+		feraiseexcept(FE_DIVBYZERO);
+		return - INFINITY;
+	}
+
+	if (x < -1) {
+		feraiseexcept(FE_INVALID);
+		return NAN;
+	}
+
 	return x;
 }
 
-- 
cgit v1.2.1