X-Git-Url: http://nsz.repo.hu/git/?a=blobdiff_plain;f=src%2Fmath%2Fexp10f.c;h=d009f0a84fdd9f15010ee83de729b86e50f69caf;hb=7590203c486d9002522019045d34ee3dee0a66f5;hp=c9521411ec5da7629596b7542ac31bfeda169572;hpb=f6819755779a084bf2f82cb90175a4d9a018de73;p=musl

diff --git a/src/math/exp10f.c b/src/math/exp10f.c
index c9521411..d009f0a8 100644
--- a/src/math/exp10f.c
+++ b/src/math/exp10f.c
@@ -1,17 +1,22 @@
 #define _GNU_SOURCE
 #include <math.h>
+#include <stdint.h>
 
 float exp10f(float x)
 {
 	static const float p10[] = {
-		1e-7, 1e-6, 1e-5, 1e-4, 1e-3, 1e-2, 1e-1,
+		1e-7f, 1e-6f, 1e-5f, 1e-4f, 1e-3f, 1e-2f, 1e-1f,
 		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];
 	}
 	return exp2(3.32192809488736234787031942948939 * x);
 }
+
+weak_alias(exp10f, pow10f);