aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
Diffstat
-rw-r--r--SHARED-FILES4
-rw-r--r--sysdeps/aarch64/libm-test-ulps4
-rw-r--r--sysdeps/arc/fpu/libm-test-ulps4
-rw-r--r--sysdeps/arc/nofpu/libm-test-ulps1
-rw-r--r--sysdeps/arm/libm-test-ulps4
-rw-r--r--sysdeps/hppa/fpu/libm-test-ulps4
-rw-r--r--sysdeps/i386/fpu/libm-test-ulps4
-rw-r--r--sysdeps/i386/i686/fpu/multiarch/libm-test-ulps4
-rw-r--r--sysdeps/ieee754/flt-32/s_sinpif.c134
-rw-r--r--sysdeps/loongarch/lp64/libm-test-ulps4
-rw-r--r--sysdeps/mips/mips64/libm-test-ulps4
-rw-r--r--sysdeps/or1k/fpu/libm-test-ulps4
-rw-r--r--sysdeps/or1k/nofpu/libm-test-ulps1
-rw-r--r--sysdeps/powerpc/fpu/libm-test-ulps4
-rw-r--r--sysdeps/riscv/nofpu/libm-test-ulps1
-rw-r--r--sysdeps/riscv/rvd/libm-test-ulps4
-rw-r--r--sysdeps/s390/fpu/libm-test-ulps4
-rw-r--r--sysdeps/sparc/fpu/libm-test-ulps4
-rw-r--r--sysdeps/x86_64/fpu/libm-test-ulps4
19 files changed, 138 insertions, 59 deletions
diff --git a/SHARED-FILES b/SHARED-FILES
index 3ce38d1542..c108f3b308 100644
--- a/SHARED-FILES
+++ b/SHARED-FILES
@@ -354,3 +354,7 @@ sysdeps/ieee754/flt-32/s_cospif.c:
(src/binary32/cospi/cospif.c in CORE-MATH)
- the code was adapted to use glibc code style and internal
functions to handle errno, overflow, and underflow.
+sysdeps/ieee754/flt-32/s_sinpif.c:
+ (src/binary32/sinpi/sinpif.c in CORE-MATH)
+ - the code was adapted to use glibc code style and internal
+ functions to handle errno, overflow, and underflow.
diff --git a/sysdeps/aarch64/libm-test-ulps b/sysdeps/aarch64/libm-test-ulps
index a15f3effa0..c6c93aa0e4 100644
--- a/sysdeps/aarch64/libm-test-ulps
+++ b/sysdeps/aarch64/libm-test-ulps
@@ -1610,7 +1610,6 @@ ldouble: 4
Function: "sinpi":
double: 2
-float: 2
ldouble: 2
Function: "sinpi_advsimd":
@@ -1619,7 +1618,6 @@ float: 2
Function: "sinpi_downward":
double: 2
-float: 2
ldouble: 2
Function: "sinpi_sve":
@@ -1628,12 +1626,10 @@ float: 2
Function: "sinpi_towardzero":
double: 2
-float: 1
ldouble: 2
Function: "sinpi_upward":
double: 2
-float: 3
ldouble: 2
Function: "tan":
diff --git a/sysdeps/arc/fpu/libm-test-ulps b/sysdeps/arc/fpu/libm-test-ulps
index f2acbf453e..65ebf6f9a0 100644
--- a/sysdeps/arc/fpu/libm-test-ulps
+++ b/sysdeps/arc/fpu/libm-test-ulps
@@ -1101,19 +1101,15 @@ double: 3
Function: "sinpi":
double: 2
-float: 2
Function: "sinpi_downward":
double: 2
-float: 2
Function: "sinpi_towardzero":
double: 2
-float: 1
Function: "sinpi_upward":
double: 2
-float: 3
Function: "tan":
double: 1
diff --git a/sysdeps/arc/nofpu/libm-test-ulps b/sysdeps/arc/nofpu/libm-test-ulps
index 8716e5d29e..3ba4f01cbf 100644
--- a/sysdeps/arc/nofpu/libm-test-ulps
+++ b/sysdeps/arc/nofpu/libm-test-ulps
@@ -265,7 +265,6 @@ double: 2
Function: "sinpi":
double: 2
-float: 2
Function: "tanh":
double: 2
diff --git a/sysdeps/arm/libm-test-ulps b/sysdeps/arm/libm-test-ulps
index 647f92944c..f887712d8e 100644
--- a/sysdeps/arm/libm-test-ulps
+++ b/sysdeps/arm/libm-test-ulps
@@ -1097,19 +1097,15 @@ double: 3
Function: "sinpi":
double: 2
-float: 2
Function: "sinpi_downward":
double: 2
-float: 2
Function: "sinpi_towardzero":
double: 2
-float: 1
Function: "sinpi_upward":
double: 2
-float: 3
Function: "tan_downward":
double: 1
diff --git a/sysdeps/hppa/fpu/libm-test-ulps b/sysdeps/hppa/fpu/libm-test-ulps
index 88f7701c0e..10f7f2ebde 100644
--- a/sysdeps/hppa/fpu/libm-test-ulps
+++ b/sysdeps/hppa/fpu/libm-test-ulps
@@ -1123,19 +1123,15 @@ double: 3
Function: "sinpi":
double: 2
-float: 2
Function: "sinpi_downward":
double: 2
-float: 2
Function: "sinpi_towardzero":
double: 2
-float: 1
Function: "sinpi_upward":
double: 2
-float: 3
Function: "tan":
double: 1
diff --git a/sysdeps/i386/fpu/libm-test-ulps b/sysdeps/i386/fpu/libm-test-ulps
index 39066956b0..77aa7155db 100644
--- a/sysdeps/i386/fpu/libm-test-ulps
+++ b/sysdeps/i386/fpu/libm-test-ulps
@@ -1691,25 +1691,21 @@ ldouble: 5
Function: "sinpi":
double: 2
-float: 2
float128: 2
ldouble: 2
Function: "sinpi_downward":
double: 2
-float: 2
float128: 2
ldouble: 2
Function: "sinpi_towardzero":
double: 2
-float: 1
float128: 2
ldouble: 2
Function: "sinpi_upward":
double: 2
-float: 3
float128: 2
ldouble: 2
diff --git a/sysdeps/i386/i686/fpu/multiarch/libm-test-ulps b/sysdeps/i386/i686/fpu/multiarch/libm-test-ulps
index a8c4723850..7168d577d8 100644
--- a/sysdeps/i386/i686/fpu/multiarch/libm-test-ulps
+++ b/sysdeps/i386/i686/fpu/multiarch/libm-test-ulps
@@ -1696,25 +1696,21 @@ ldouble: 5
Function: "sinpi":
double: 2
-float: 2
float128: 2
ldouble: 2
Function: "sinpi_downward":
double: 2
-float: 2
float128: 2
ldouble: 3
Function: "sinpi_towardzero":
double: 2
-float: 1
float128: 2
ldouble: 2
Function: "sinpi_upward":
double: 2
-float: 3
float128: 2
ldouble: 2
diff --git a/sysdeps/ieee754/flt-32/s_sinpif.c b/sysdeps/ieee754/flt-32/s_sinpif.c
new file mode 100644
index 0000000000..772ea428e2
--- /dev/null
+++ b/sysdeps/ieee754/flt-32/s_sinpif.c
@@ -0,0 +1,134 @@
+/* Correctly-rounded sine of binary32 value for angles in half-revolutions
+
+Copyright (c) 2022-2025 Alexei Sibidanov.
+
+The original version of this file was copied from the CORE-MATH
+project (src/binary32/sinpi/sinpif.c, revision f786e13).
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.
+*/
+
+#include <math.h>
+#include <stdint.h>
+#include <errno.h>
+#include <fenv.h>
+#include <libm-alias-float.h>
+#include "math_config.h"
+
+float
+__sinpif (float x)
+{
+ static const double sn[] =
+ {
+ 0x1.921fb54442d0fp-37, -0x1.4abbce6102b94p-112, 0x1.4669fa3c58463p-189
+ };
+ static const double cn[] =
+ {
+ -0x1.3bd3cc9be45cfp-74, 0x1.03c1f08088742p-150, -0x1.55d1e5eff55a5p-228
+ };
+ static const double S[] =
+ {
+ 0x0p+0, 0x1.91f65f10dd814p-5, 0x1.917a6bc29b42cp-4,
+ 0x1.2c8106e8e613ap-3, 0x1.8f8b83c69a60bp-3, 0x1.f19f97b215f1bp-3,
+ 0x1.294062ed59f06p-2, 0x1.58f9a75ab1fddp-2, 0x1.87de2a6aea963p-2,
+ 0x1.b5d1009e15ccp-2, 0x1.e2b5d3806f63bp-2, 0x1.073879922ffeep-1,
+ 0x1.1c73b39ae68c8p-1, 0x1.30ff7fce17035p-1, 0x1.44cf325091dd6p-1,
+ 0x1.57d69348cecap-1, 0x1.6a09e667f3bcdp-1, 0x1.7b5df226aafafp-1,
+ 0x1.8bc806b151741p-1, 0x1.9b3e047f38741p-1, 0x1.a9b66290ea1a3p-1,
+ 0x1.b728345196e3ep-1, 0x1.c38b2f180bdb1p-1, 0x1.ced7af43cc773p-1,
+ 0x1.d906bcf328d46p-1, 0x1.e212104f686e5p-1, 0x1.e9f4156c62ddap-1,
+ 0x1.f0a7efb9230d7p-1, 0x1.f6297cff75cbp-1, 0x1.fa7557f08a517p-1,
+ 0x1.fd88da3d12526p-1, 0x1.ff621e3796d7ep-1, 0x1p+0,
+ 0x1.ff621e3796d7ep-1, 0x1.fd88da3d12526p-1, 0x1.fa7557f08a517p-1,
+ 0x1.f6297cff75cbp-1, 0x1.f0a7efb9230d7p-1, 0x1.e9f4156c62ddap-1,
+ 0x1.e212104f686e5p-1, 0x1.d906bcf328d46p-1, 0x1.ced7af43cc773p-1,
+ 0x1.c38b2f180bdb1p-1, 0x1.b728345196e3ep-1, 0x1.a9b66290ea1a3p-1,
+ 0x1.9b3e047f38741p-1, 0x1.8bc806b151741p-1, 0x1.7b5df226aafafp-1,
+ 0x1.6a09e667f3bcdp-1, 0x1.57d69348cecap-1, 0x1.44cf325091dd6p-1,
+ 0x1.30ff7fce17035p-1, 0x1.1c73b39ae68c8p-1, 0x1.073879922ffeep-1,
+ 0x1.e2b5d3806f63bp-2, 0x1.b5d1009e15ccp-2, 0x1.87de2a6aea963p-2,
+ 0x1.58f9a75ab1fddp-2, 0x1.294062ed59f06p-2, 0x1.f19f97b215f1bp-3,
+ 0x1.8f8b83c69a60bp-3, 0x1.2c8106e8e613ap-3, 0x1.917a6bc29b42cp-4,
+ 0x1.91f65f10dd814p-5, 0x0p+0, -0x1.91f65f10dd814p-5,
+ -0x1.917a6bc29b42cp-4, -0x1.2c8106e8e613ap-3, -0x1.8f8b83c69a60bp-3,
+ -0x1.f19f97b215f1bp-3, -0x1.294062ed59f06p-2, -0x1.58f9a75ab1fddp-2,
+ -0x1.87de2a6aea963p-2, -0x1.b5d1009e15ccp-2, -0x1.e2b5d3806f63bp-2,
+ -0x1.073879922ffeep-1, -0x1.1c73b39ae68c8p-1, -0x1.30ff7fce17035p-1,
+ -0x1.44cf325091dd6p-1, -0x1.57d69348cecap-1, -0x1.6a09e667f3bcdp-1,
+ -0x1.7b5df226aafafp-1, -0x1.8bc806b151741p-1, -0x1.9b3e047f38741p-1,
+ -0x1.a9b66290ea1a3p-1, -0x1.b728345196e3ep-1, -0x1.c38b2f180bdb1p-1,
+ -0x1.ced7af43cc773p-1, -0x1.d906bcf328d46p-1, -0x1.e212104f686e5p-1,
+ -0x1.e9f4156c62ddap-1, -0x1.f0a7efb9230d7p-1, -0x1.f6297cff75cbp-1,
+ -0x1.fa7557f08a517p-1, -0x1.fd88da3d12526p-1, -0x1.ff621e3796d7ep-1,
+ -0x1p+0, -0x1.ff621e3796d7ep-1, -0x1.fd88da3d12526p-1,
+ -0x1.fa7557f08a517p-1, -0x1.f6297cff75cbp-1, -0x1.f0a7efb9230d7p-1,
+ -0x1.e9f4156c62ddap-1, -0x1.e212104f686e5p-1, -0x1.d906bcf328d46p-1,
+ -0x1.ced7af43cc773p-1, -0x1.c38b2f180bdb1p-1, -0x1.b728345196e3ep-1,
+ -0x1.a9b66290ea1a3p-1, -0x1.9b3e047f38741p-1, -0x1.8bc806b151741p-1,
+ -0x1.7b5df226aafafp-1, -0x1.6a09e667f3bcdp-1, -0x1.57d69348cecap-1,
+ -0x1.44cf325091dd6p-1, -0x1.30ff7fce17035p-1, -0x1.1c73b39ae68c8p-1,
+ -0x1.073879922ffeep-1, -0x1.e2b5d3806f63bp-2, -0x1.b5d1009e15ccp-2,
+ -0x1.87de2a6aea963p-2, -0x1.58f9a75ab1fddp-2, -0x1.294062ed59f06p-2,
+ -0x1.f19f97b215f1bp-3, -0x1.8f8b83c69a60bp-3, -0x1.2c8106e8e613ap-3,
+ -0x1.917a6bc29b42cp-4, -0x1.91f65f10dd814p-5
+ };
+
+ uint32_t ix = asuint (x);
+ int32_t e = (ix >> 23) & 0xff;
+ if (__glibc_unlikely (e == 0xff))
+ {
+ if (!(ix << 9))
+ return __math_invalidf (x);
+ return x + x; /* nan */
+ }
+ int32_t m = (ix & ~0u >> 9) | 1 << 23, sgn = ix;
+ sgn >>= 31;
+ m = (m ^ sgn) - sgn;
+ int32_t s = 143 - e;
+ if (__glibc_unlikely (s < 0))
+ {
+ if (__glibc_unlikely (s < -6))
+ return copysignf (0.0f, x);
+ int32_t iq = m << (-s - 1);
+ iq &= 127;
+ if (iq == 0 || iq == 64)
+ return copysignf (0.0f, x);
+ return S[iq];
+ }
+ else if (__glibc_unlikely (s > 30))
+ {
+ double z = x, z2 = z * z;
+ return z * (0x1.921fb54442d18p+1 + z2 * (-0x1.4abbce625be53p+2));
+ }
+ int32_t si = 25 - s;
+ if (__glibc_unlikely (si >= 0 && (m << si) == 0))
+ return copysignf (0.0f, x);
+
+ int32_t k = m << (31 - s);
+ double z = k, z2 = z * z;
+ double fs = sn[0] + z2 * (sn[1] + z2 * sn[2]);
+ double fc = cn[0] + z2 * (cn[1] + z2 * cn[2]);
+ uint32_t iq = m >> s;
+ iq = (iq + 1) >> 1;
+ uint32_t is = iq & 127, ic = (iq + 32) & 127;
+ double ts = S[is], tc = S[ic];
+ double r = ts + (ts * z2) * fc + (tc * z) * fs;
+ return r;
+}
+libm_alias_float (__sinpi, sinpi)
diff --git a/sysdeps/loongarch/lp64/libm-test-ulps b/sysdeps/loongarch/lp64/libm-test-ulps
index fe84c60913..4fadba43c2 100644
--- a/sysdeps/loongarch/lp64/libm-test-ulps
+++ b/sysdeps/loongarch/lp64/libm-test-ulps
@@ -1389,22 +1389,18 @@ ldouble: 4
Function: "sinpi":
double: 2
-float: 2
ldouble: 2
Function: "sinpi_downward":
double: 2
-float: 2
ldouble: 2
Function: "sinpi_towardzero":
double: 2
-float: 1
ldouble: 2
Function: "sinpi_upward":
double: 2
-float: 3
ldouble: 2
Function: "tan":
diff --git a/sysdeps/mips/mips64/libm-test-ulps b/sysdeps/mips/mips64/libm-test-ulps
index ddc78d0239..5177b54557 100644
--- a/sysdeps/mips/mips64/libm-test-ulps
+++ b/sysdeps/mips/mips64/libm-test-ulps
@@ -1402,22 +1402,18 @@ ldouble: 4
Function: "sinpi":
double: 2
-float: 2
ldouble: 2
Function: "sinpi_downward":
double: 2
-float: 2
ldouble: 2
Function: "sinpi_towardzero":
double: 2
-float: 1
ldouble: 2
Function: "sinpi_upward":
double: 2
-float: 3
ldouble: 2
Function: "tan":
diff --git a/sysdeps/or1k/fpu/libm-test-ulps b/sysdeps/or1k/fpu/libm-test-ulps
index 884b4cc361..1fb4ec57c0 100644
--- a/sysdeps/or1k/fpu/libm-test-ulps
+++ b/sysdeps/or1k/fpu/libm-test-ulps
@@ -1082,19 +1082,15 @@ double: 3
Function: "sinpi":
double: 2
-float: 2
Function: "sinpi_downward":
double: 2
-float: 2
Function: "sinpi_towardzero":
double: 2
-float: 1
Function: "sinpi_upward":
double: 2
-float: 3
Function: "tan_downward":
double: 1
diff --git a/sysdeps/or1k/nofpu/libm-test-ulps b/sysdeps/or1k/nofpu/libm-test-ulps
index aec66e0fa3..aff11b5148 100644
--- a/sysdeps/or1k/nofpu/libm-test-ulps
+++ b/sysdeps/or1k/nofpu/libm-test-ulps
@@ -991,7 +991,6 @@ double: 3
Function: "sinpi":
double: 2
-float: 2
Function: "tan_downward":
double: 1
diff --git a/sysdeps/powerpc/fpu/libm-test-ulps b/sysdeps/powerpc/fpu/libm-test-ulps
index bdf0c98dc7..e59c3e47ef 100644
--- a/sysdeps/powerpc/fpu/libm-test-ulps
+++ b/sysdeps/powerpc/fpu/libm-test-ulps
@@ -1779,25 +1779,21 @@ ldouble: 6
Function: "sinpi":
double: 2
-float: 2
float128: 2
ldouble: 2
Function: "sinpi_downward":
double: 2
-float: 2
float128: 2
ldouble: 5
Function: "sinpi_towardzero":
double: 2
-float: 1
float128: 2
ldouble: 8
Function: "sinpi_upward":
double: 2
-float: 3
float128: 2
ldouble: 8
diff --git a/sysdeps/riscv/nofpu/libm-test-ulps b/sysdeps/riscv/nofpu/libm-test-ulps
index 08af2495f3..2545d0e166 100644
--- a/sysdeps/riscv/nofpu/libm-test-ulps
+++ b/sysdeps/riscv/nofpu/libm-test-ulps
@@ -1271,7 +1271,6 @@ ldouble: 4
Function: "sinpi":
double: 2
-float: 2
ldouble: 2
Function: "tan":
diff --git a/sysdeps/riscv/rvd/libm-test-ulps b/sysdeps/riscv/rvd/libm-test-ulps
index 6644e38ebc..94534a4f80 100644
--- a/sysdeps/riscv/rvd/libm-test-ulps
+++ b/sysdeps/riscv/rvd/libm-test-ulps
@@ -1401,22 +1401,18 @@ ldouble: 4
Function: "sinpi":
double: 2
-float: 2
ldouble: 2
Function: "sinpi_downward":
double: 2
-float: 2
ldouble: 2
Function: "sinpi_towardzero":
double: 2
-float: 1
ldouble: 2
Function: "sinpi_upward":
double: 2
-float: 3
ldouble: 2
Function: "tan":
diff --git a/sysdeps/s390/fpu/libm-test-ulps b/sysdeps/s390/fpu/libm-test-ulps
index 6318760eb5..2c89048b56 100644
--- a/sysdeps/s390/fpu/libm-test-ulps
+++ b/sysdeps/s390/fpu/libm-test-ulps
@@ -1387,22 +1387,18 @@ ldouble: 4
Function: "sinpi":
double: 2
-float: 2
ldouble: 2
Function: "sinpi_downward":
double: 2
-float: 2
ldouble: 2
Function: "sinpi_towardzero":
double: 2
-float: 1
ldouble: 2
Function: "sinpi_upward":
double: 2
-float: 3
ldouble: 2
Function: "tan":
diff --git a/sysdeps/sparc/fpu/libm-test-ulps b/sysdeps/sparc/fpu/libm-test-ulps
index 2c319f8ae2..3af2355545 100644
--- a/sysdeps/sparc/fpu/libm-test-ulps
+++ b/sysdeps/sparc/fpu/libm-test-ulps
@@ -1402,22 +1402,18 @@ ldouble: 4
Function: "sinpi":
double: 2
-float: 2
ldouble: 2
Function: "sinpi_downward":
double: 2
-float: 2
ldouble: 2
Function: "sinpi_towardzero":
double: 2
-float: 1
ldouble: 2
Function: "sinpi_upward":
double: 2
-float: 3
ldouble: 2
Function: "tan":
diff --git a/sysdeps/x86_64/fpu/libm-test-ulps b/sysdeps/x86_64/fpu/libm-test-ulps
index e2cf3e04b6..f6da5ba186 100644
--- a/sysdeps/x86_64/fpu/libm-test-ulps
+++ b/sysdeps/x86_64/fpu/libm-test-ulps
@@ -2197,25 +2197,21 @@ float: 1
Function: "sinpi":
double: 2
-float: 2
float128: 2
ldouble: 2
Function: "sinpi_downward":
double: 2
-float: 2
float128: 2
ldouble: 3
Function: "sinpi_towardzero":
double: 2
-float: 1
float128: 2
ldouble: 2
Function: "sinpi_upward":
double: 2
-float: 3
float128: 2
ldouble: 2
generated by cgit v1.2.3 (git 2.25.1) at 2026-01-07 14:52:17 +0000