diff options
| author | Thomas Schwinge <thomas@codesourcery.com> | 2013-04-02 13:51:02 +0200 |
|---|---|---|
| committer | Thomas Schwinge <thomas@codesourcery.com> | 2013-04-02 13:51:02 +0200 |
| commit | 572676160d5639edc0ecb663147bd291841458d1 (patch) | |
| tree | 26abea75b21e81f568075075249aa3dbedad20c7 | |
| parent | 60c414c346a1d5ef0510ffbdc0ab75f288ee4d3f (diff) | |
| download | glibc-572676160d5639edc0ecb663147bd291841458d1.tar.xz glibc-572676160d5639edc0ecb663147bd291841458d1.zip | |
New <math.h> macro named issignaling to check for a signaling NaN (sNaN).
It is based on draft TS 18661 and currently enabled as a GNU extension.
49 files changed, 692 insertions, 16 deletions
@@ -1,3 +1,35 @@ +2013-04-02 Thomas Schwinge <thomas@codesourcery.com> + + * Versions.def (libm): Add GLIBC_2.18. + * include/math.h [!_ISOMAC] (__issignaling, __issignalingf): Add + hidden libm prototypes. + [!_ISOMAC && !__NO_LONG_DOUBLE_MATH] (__issignalingl): Likewise. + * math/Makefile (libm-calls): Add s_issignaling. + * math/Versions (libm: GLIBC_2.18): Add __issignaling, + __issignalingf, __issignalingl. Adjust all libm.abilist files. + * math/bits/mathcalls.h [__USE_GNU] (__issignaling): New prototype + declaration. + * math/math.h [__USE_GNU] (issignaling): New macro. + * sysdeps/ieee754/flt-32/s_issignalingf.c: New file. + * sysdeps/ieee754/dbl-64/s_issignaling.c: Likewise. + * sysdeps/ieee754/dbl-64/wordsize-64/s_issignaling.c: Likewise. + * sysdeps/ieee754/ldbl-96/s_issignalingl.c: Likewise. + * sysdeps/ieee754/ldbl-128/s_issignalingl.c: Likewise. + * sysdeps/ieee754/ldbl-128ibm/s_issignalingl.c: Likewise. + * manual/arith.texi (issignaling): New section. + * manual/libm-err-tab.pl (@all_functions): Update comment. + * math/gen-libm-test.pl (parse_args): Apply special handling for + issignaling. + * math/libm-test.inc (print_float, issignaling_test): New + functions. + (check_float_internal): Add issignaling checks. + * sysdeps/generic/math-tests.h (SNAN_TESTS_TYPE_CAST): Provide + default definition. + * sysdeps/powerpc/math-tests.h: New file. + * math/basic-test.c (TEST_FUNC, TEST_TRUNC): Add issignaling + tests. + * math/test-snan.c (TEST_FUNC): Likewise. + 2013-03-30 David S. Miller <davem@davemloft.net> * po/de.po: Update from translation team. @@ -26,6 +26,8 @@ Version 2.18 * Added a benchmark framework to track performance of functions in glibc. +* New <math.h> macro named issignaling to check for a signaling NaN (sNaN). + It is based on draft TS 18661 and currently enabled as a GNU extension. Version 2.17 diff --git a/Versions.def b/Versions.def index 7c7d1f89d3..04622397fd 100644 --- a/Versions.def +++ b/Versions.def @@ -58,6 +58,7 @@ libm { GLIBC_2.3.4 GLIBC_2.4 GLIBC_2.15 + GLIBC_2.18 } libnsl { GLIBC_2.0 diff --git a/include/math.h b/include/math.h index 2b040c67fd..cbbce59501 100644 --- a/include/math.h +++ b/include/math.h @@ -23,11 +23,14 @@ hidden_proto (__isnanl) libm_hidden_proto (__fpclassify) libm_hidden_proto (__fpclassifyf) +libm_hidden_proto (__issignaling) +libm_hidden_proto (__issignalingf) libm_hidden_proto (__exp) libm_hidden_proto (__expf) # ifndef __NO_LONG_DOUBLE_MATH libm_hidden_proto (__fpclassifyl) +libm_hidden_proto (__issignalingl) libm_hidden_proto (__expl) libm_hidden_proto (__expm1l) # endif diff --git a/manual/arith.texi b/manual/arith.texi index d060ff9463..ce8844ec2d 100644 --- a/manual/arith.texi +++ b/manual/arith.texi @@ -387,6 +387,14 @@ to @end smallexample @end deftypefn +@comment math.h +@comment GNU +@deftypefn {Macro} int issignaling (@emph{float-type} @var{x}) +This macro returns a nonzero value if @var{x} is a signaling NaN +(sNaN). It is based on draft TS 18661 and currently enabled as a GNU +extension. +@end deftypefn + Another set of floating-point classification functions was provided by BSD. @Theglibc{} also supports these functions; however, we recommend that you use the ISO C99 macros in new code. Those are standard diff --git a/manual/libm-err-tab.pl b/manual/libm-err-tab.pl index 5aa26039ac..15646693fe 100755 --- a/manual/libm-err-tab.pl +++ b/manual/libm-err-tab.pl @@ -77,7 +77,9 @@ use vars qw (%results @all_floats %suffices @all_functions); "remainder", "remquo", "rint", "round", "scalb", "scalbn", "scalbln", "sin", "sincos", "sinh", "sqrt", "tan", "tanh", "tgamma", "trunc", "y0", "y1", "yn" ); -# "fpclassify", "isfinite", "isnormal", "signbit" are not tabulated +# fpclassify, isnormal, isfinite, isinf, isnan, issignaling, signbit, +# isgreater, isgreaterequal, isless, islessequal, islessgreater, isunordered +# are not tabulated. if ($#ARGV == 0) { $sources = $ARGV[0]; diff --git a/math/Makefile b/math/Makefile index e216dfb96a..f396ba206a 100644 --- a/math/Makefile +++ b/math/Makefile @@ -58,7 +58,7 @@ libm-calls = e_acos e_acosh e_asin e_atan2 e_atanh e_cosh e_exp e_fmod \ s_catan s_casin s_ccos s_csin s_ctan s_ctanh s_cacos \ s_casinh s_cacosh s_catanh s_csqrt s_cpow s_cproj s_clog10 \ s_fma s_lrint s_llrint s_lround s_llround e_exp10 w_log2 \ - s_isinf_ns $(calls:s_%=m_%) x2y2m1 k_casinh + s_isinf_ns s_issignaling $(calls:s_%=m_%) x2y2m1 k_casinh include ../Makeconfig diff --git a/math/Versions b/math/Versions index 09888510b1..513ab14244 100644 --- a/math/Versions +++ b/math/Versions @@ -198,4 +198,7 @@ libm { __gamma_r_finite; __gammaf_r_finite; __gammal_r_finite; __exp_finite; __expf_finite; __expl_finite; } + GLIBC_2.18 { + __issignaling; __issignalingf; __issignalingl; + } } diff --git a/math/basic-test.c b/math/basic-test.c index ffead2ec0f..9e9b848ee5 100644 --- a/math/basic-test.c +++ b/math/basic-test.c @@ -20,6 +20,9 @@ #include <float.h> #include <stdio.h> +#include <math-tests.h> + + static int errors = 0; @@ -39,6 +42,10 @@ NAME (void) \ /* Variables are declared volatile to forbid some compiler \ optimizations. */ \ volatile FLOAT Inf_var, qNaN_var, zero_var, one_var; \ + /* A sNaN is only guaranteed to be representable in variables with */ \ + /* static (or thread-local) storage duration. */ \ + static volatile FLOAT sNaN_var = __builtin_nans ## SUFFIX (""); \ + static volatile FLOAT minus_sNaN_var = -__builtin_nans ## SUFFIX (""); \ FLOAT x1, x2; \ \ zero_var = 0.0; \ @@ -49,6 +56,8 @@ NAME (void) \ (void) &zero_var; \ (void) &one_var; \ (void) &qNaN_var; \ + (void) &sNaN_var; \ + (void) &minus_sNaN_var; \ (void) &Inf_var; \ \ \ @@ -56,16 +65,41 @@ NAME (void) \ check (#FLOAT " isinf (-inf) == -1", isinf (-Inf_var) == -1); \ check (#FLOAT " !isinf (1)", !(isinf (one_var))); \ check (#FLOAT " !isinf (qNaN)", !(isinf (qNaN_var))); \ + if (SNAN_TESTS (FLOAT)) \ + check (#FLOAT " !isinf (sNaN)", !(isinf (sNaN_var))); \ \ check (#FLOAT " isnan (qNaN)", isnan (qNaN_var)); \ + if (SNAN_TESTS (FLOAT)) \ + check (#FLOAT " isnan (sNaN)", isnan (sNaN_var)); \ check (#FLOAT " isnan (-qNaN)", isnan (-qNaN_var)); \ + if (SNAN_TESTS (FLOAT)) \ + check (#FLOAT " isnan (-sNaN)", isnan (minus_sNaN_var)); \ check (#FLOAT " !isnan (1)", !(isnan (one_var))); \ check (#FLOAT " !isnan (inf)", !(isnan (Inf_var))); \ \ + check (#FLOAT " !issignaling (qNaN)", !(issignaling (qNaN_var))); \ + if (SNAN_TESTS (FLOAT)) \ + check (#FLOAT " issignaling (sNaN)", issignaling (sNaN_var)); \ + check (#FLOAT " !issignaling (-qNaN)", !(issignaling (-qNaN_var))); \ + if (SNAN_TESTS (FLOAT)) \ + check (#FLOAT " issignaling (-sNaN)", issignaling (minus_sNaN_var)); \ + check (#FLOAT " !issignaling (1)", !(issignaling (one_var))); \ + check (#FLOAT " !issignaling (inf)", !(issignaling (Inf_var))); \ + \ check (#FLOAT " inf == inf", Inf_var == Inf_var); \ check (#FLOAT " -inf == -inf", -Inf_var == -Inf_var); \ check (#FLOAT " inf != -inf", Inf_var != -Inf_var); \ check (#FLOAT " qNaN != qNaN", qNaN_var != qNaN_var); \ + if (SNAN_TESTS (FLOAT)) \ + check (#FLOAT " sNaN != sNaN", sNaN_var != sNaN_var); \ + if (SNAN_TESTS (FLOAT)) \ + check (#FLOAT " qNaN != sNaN", qNaN_var != sNaN_var); \ + if (SNAN_TESTS (FLOAT)) \ + check (#FLOAT " -sNaN != -sNaN", minus_sNaN_var != minus_sNaN_var); \ + if (SNAN_TESTS (FLOAT)) \ + check (#FLOAT " sNaN != -sNaN", sNaN_var != minus_sNaN_var); \ + if (SNAN_TESTS (FLOAT)) \ + check (#FLOAT " qNaN != -sNaN", qNaN_var != minus_sNaN_var); \ \ /* \ the same tests but this time with NAN from <bits/nan.h> \ @@ -76,6 +110,11 @@ NAME (void) \ check (#FLOAT " !isinf (NAN)", !(isinf (NAN))); \ check (#FLOAT " !isinf (-NAN)", !(isinf (-NAN))); \ check (#FLOAT " NAN != NAN", NAN != NAN); \ + check (#FLOAT " NAN != qNaN", NAN != qNaN_var); \ + if (SNAN_TESTS (FLOAT)) \ + check (#FLOAT " NAN != sNaN", NAN != sNaN_var); \ + if (SNAN_TESTS (FLOAT)) \ + check (#FLOAT " NAN != -sNaN", NAN != minus_sNaN_var); \ \ /* \ And again with the value returned by the `nan' function. \ @@ -86,6 +125,12 @@ NAME (void) \ check (#FLOAT " !isinf (-nan (\"\"))", !(isinf (-nan ## SUFFIX ("")))); \ check (#FLOAT " nan (\"\") != nan (\"\")", \ nan ## SUFFIX ("") != nan ## SUFFIX ("")); \ + check (#FLOAT " nan (\"\") != qNaN", nan ## SUFFIX ("") != qNaN_var); \ + if (SNAN_TESTS (FLOAT)) \ + check (#FLOAT " nan (\"\") != sNaN", nan ## SUFFIX ("") != sNaN_var); \ + if (SNAN_TESTS (FLOAT)) \ + check (#FLOAT " nan (\"\") != -sNaN", \ + nan ## SUFFIX ("") != minus_sNaN_var); \ \ /* test if EPSILON is ok */ \ x1 = 1.0; \ @@ -108,6 +153,9 @@ void \ NAME (void) \ { \ volatile DOUBLE Inf_var, qNaN_var, zero_var, one_var; \ + /* A sNaN is only guaranteed to be representable in variables with */ \ + /* static (or thread-local) storage duration. */ \ + static volatile DOUBLE sNaN_var = __builtin_nans ## SUFFIX (""); \ FLOAT x1, x2; \ \ zero_var = 0.0; \ @@ -116,10 +164,25 @@ NAME (void) \ Inf_var = one_var / zero_var; \ \ (void) &qNaN_var; \ + (void) &sNaN_var; \ (void) &Inf_var; \ \ x1 = (FLOAT) qNaN_var; \ - check (" "#FLOAT" x = ("#FLOAT") ("#DOUBLE") qNaN", isnan (x1) != 0); \ + check (" "#FLOAT" x = ("#FLOAT") ("#DOUBLE") qNaN, isnan", isnan (x1)); \ + check (" "#FLOAT" x = ("#FLOAT") ("#DOUBLE") qNaN, !issignaling", \ + !issignaling (x1)); \ + if (SNAN_TESTS (FLOAT)) \ + { \ + x1 = (FLOAT) sNaN_var; \ + check (" "#FLOAT" x = ("#FLOAT") ("#DOUBLE") sNaN, isnan", isnan (x1)); \ + if (SNAN_TESTS_TYPE_CAST) |