Add more tests of pow.

Prompted by a gcc-patches discussion, this patch adds tests of pow for
the cases where pow (x, 0.5) is required to return a different result
from sqrt (x), as those cases were previously missing from the tests
(although they worked correctly).

Tested for x86_64 and x86.

	* math/auto-libm-test-in: Add another test of pow.
	* math/auto-libm-test-out: Regenerated.
	* math/libm-test.inc (pow_test_data): Add another test.

diff --git a/ChangeLog b/ChangeLog
index 833cdcc01f0f556e307089f11e4e0ae0274a0a44..8a0882c661313962b58fc1a1711524e308334a67 100644 (file)
--- a/ChangeLog
+++ b/ChangeLog
@@ -1,3 +1,9 @@
+2015-11-10  Joseph Myers  <joseph@codesourcery.com>
+
+       * math/auto-libm-test-in: Add another test of pow.
+       * math/auto-libm-test-out: Regenerated.
+       * math/libm-test.inc (pow_test_data): Add another test.
+
 2015-11-10  Adhemerval Zanella  <adhemerval.zanella@linaro.org>
 
        [BZ #19219]

diff --git a/math/auto-libm-test-in b/math/auto-libm-test-in
index 0baca59e88413a470e30e0f99dbf7713ee8d9974..f1abcba578df51eb9a3fa6e478dd448fbcac02e2 100644 (file)
--- a/math/auto-libm-test-in
+++ b/math/auto-libm-test-in
@@ -3254,6 +3254,7 @@ pow 0.0 max
 pow 0.0 min_subnorm
 
 # pow (-0, y) == +0 for y > 0 and not an odd integer.
+pow -0 0.5
 pow -0 4
 pow -0 0x1p24
 pow -0 0x1p127

diff --git a/math/auto-libm-test-out b/math/auto-libm-test-out
index 22c4b4188130f381a3ef7a2fc8c2a004ef1a88cf..9d3819ca0125fa4969e156824b97e0d93027cabc 100644 (file)
--- a/math/auto-libm-test-out
+++ b/math/auto-libm-test-out
@@ -219802,6 +219802,31 @@ pow 0.0 min_subnorm
 = pow tonearest ldbl-128 0x0p+0L 0x4p-16496L : 0x0p+0L : inexact-ok
 = pow towardzero ldbl-128 0x0p+0L 0x4p-16496L : 0x0p+0L : inexact-ok
 = pow upward ldbl-128 0x0p+0L 0x4p-16496L : 0x0p+0L : inexact-ok
+pow -0 0.5
+= pow downward flt-32 -0x0p+0f 0x8p-4f : 0x0p+0f : inexact-ok
+= pow tonearest flt-32 -0x0p+0f 0x8p-4f : 0x0p+0f : inexact-ok
+= pow towardzero flt-32 -0x0p+0f 0x8p-4f : 0x0p+0f : inexact-ok
+= pow upward flt-32 -0x0p+0f 0x8p-4f : 0x0p+0f : inexact-ok
+= pow downward dbl-64 -0x0p+0 0x8p-4 : 0x0p+0 : inexact-ok
+= pow tonearest dbl-64 -0x0p+0 0x8p-4 : 0x0p+0 : inexact-ok
+= pow towardzero dbl-64 -0x0p+0 0x8p-4 : 0x0p+0 : inexact-ok
+= pow upward dbl-64 -0x0p+0 0x8p-4 : 0x0p+0 : inexact-ok
+= pow downward ldbl-96-intel -0x0p+0L 0x8p-4L : 0x0p+0L : inexact-ok
+= pow tonearest ldbl-96-intel -0x0p+0L 0x8p-4L : 0x0p+0L : inexact-ok
+= pow towardzero ldbl-96-intel -0x0p+0L 0x8p-4L : 0x0p+0L : inexact-ok
+= pow upward ldbl-96-intel -0x0p+0L 0x8p-4L : 0x0p+0L : inexact-ok
+= pow downward ldbl-96-m68k -0x0p+0L 0x8p-4L : 0x0p+0L : inexact-ok
+= pow tonearest ldbl-96-m68k -0x0p+0L 0x8p-4L : 0x0p+0L : inexact-ok
+= pow towardzero ldbl-96-m68k -0x0p+0L 0x8p-4L : 0x0p+0L : inexact-ok
+= pow upward ldbl-96-m68k -0x0p+0L 0x8p-4L : 0x0p+0L : inexact-ok
+= pow downward ldbl-128 -0x0p+0L 0x8p-4L : 0x0p+0L : inexact-ok
+= pow tonearest ldbl-128 -0x0p+0L 0x8p-4L : 0x0p+0L : inexact-ok
+= pow towardzero ldbl-128 -0x0p+0L 0x8p-4L : 0x0p+0L : inexact-ok
+= pow upward ldbl-128 -0x0p+0L 0x8p-4L : 0x0p+0L : inexact-ok
+= pow downward ldbl-128ibm -0x0p+0L 0x8p-4L : 0x0p+0L : inexact-ok
+= pow tonearest ldbl-128ibm -0x0p+0L 0x8p-4L : 0x0p+0L : inexact-ok
+= pow towardzero ldbl-128ibm -0x0p+0L 0x8p-4L : 0x0p+0L : inexact-ok
+= pow upward ldbl-128ibm -0x0p+0L 0x8p-4L : 0x0p+0L : inexact-ok
 pow -0 4
 = pow downward flt-32 -0x0p+0f 0x4p+0f : 0x0p+0f : inexact-ok
 = pow tonearest flt-32 -0x0p+0f 0x4p+0f : 0x0p+0f : inexact-ok

diff --git a/math/libm-test.inc b/math/libm-test.inc
index 669b060dc9f12f1f2652003839877c9c433c0681..a09cf90826828866ad22e8f014f0ad189f8e12de 100644 (file)
--- a/math/libm-test.inc
+++ b/math/libm-test.inc
@@ -10285,6 +10285,7 @@ static const struct test_ff_f_data pow_test_data[] =
 #endif
 
     /* pow (-inf, y) == +inf for y > 0 and not an odd integer.  */
+    TEST_ff_f (pow, minus_infty, 0.5, plus_infty, ERRNO_UNCHANGED),
     TEST_ff_f (pow, minus_infty, 28, plus_infty, ERRNO_UNCHANGED),
     TEST_ff_f (pow, minus_infty, 0x1p24, plus_infty, ERRNO_UNCHANGED),
     TEST_ff_f (pow, minus_infty, 0x1p127, plus_infty, ERRNO_UNCHANGED),