summaryrefslogtreecommitdiff
path: root/backend/src
diff options
context:
space:
mode:
Diffstat (limited to 'backend/src')
-rw-r--r--backend/src/libocl/tmpl/ocl_math_common.tmpl.cl50
1 files changed, 10 insertions, 40 deletions
diff --git a/backend/src/libocl/tmpl/ocl_math_common.tmpl.cl b/backend/src/libocl/tmpl/ocl_math_common.tmpl.cl
index c0ab2514..6026629b 100644
--- a/backend/src/libocl/tmpl/ocl_math_common.tmpl.cl
+++ b/backend/src/libocl/tmpl/ocl_math_common.tmpl.cl
@@ -201,38 +201,19 @@ OVERLOADABLE float inline __gen_ocl_internal_log_valid(float x) {
k += (i>>23);
f = x - 1.0f;
fsq = f * f;
-
- if((0x007fffff & (15 + ix)) < 16) { /* |f| < 2**-20 */
- R = fsq * (0.5f - 0.33333333333333333f * f);
- return k * ln2_hi + k * ln2_lo + f - R;
- }
-
- s = f / (2.0f + f);
+ s = mad(-2.0f, 1.0f / (2.0f + f), 1.0f);
z = s * s;
- i = ix - (0x6147a << 3);
w = z * z;
- j = (0x6b851 << 3) - ix;
- t1= w * mad(w, Lg4, Lg2);
- t2= z * mad(w, Lg3, Lg1);
- i |= j;
- R = t2 + t1;
- partial = (i > 0) ? -mad(s, 0.5f * fsq, -0.5f * fsq) : (s * f);
-
- return mad(s, R, f) - partial + k * ln2_hi + k * ln2_lo;;
+ t1 = w * mad(w, Lg4, Lg2);
+ R = mad(z, mad(w, Lg3, Lg1), t1);
+ w = 0.5f * fsq;
+ partial = -mad(s, w, -w);
+ return mad(k, ln2_lo, mad(k, ln2_hi, mad(s, R, f) - partial));
}
OVERLOADABLE float __gen_ocl_internal_log(float x)
{
- union { unsigned int i; float f; } u;
- u.f = x;
- int ix = u.i;
-
- if (ix < 0 )
- return NAN; /* log(-#) = NaN */
- if (ix >= 0x7f800000)
- return NAN;
-
- return __gen_ocl_internal_log_valid(x);
+ return __gen_ocl_internal_log_valid(x) + (x - x);
}
OVERLOADABLE float __gen_ocl_internal_log10(float x)
@@ -244,12 +225,10 @@ OVERLOADABLE float __gen_ocl_internal_log10(float x)
log10_2lo = 7.9034151668e-07; /* 0x355427db */
float y, z;
- int i, k, hx;
+ int i, k;
+ unsigned int hx;
u.f = x; hx = u.i;
-
- if (hx<0)
- return NAN; /* log(-#) = NaN */
if (hx >= 0x7f800000)
return NAN;
@@ -267,17 +246,8 @@ OVERLOADABLE float __gen_ocl_internal_log2(float x)
{
const float zero = 0.0,
invln2 = 0x1.715476p+0f;
- int ix;
-
- union { float f; int i; } u;
- u.f = x; ix = u.i;
-
- if (ix < 0)
- return NAN; /** log(-#) = NaN */
- if (ix >= 0x7f800000)
- return NAN;
- return invln2 * __gen_ocl_internal_log_valid(x);
+ return invln2 * __gen_ocl_internal_log_valid(x) + (x - x);
}