X-Git-Url: http://nsz.repo.hu/git/?a=blobdiff_plain;f=src%2Fmath%2Fexp10f.c;h=d009f0a84fdd9f15010ee83de729b86e50f69caf;hb=0c277ff156749628c678257f878d3a6edb5868de;hp=5fd1af9c7d8180ec9e5eea88a44be17d9c527741;hpb=a764db9a086ca036d4fc9181f96d19ab312a6560;p=musl

diff --git a/src/math/exp10f.c b/src/math/exp10f.c
index 5fd1af9c..d009f0a8 100644
--- a/src/math/exp10f.c
+++ b/src/math/exp10f.c
@@ -1,6 +1,6 @@
 #define _GNU_SOURCE
 #include <math.h>
-#include "libc.h"
+#include <stdint.h>
 
 float exp10f(float x)
 {
@@ -9,7 +9,9 @@ float exp10f(float x)
 		1, 1e1, 1e2, 1e3, 1e4, 1e5, 1e6, 1e7
 	};
 	float n, y = modff(x, &n);
-	if (fabsf(n) < 8) {
+	union {float f; uint32_t i;} u = {n};
+	/* fabsf(n) < 8 without raising invalid on nan */
+	if ((u.i>>23 & 0xff) < 0x7f+3) {
 		if (!y) return p10[(int)n+7];
 		y = exp2f(3.32192809488736234787031942948939f * y);
 		return y * p10[(int)n+7];