Added sin(pi*x) cos(pi*x) and tan(pi*x)

Author: ZERICO2005

src/libc/include/__math_def.h

@@ -110,6 +110,10 @@ double      cos(double);
 float       cosf(float);
 long double cosl(long double);
 
+double      cospi(double);
+float       cospif(float);
+long double cospil(long double);
+
 double      cosh(double);
 float       coshf(float);
 long double coshl(long double);
@@ -289,6 +293,10 @@ double      sin(double);
 float       sinf(float);
 long double sinl(long double);
 
+double      sinpi(double);
+float       sinpif(float);
+long double sinpil(long double);
+
 double      sinh(double);
 float       sinhf(float);
 long double sinhl(long double);
@@ -301,6 +309,10 @@ double      tan(double);
 float       tanf(float);
 long double tanl(long double);
 
+double      tanpi(double);
+float       tanpif(float);
+long double tanpil(long double);
+
 double      tanh(double);
 float       tanhf(float);
 long double tanhl(long double);

src/libc/include/tgmath.h

@@ -71,6 +71,12 @@
     float: cosf \
 )(x)
 
+#define cospi(x) _Generic(__tgmath_promote(x), \
+    long double: cospil, \
+    default: cospi, \
+    float: cospif \
+)(x)
+
 #define cosh(x) _Generic(__tgmath_promote(x), \
     long double: coshl, \
     default: cosh, \
@@ -227,6 +233,12 @@
     float: sinf \
 )(x)
 
+#define sinpi(x) _Generic(__tgmath_promote(x), \
+    long double: sinpil, \
+    default: sinpi, \
+    float: sinpif \
+)(x)
+
 #define sinh(x) _Generic(__tgmath_promote(x), \
     long double: sinhl, \
     default: sinh, \
@@ -245,6 +257,12 @@
     float: tanf \
 )(x)
 
+#define tanpi(x) _Generic(__tgmath_promote(x), \
+    long double: tanpil, \
+    default: tanpi, \
+    float: tanpif \
+)(x)
+
 #define tanh(x) _Generic(__tgmath_promote(x), \
     long double: tanhl, \
     default: tanh, \

src/libc/trigpif.c

@@ -0,0 +1,23 @@
+#include <math.h>
+#include "__float32_constants.h"
+
+static float reduce_f32(float x) {
+    // x is (-2.0, +2.0)
+    return (x - 2.0f * truncf(x * 0.5f)) * F32_PI;
+}
+
+float sinpif(float x) {
+    return sinf(reduce_f32(x));
+}
+
+float cospif(float x) {
+    return cosf(reduce_f32(x));
+}
+
+float tanpif(float x) {
+    return tanf(reduce_f32(x));
+}
+
+double sinpi(double) __attribute__((alias("sinpif")));
+double cospi(double) __attribute__((alias("cospif")));
+double tanpi(double) __attribute__((alias("tanpif")));

src/libc/trigpil.c

@@ -0,0 +1,19 @@
+#include <math.h>
+#include "__float64_constants.h"
+
+static long double reduce_f64(long double x) {
+    // x is (-2.0, +2.0)
+    return (x - 2.0L * truncl(x * 0.5L)) * F64_PI;
+}
+
+long double sinpil(long double x) {
+    return sinl(reduce_f64(x));
+}
+
+long double cospil(long double x) {
+    return cosl(reduce_f64(x));
+}
+
+long double tanpil(long double x) {
+    return tanl(reduce_f64(x));
+}

src/libcxx/include/cmath

@@ -209,6 +209,15 @@ template<typename _Tp> inline constexpr
 __cmath_enable_if_t<__cmath_is_integral_v<_Tp>, double>
 cos(_Tp __x) { return cos(__x); }
 
+using ::cospi;
+using ::cospif;
+using ::cospil;
+inline float cospi(float __x) { return cospif(__x); }
+inline long double cospi(long double __x) { return cospil(__x); }
+template<typename _Tp> inline
+__cmath_enable_if_t<__cmath_is_integral_v<_Tp>, double>
+cospi(_Tp __x) { return cospi(__x); }
+
 using ::cosh;
 using ::coshf;
 using ::coshl;
@@ -594,6 +603,15 @@ template<typename _Tp> inline constexpr
 __cmath_enable_if_t<__cmath_is_integral_v<_Tp>, double>
 sin(_Tp __x) { return sin(__x); }
 
+using ::sinpi;
+using ::sinpif;
+using ::sinpil;
+inline float sinpi(float __x) { return sinpif(__x); }
+inline long double sinpi(long double __x) { return sinpil(__x); }
+template<typename _Tp> inline
+__cmath_enable_if_t<__cmath_is_integral_v<_Tp>, double>
+sinpi(_Tp __x) { return sinpi(__x); }
+
 using ::sinh;
 using ::sinhf;
 using ::sinhl;
@@ -621,6 +639,15 @@ template<typename _Tp> inline constexpr
 __cmath_enable_if_t<__cmath_is_integral_v<_Tp>, double>
 tan(_Tp __x) { return tan(__x); }
 
+using ::tanpi;
+using ::tanpif;
+using ::tanpil;
+inline float tanpi(float __x) { return tanpif(__x); }
+inline long double tanpi(long double __x) { return tanpil(__x); }
+template<typename _Tp> inline
+__cmath_enable_if_t<__cmath_is_integral_v<_Tp>, double>
+tanpi(_Tp __x) { return tanpi(__x); }
+
 using ::tanh;
 using ::tanhf;
 using ::tanhl;