diff options
51 files changed, 9778 insertions, 1 deletions
@@ -31,6 +31,13 @@ Major new features: * The iconv program now supports converting files in place. The program automatically uses a temporary file if required. +* The following ISO C23 function families (introduced in TS + 18661-4:2015) are now supported in <math.h>. Each family includes + functions for float, double, long double, _FloatN and _FloatNx, and a + type-generic macro in <tgmath.h>. + + - Trigonometric functions: cospi. + * The GNU C Library now supports a feature test macro _ISOC2Y_SOURCE to enable features from the draft ISO C2Y standard. Only some features from this draft standard are supported by the GNU C Library, and as the draft diff --git a/manual/math.texi b/manual/math.texi index e4f9122609..b818db2c1b 100644 --- a/manual/math.texi +++ b/manual/math.texi @@ -240,6 +240,19 @@ All these functions, including the @code{_Float@var{N}} and should be prepared to cope with their absence. @end deftypefun +@deftypefun double cospi (double @var{x}) +@deftypefunx float cospif (float @var{x}) +@deftypefunx {long double} cospil (long double @var{x}) +@deftypefunx _FloatN cospifN (_Float@var{N} @var{x}) +@deftypefunx _FloatNx cospifNx (_Float@var{N}x @var{x}) +@standards{TS 18661-4:2015, math.h} +@safety{@prelim{}@mtsafe{}@assafe{}@acsafe{}} +These functions return the cosine of pi multiplied by @var{x}. The +return value is in the range @code{-1} to @code{1}. + +The @code{cospi} functions are from TS 18661-4:2015. +@end deftypefun + @cindex complex trigonometric functions @w{ISO C99} defines variants of the trig functions which work on diff --git a/math/Makefile b/math/Makefile index 2d5e016496..602797584e 100644 --- a/math/Makefile +++ b/math/Makefile @@ -105,6 +105,7 @@ gen-libm-calls = \ s_cexpF \ s_clog10F \ s_clogF \ + s_cospiF \ s_cpowF \ s_cprojF \ s_csinF \ @@ -635,6 +636,7 @@ libm-test-funcs-auto = \ clog10 \ cos \ cosh \ + cospi \ cpow \ csin \ csinh \ @@ -932,6 +934,7 @@ tgmath3-macros = \ copysign \ cos \ cosh \ + cospi \ cproj \ creal \ erf \ @@ -1324,6 +1327,7 @@ CFLAGS-s_conj.c += -fno-builtin-conjl CFLAGS-s_copysign.c += -fno-builtin-copysignl CFLAGS-s_cos.c += -fno-builtin-cosl CFLAGS-w_cosh.c += -fno-builtin-coshl +CFLAGS-s_cospi.c += -fno-builtin-cospil CFLAGS-s_cpow.c += -fno-builtin-cpowl CFLAGS-s_cproj.c += -fno-builtin-cprojl CFLAGS-s_creal.c += -fno-builtin-creall @@ -1462,6 +1466,7 @@ CFLAGS-s_conj.c += -fno-builtin-conjf32x -fno-builtin-conjf64 CFLAGS-s_copysign.c += -fno-builtin-copysignf32x -fno-builtin-copysignf64 CFLAGS-s_cos.c += -fno-builtin-cosf32x -fno-builtin-cosf64 CFLAGS-w_cosh.c += -fno-builtin-coshf32x -fno-builtin-coshf64 +CFLAGS-s_cospi.c += -fno-builtin-cospif32x -fno-builtin-cospif64 CFLAGS-s_cpow.c += -fno-builtin-cpowf32x -fno-builtin-cpowf64 CFLAGS-s_cproj.c += -fno-builtin-cprojf32x -fno-builtin-cprojf64 CFLAGS-s_creal.c += -fno-builtin-crealf32x -fno-builtin-crealf64 @@ -1586,6 +1591,7 @@ CFLAGS-s_conjf.c += -fno-builtin-conjf32 CFLAGS-s_copysignf.c += -fno-builtin-copysignf32 CFLAGS-s_cosf.c += -fno-builtin-cosf32 CFLAGS-w_coshf.c += -fno-builtin-coshf32 +CFLAGS-s_cospif.c += -fno-builtin-cospif32 CFLAGS-s_cpowf.c += -fno-builtin-cpowf32 CFLAGS-s_cprojf.c += -fno-builtin-cprojf32 CFLAGS-s_crealf.c += -fno-builtin-crealf32 diff --git a/math/Versions b/math/Versions index 49d61f212f..72dff5ff5b 100644 --- a/math/Versions +++ b/math/Versions @@ -653,4 +653,10 @@ libm { log10p1f64x; log10p1f128; logp1f64x; logp1f128; } + GLIBC_2.41 { + # Functions not involving _Float64x or _Float128, for all configurations. + cospi; cospif; cospil; cospif32; cospif64; cospif32x; + # Functions involving _Float64x or _Float128, for some configurations. + cospif64x; cospif128; + } } diff --git a/math/auto-libm-test-in b/math/auto-libm-test-in index b570e88e78..04fe31b863 100644 --- a/math/auto-libm-test-in +++ b/math/auto-libm-test-in @@ -4463,6 +4463,97 @@ cos 0x2.3c6ef4p-12 # the next value generates larger error bounds on x86_64 (binary128) cos 0xe.6672d458b05edf50af4fab1a42p+40 +cospi 0 +cospi -0 +cospi min +cospi -min +cospi min_subnorm +cospi -min_subnorm +cospi max +cospi -max +cospi 1 +cospi -1 +cospi 2 +cospi -2 +cospi 3 +cospi -3 +cospi 0.1 +cospi -0.1 +cospi 0.2 +cospi -0.2 +cospi 0.3 +cospi -0.3 +cospi 0.4 +cospi -0.4 +cospi 0.5 +cospi -0.5 +cospi 0.6 +cospi -0.6 +cospi 0.7 +cospi -0.7 +cospi 0.8 +cospi -0.8 +cospi 0.9 +cospi -0.9 +cospi 1.1 +cospi -1.1 +cospi 1.2 +cospi -1.2 +cospi 1.3 +cospi -1.3 +cospi 1.4 +cospi -1.4 +cospi 1.5 +cospi -1.5 +cospi 2.5 +cospi -2.5 +cospi 3.5 +cospi -3.5 +cospi 4.5 +cospi -4.5 +cospi 0.24 +cospi -0.24 +cospi 0.49 +cospi -0.49 +cospi 0.51 +cospi -0.51 +cospi 0.76 +cospi -0.76 +cospi 1.24 +cospi -1.24 +cospi 1.49 +cospi -1.49 +cospi 1.51 +cospi -1.51 +cospi 1.76 +cospi -1.76 +cospi 100.24 +cospi -100.24 +cospi 100.49 +cospi -100.49 +cospi 100.51 +cospi -100.51 +cospi 100.76 +cospi -100.76 +cospi 1234.56 +cospi -1234.56 +cospi 0x1fffff.6p0 +cospi -0x1fffff.6p0 +cospi 0x1fffff.ap0 +cospi -0x1fffff.ap0 +cospi 0x3ffffffffffff.6p0 +cospi -0x3ffffffffffff.6p0 +cospi 0x3ffffffffffff.ap0 +cospi -0x3ffffffffffff.ap0 +cospi 0x1fffffffffffffff.6p0 +cospi -0x1fffffffffffffff.6p0 +cospi 0x1fffffffffffffff.ap0 +cospi -0x1fffffffffffffff.ap0 +cospi 0x3fffffffffffffffffffffffffff.6p0 +cospi -0x3fffffffffffffffffffffffffff.6p0 +cospi 0x3fffffffffffffffffffffffffff.ap0 +cospi -0x3fffffffffffffffffffffffffff.ap0 + cosh 0 cosh -0 cosh 0.75 diff --git a/math/auto-libm-test-out-cospi b/math/auto-libm-test-out-cospi new file mode 100644 index 0000000000..cc6db07614 --- /dev/null +++ b/math/auto-libm-test-out-cospi @@ -0,0 +1,9258 @@ +cospi 0 += cospi downward binary32 0x0p+0 : 0x1p+0 : inexact-ok += cospi tonearest binary32 0x0p+0 : 0x1p+0 : inexact-ok += cospi towardzero binary32 0x0p+0 : 0x1p+0 : inexact-ok += cospi upward binary32 0x0p+0 : 0x1p+0 : inexact-ok += cospi downward binary64 0x0p+0 : 0x1p+0 : inexact-ok += cospi tonearest binary64 0x0p+0 : 0x1p+0 : inexact-ok += cospi towardzero binary64 0x0p+0 : 0x1p+0 : inexact-ok += cospi upward binary64 0x0p+0 : 0x1p+0 : inexact-ok += cospi downward intel96 0x0p+0 : 0x1p+0 : inexact-ok += cospi tonearest intel96 0x0p+0 : 0x1p+0 : inexact-ok += cospi towardzero intel96 0x0p+0 : 0x1p+0 : inexact-ok += cospi upward intel96 0x0p+0 : 0x1p+0 : inexact-ok += cospi downward m68k96 0x0p+0 : 0x1p+0 : inexact-ok += cospi tonearest m68k96 0x0p+0 : 0x1p+0 : inexact-ok += cospi towardzero m68k96 0x0p+0 : 0x1p+0 : inexact-ok += cospi upward m68k96 0x0p+0 : 0x1p+0 : inexact-ok += cospi downward binary128 0x0p+0 : 0x1p+0 : inexact-ok += cospi tonearest binary128 0x0p+0 : 0x1p+0 : inexact-ok += cospi towardzero binary128 0x0p+0 : 0x1p+0 : inexact-ok += cospi upward binary128 0x0p+0 : 0x1p+0 : inexact-ok += cospi downward ibm128 0x0p+0 : 0x1p+0 : inexact-ok += cospi tonearest ibm128 0x0p+0 : 0x1p+0 : inexact-ok += cospi towardzero ibm128 0x0p+0 : 0x1p+0 : inexact-ok += cospi upward ibm128 0x0p+0 : 0x1p+0 : inexact-ok +cospi -0 += cospi downward binary32 -0x0p+0 : 0x1p+0 : inexact-ok += cospi tonearest binary32 -0x0p+0 : 0x1p+0 : inexact-ok += cospi towardzero binary32 -0x0p+0 : 0x1p+0 : inexact-ok += cospi upward binary32 -0x0p+0 : 0x1p+0 : inexact-ok += cospi downward binary64 -0x0p+0 : 0x1p+0 : inexact-ok += cospi tonearest binary64 -0x0p+0 : 0x1p+0 : inexact-ok += cospi towardzero binary64 -0x0p+0 : 0x1p+0 : inexact-ok += cospi upward binary64 -0x0p+0 : 0x1p+0 : inexact-ok += cospi downward intel96 -0x0p+0 : 0x1p+0 : inexact-ok += cospi tonearest intel96 -0x0p+0 : 0x1p+0 : inexact-ok += cospi towardzero intel96 -0x0p+0 : 0x1p+0 : inexact-ok += cospi upward intel96 -0x0p+0 : 0x1p+0 : inexact-ok += cospi downward m68k96 -0x0p+0 : 0x1p+0 : inexact-ok += cospi tonearest m68k96 -0x0p+0 : 0x1p+0 : inexact-ok += cospi towardzero m68k96 -0x0p+0 : 0x1p+0 : inexact-ok += cospi upward m68k96 -0x0p+0 : 0x1p+0 : inexact-ok += cospi downward binary128 -0x0p+0 : 0x1p+0 : inexact-ok += cospi tonearest binary128 -0x0p+0 : 0x1p+0 : inexact-ok += cospi towardzero binary128 -0x0p+0 : 0x1p+0 : inexact-ok += cospi upward binary128 -0x0p+0 : 0x1p+0 : inexact-ok += cospi downward ibm128 -0x0p+0 : 0x1p+0 : inexact-ok += cospi tonearest ibm128 -0x0p+0 : 0x1p+0 : inexact-ok += cospi towardzero ibm128 -0x0p+0 : 0x1p+0 : inexact-ok += cospi upward ibm128 -0x0p+0 : 0x1p+0 : inexact-ok +cospi min += cospi downward binary32 0x4p-128 : 0xf.fffffp-4 : inexact-ok += cospi tonearest binary32 0x4p-128 : 0x1p+0 : inexact-ok += cospi towardzero binary32 0x4p-128 : 0xf.fffffp-4 : inexact-ok += cospi upward binary32 0x4p-128 : 0x1p+0 : inexact-ok += cospi downward binary64 0x4p-128 : 0xf.ffffffffffff8p-4 : inexact-ok += cospi tonearest binary64 0x4p-128 : 0x1p+0 : inexact-ok += cospi towardzero binary64 0x4p-128 : 0xf.ffffffffffff8p-4 : inexact-ok += cospi upward binary64 0x4p-128 : 0x1p+0 : inexact-ok += cospi downward intel96 0x4p-128 : 0xf.fffffffffffffffp-4 : inexact-ok += cospi tonearest intel96 0x4p-128 : 0x1p+0 : inexact-ok += cospi towardzero intel96 0x4p-128 : 0xf.fffffffffffffffp-4 : inexact-ok += cospi upward intel96 0x4p-128 : 0x1p+0 : inexact-ok += cospi downward m68k96 0x4p-128 : 0xf.fffffffffffffffp-4 : inexact-ok += cospi tonearest m68k96 0x4p-128 : 0x1p+0 : inexact-ok += cospi towardzero m68k96 0x4p-128 : 0xf.fffffffffffffffp-4 : inexact-ok += cospi upward m68k96 0x4p-128 : 0x1p+0 : inexact-ok += cospi downward binary128 0x4p-128 : 0xf.fffffffffffffffffffffffffff8p-4 : inexact-ok += cospi tonearest binary128 0x4p-128 : 0x1p+0 : inexact-ok += cospi towardzero binary128 0x4p-128 : 0xf.fffffffffffffffffffffffffff8p-4 : inexact-ok += cospi upward binary128 0x4p-128 : 0x1p+0 : inexact-ok += cospi downward ibm128 0x4p-128 : 0xf.fffffffffffffffffffffffffcp-4 : inexact-ok += cospi tonearest ibm128 0x4p-128 : 0x1p+0 : inexact-ok += cospi towardzero ibm128 0x4p-128 : 0xf.fffffffffffffffffffffffffcp-4 : inexact-ok += cospi upward ibm128 0x4p-128 : 0x1p+0 : inexact-ok += cospi downward binary64 0x4p-1024 : 0xf.ffffffffffff8p-4 : inexact-ok += cospi tonearest binary64 0x4p-1024 : 0x1p+0 : inexact-ok += cospi towardzero binary64 0x4p-1024 : 0xf.ffffffffffff8p-4 : inexact-ok += cospi upward binary64 0x4p-1024 : 0x1p+0 : inexact-ok += cospi downward intel96 0x4p-1024 : 0xf.fffffffffffffffp-4 : inexact-ok += cospi tonearest intel96 0x4p-1024 : 0x1p+0 : inexact-ok += cospi towardzero intel96 0x4p-1024 : 0xf.fffffffffffffffp-4 : inexact-ok += cospi upward intel96 0x4p-1024 : 0x1p+0 : inexact-ok += cospi downward m68k96 0x4p-1024 : 0xf.fffffffffffffffp-4 : inexact-ok += cospi tonearest m68k96 0x4p-1024 : 0x1p+0 : inexact-ok += cospi towardzero m68k96 0x4p-1024 : 0xf.fffffffffffffffp-4 : inexact-ok += cospi upward m68k96 0x4p-1024 : 0x1p+0 : inexact-ok += cospi downward binary128 0x4p-1024 : 0xf.fffffffffffffffffffffffffff8p-4 : inexact-ok += cospi tonearest binary128 0x4p-1024 : 0x1p+0 : inexact-ok += cospi towardzero binary128 0x4p-1024 : 0xf.fffffffffffffffffffffffffff8p-4 : inexact-ok += cospi upward binary128 0x4p-1024 : 0x1p+0 : inexact-ok += cospi downward ibm128 0x4p-1024 : 0xf.fffffffffffffffffffffffffcp-4 : inexact-ok += cospi tonearest ibm128 0x4p-1024 : 0x1p+0 : inexact-ok += cospi towardzero ibm128 0x4p-1024 : 0xf.fffffffffffffffffffffffffcp-4 : inexact-ok += cospi upward ibm128 0x4p-1024 : 0x1p+0 : inexact-ok += cospi downward intel96 0x4p-16 |