fix pow(+-0, inf) case http://www.open-std.org/jtc1/sc22/wg14/www/docs/n1515.pdf

[libc-test] / src / math / special / powf.h

diff --git a/src/math/special/powf.h b/src/math/special/powf.h
index bec842d..8c7d72d 100644 (file)
--- a/src/math/special/powf.h
+++ b/src/math/special/powf.h
@@ -11,7 +11,7 @@ T(RN,          0x0p+0,         -0x1p+0,             inf,          0x0p+0, DIVBYZ
 T(RN,          0x0p+0,         -0x1p+1,             inf,          0x0p+0, DIVBYZERO)
 T(RN,          0x0p+0,       -0x1.8p+1,             inf,          0x0p+0, DIVBYZERO)
 T(RN,          0x0p+0,         -0x1p+2,             inf,          0x0p+0, DIVBYZERO)
-T(RN,          0x0p+0,            -inf,             inf,          0x0p+0, DIVBYZERO)
+T(RN,          0x0p+0,            -inf,             inf,          0x0p+0, 0)
 T(RN,         -0x0p+0,             nan,             nan,          0x0p+0, 0)
 T(RN,         -0x0p+0,             inf,          0x0p+0,          0x0p+0, 0)
 T(RN,         -0x0p+0,        0x1.8p+1,         -0x0p+0,          0x0p+0, 0)
@@ -25,7 +25,7 @@ T(RN,         -0x0p+0,         -0x1p+0,            -inf,          0x0p+0, DIVBYZ
 T(RN,         -0x0p+0,         -0x1p+1,             inf,          0x0p+0, DIVBYZERO)
 T(RN,         -0x0p+0,       -0x1.8p+1,            -inf,          0x0p+0, DIVBYZERO)
 T(RN,         -0x0p+0,         -0x1p+2,             inf,          0x0p+0, DIVBYZERO)
-T(RN,         -0x0p+0,            -inf,             inf,          0x0p+0, DIVBYZERO)
+T(RN,         -0x0p+0,            -inf,             inf,          0x0p+0, 0)
 T(RN,             nan,          0x0p+0,          0x1p+0,          0x0p+0, 0)
 T(RN,             inf,          0x0p+0,          0x1p+0,          0x0p+0, 0)
 T(RN,            -inf,          0x0p+0,          0x1p+0,          0x0p+0, 0)