aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
Diffstat
-rw-r--r--libclc/generic/lib/math/ep_log.cl3
-rw-r--r--libclc/generic/lib/math/ep_log.h3
-rw-r--r--libclc/generic/lib/math/modf.inc18
-rw-r--r--libclc/generic/lib/math/sincos_helpers.cl24
-rw-r--r--libclc/generic/lib/math/sincos_helpers.h12
5 files changed, 35 insertions, 25 deletions
diff --git a/libclc/generic/lib/math/ep_log.cl b/libclc/generic/lib/math/ep_log.cl
index 65db94a..f0b5d3f 100644
--- a/libclc/generic/lib/math/ep_log.cl
+++ b/libclc/generic/lib/math/ep_log.cl
@@ -38,7 +38,8 @@
#define LF1 1.24999999978138668903e-02
#define LF2 2.23219810758559851206e-03
-_CLC_DEF void __clc_ep_log(double x, int *xexp, double *r1, double *r2) {
+_CLC_DEF void __clc_ep_log(double x, private int *xexp, private double *r1,
+ private double *r2) {
// Computes natural log(x). Algorithm based on:
// Ping-Tak Peter Tang
// "Table-driven implementation of the logarithm function in IEEE
diff --git a/libclc/generic/lib/math/ep_log.h b/libclc/generic/lib/math/ep_log.h
index 414e623..3176cfe 100644
--- a/libclc/generic/lib/math/ep_log.h
+++ b/libclc/generic/lib/math/ep_log.h
@@ -26,6 +26,7 @@
#pragma OPENCL EXTENSION cl_khr_fp64 : enable
-_CLC_DECL void __clc_ep_log(double x, int *xexp, double *r1, double *r2);
+_CLC_DECL void __clc_ep_log(double x, private int *xexp, private double *r1,
+ private double *r2);
#endif
diff --git a/libclc/generic/lib/math/modf.inc b/libclc/generic/lib/math/modf.inc
index 1ffc6d9e..ff7ef30 100644
--- a/libclc/generic/lib/math/modf.inc
+++ b/libclc/generic/lib/math/modf.inc
@@ -28,18 +28,20 @@
#define ZERO 0.0h
#endif
-_CLC_OVERLOAD _CLC_DEF __CLC_GENTYPE modf(__CLC_GENTYPE x, __CLC_GENTYPE *iptr) {
+_CLC_OVERLOAD _CLC_DEF __CLC_GENTYPE modf(__CLC_GENTYPE x,
+ private __CLC_GENTYPE *iptr) {
*iptr = trunc(x);
return copysign(isinf(x) ? ZERO : x - *iptr, x);
}
-#define MODF_DEF(addrspace) \
- _CLC_OVERLOAD _CLC_DEF __CLC_GENTYPE modf(__CLC_GENTYPE x, addrspace __CLC_GENTYPE *iptr) { \
- __CLC_GENTYPE private_iptr; \
- __CLC_GENTYPE ret = modf(x, &private_iptr); \
- *iptr = private_iptr; \
- return ret; \
-}
+#define MODF_DEF(addrspace) \
+ _CLC_OVERLOAD _CLC_DEF __CLC_GENTYPE modf(__CLC_GENTYPE x, \
+ addrspace __CLC_GENTYPE *iptr) { \
+ __CLC_GENTYPE private_iptr; \
+ __CLC_GENTYPE ret = modf(x, &private_iptr); \
+ *iptr = private_iptr; \
+ return ret; \
+ }
MODF_DEF(local);
MODF_DEF(global);
diff --git a/libclc/generic/lib/math/sincos_helpers.cl b/libclc/generic/lib/math/sincos_helpers.cl
index 22f2bf61..441bad2 100644
--- a/libclc/generic/lib/math/sincos_helpers.cl
+++ b/libclc/generic/lib/math/sincos_helpers.cl
@@ -119,8 +119,8 @@ _CLC_DEF float __clc_tanf_piby4(float x, int regn) {
return regn & 1 ? tr : t;
}
-_CLC_DEF void __clc_fullMulS(float *hi, float *lo, float a, float b, float bh,
- float bt) {
+_CLC_DEF void __clc_fullMulS(private float *hi, private float *lo, float a,
+ float b, float bh, float bt) {
if (HAVE_HW_FMA32()) {
float ph = a * b;
*hi = ph;
@@ -136,7 +136,7 @@ _CLC_DEF void __clc_fullMulS(float *hi, float *lo, float a, float b, float bh,
}
}
-_CLC_DEF float __clc_removePi2S(float *hi, float *lo, float x) {
+_CLC_DEF float __clc_removePi2S(private float *hi, private float *lo, float x) {
// 72 bits of pi/2
const float fpiby2_1 = (float)0xC90FDA / 0x1.0p+23f;
const float fpiby2_1_h = (float)0xC90 / 0x1.0p+11f;
@@ -174,7 +174,8 @@ _CLC_DEF float __clc_removePi2S(float *hi, float *lo, float x) {
return fnpi2;
}
-_CLC_DEF int __clc_argReductionSmallS(float *r, float *rr, float x) {
+_CLC_DEF int __clc_argReductionSmallS(private float *r, private float *rr,
+ float x) {
float fnpi2 = __clc_removePi2S(r, rr, x);
return (int)fnpi2 & 0x3;
}
@@ -188,7 +189,8 @@ _CLC_DEF int __clc_argReductionSmallS(float *r, float *rr, float x) {
HI = __clc_mul_hi(A, B); \
HI += LO < C
-_CLC_DEF int __clc_argReductionLargeS(float *r, float *rr, float x) {
+_CLC_DEF int __clc_argReductionLargeS(private float *r, private float *rr,
+ float x) {
int xe = (int)(as_uint(x) >> 23) - 127;
uint xm = 0x00800000U | (as_uint(x) & 0x7fffffU);
@@ -330,7 +332,7 @@ _CLC_DEF int __clc_argReductionLargeS(float *r, float *rr, float x) {
return ((i >> 1) + (i & 1)) & 0x3;
}
-_CLC_DEF int __clc_argReductionS(float *r, float *rr, float x) {
+_CLC_DEF int __clc_argReductionS(private float *r, private float *rr, float x) {
if (x < 0x1.0p+23f)
return __clc_argReductionSmallS(r, rr, x);
else
@@ -342,8 +344,9 @@ _CLC_DEF int __clc_argReductionS(float *r, float *rr, float x) {
#pragma OPENCL EXTENSION cl_khr_fp64 : enable
// Reduction for medium sized arguments
-_CLC_DEF void __clc_remainder_piby2_medium(double x, double *r, double *rr,
- int *regn) {
+_CLC_DEF void __clc_remainder_piby2_medium(double x, private double *r,
+ private double *rr,
+ private int *regn) {
// How many pi/2 is x a multiple of?
const double two_by_pi = 0x1.45f306dc9c883p-1;
double dnpi2 = __clc_trunc(fma(x, two_by_pi, 0.5));
@@ -387,8 +390,9 @@ _CLC_DEF void __clc_remainder_piby2_medium(double x, double *r, double *rr,
// Return value "regn" tells how many lots of pi/2 were subtracted
// from x to put it in the range [-pi/4,pi/4], mod 4.
-_CLC_DEF void __clc_remainder_piby2_large(double x, double *r, double *rr,
- int *regn) {
+_CLC_DEF void __clc_remainder_piby2_large(double x, private double *r,
+ private double *rr,
+ private int *regn) {
long ux = as_long(x);
int e = (int)(ux >> 52) - 1023;
diff --git a/libclc/generic/lib/math/sincos_helpers.h b/libclc/generic/lib/math/sincos_helpers.h
index 6dbca73..c7981e5 100644
--- a/libclc/generic/lib/math/sincos_helpers.h
+++ b/libclc/generic/lib/math/sincos_helpers.h
@@ -26,16 +26,18 @@
_CLC_DECL float __clc_sinf_piby4(float x, float y);
_CLC_DECL float __clc_cosf_piby4(float x, float y);
_CLC_DECL float __clc_tanf_piby4(float x, int y);
-_CLC_DECL int __clc_argReductionS(float *r, float *rr, float x);
+_CLC_DECL int __clc_argReductionS(private float *r, private float *rr, float x);
#ifdef cl_khr_fp64
#pragma OPENCL EXTENSION cl_khr_fp64 : enable
-_CLC_DECL void __clc_remainder_piby2_medium(double x, double *r, double *rr,
- int *regn);
-_CLC_DECL void __clc_remainder_piby2_large(double x, double *r, double *rr,
- int *regn);
+_CLC_DECL void __clc_remainder_piby2_medium(double x, private double *r,
+ private double *rr,
+ private int *regn);
+_CLC_DECL void __clc_remainder_piby2_large(double x, private double *r,
+ private double *rr,
+ private int *regn);
_CLC_DECL double2 __clc_sincos_piby4(double x, double xx);
#endif
generated by cgit v1.1 at 2025-08-02 15:32:15 +0000