outline details from C18 annex F

1 files changed, 17 insertions, 0 deletions

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;
 }