Author: jvesely Date: Wed May 6 16:59:29 2015 New Revision: 236648 URL: http://llvm.org/viewvc/llvm-project?rev=236648&view=rev Log: Move ldexp soft implementation to a separate file
Signed-off-by: Jan Vesely <[email protected]> Reviewed-by: Tom Stellard <[email protected]> Added: libclc/trunk/generic/include/math/clc_ldexp.h libclc/trunk/generic/lib/math/clc_ldexp.cl - copied, changed from r236647, libclc/trunk/generic/lib/math/ldexp.cl Modified: libclc/trunk/generic/lib/SOURCES libclc/trunk/generic/lib/math/ldexp.cl Added: libclc/trunk/generic/include/math/clc_ldexp.h URL: http://llvm.org/viewvc/llvm-project/libclc/trunk/generic/include/math/clc_ldexp.h?rev=236648&view=auto ============================================================================== --- libclc/trunk/generic/include/math/clc_ldexp.h (added) +++ libclc/trunk/generic/include/math/clc_ldexp.h Wed May 6 16:59:29 2015 @@ -0,0 +1,6 @@ +_CLC_DEF _CLC_OVERLOAD float __clc_ldexp(float, int); + +#ifdef cl_khr_fp64 + #pragma OPENCL EXTENSION cl_khr_fp64 : enable + _CLC_DEF _CLC_OVERLOAD float __clc_ldexp(double, int); +#endif Modified: libclc/trunk/generic/lib/SOURCES URL: http://llvm.org/viewvc/llvm-project/libclc/trunk/generic/lib/SOURCES?rev=236648&r1=236647&r2=236648&view=diff ============================================================================== --- libclc/trunk/generic/lib/SOURCES (original) +++ libclc/trunk/generic/lib/SOURCES Wed May 6 16:59:29 2015 @@ -84,6 +84,7 @@ math/fmod.cl math/fract.cl math/half_sqrt.cl math/hypot.cl +math/clc_ldexp.cl math/ldexp.cl math/log10.cl math/log1p.cl Copied: libclc/trunk/generic/lib/math/clc_ldexp.cl (from r236647, libclc/trunk/generic/lib/math/ldexp.cl) URL: http://llvm.org/viewvc/llvm-project/libclc/trunk/generic/lib/math/clc_ldexp.cl?p2=libclc/trunk/generic/lib/math/clc_ldexp.cl&p1=libclc/trunk/generic/lib/math/ldexp.cl&r1=236647&r2=236648&rev=236648&view=diff ============================================================================== --- libclc/trunk/generic/lib/math/ldexp.cl (original) +++ libclc/trunk/generic/lib/math/clc_ldexp.cl Wed May 6 16:59:29 2015 @@ -25,7 +25,7 @@ #include "../clcmacro.h" #include "math.h" -_CLC_DEF _CLC_OVERLOAD float ldexp(float x, int n) { +_CLC_DEF _CLC_OVERLOAD float __clc_ldexp(float x, int n) { if (!__clc_fp32_subnormals_supported()) { @@ -92,14 +92,11 @@ _CLC_DEF _CLC_OVERLOAD float ldexp(float return val_f; } -// This defines all the ldexp(floatN, intN) variants. -_CLC_BINARY_VECTORIZE(_CLC_OVERLOAD _CLC_DEF, float, ldexp, float, int) - #ifdef cl_khr_fp64 #pragma OPENCL EXTENSION cl_khr_fp64 : enable -_CLC_DEF _CLC_OVERLOAD double ldexp(double x, int n) { +_CLC_DEF _CLC_OVERLOAD double __clc_ldexp(double x, int n) { long l = as_ulong(x); int e = (l >> 52) & 0x7ff; long s = l & 0x8000000000000000; @@ -128,11 +125,4 @@ _CLC_DEF _CLC_OVERLOAD double ldexp(doub return mr; } -// This defines all the ldexp(doubleN, intN) variants. -_CLC_BINARY_VECTORIZE(_CLC_OVERLOAD _CLC_DEF, double, ldexp, double, int) - #endif - -// This defines all the ldexp(GENTYPE, int) variants -#define __CLC_BODY <ldexp.inc> -#include <clc/math/gentype.inc> Modified: libclc/trunk/generic/lib/math/ldexp.cl URL: http://llvm.org/viewvc/llvm-project/libclc/trunk/generic/lib/math/ldexp.cl?rev=236648&r1=236647&r2=236648&view=diff ============================================================================== --- libclc/trunk/generic/lib/math/ldexp.cl (original) +++ libclc/trunk/generic/lib/math/ldexp.cl Wed May 6 16:59:29 2015 @@ -24,112 +24,15 @@ #include "config.h" #include "../clcmacro.h" #include "math.h" +#include "math/clc_ldexp.h" -_CLC_DEF _CLC_OVERLOAD float ldexp(float x, int n) { - - if (!__clc_fp32_subnormals_supported()) { - - // This treats subnormals as zeros - int i = as_int(x); - int e = (i >> 23) & 0xff; - int m = i & 0x007fffff; - int s = i & 0x80000000; - int v = add_sat(e, n); - v = clamp(v, 0, 0xff); - int mr = e == 0 | v == 0 | v == 0xff ? 0 : m; - int c = e == 0xff; - mr = c ? m : mr; - int er = c ? e : v; - er = e ? er : e; - return as_float( s | (er << 23) | mr ); - } - - /* supports denormal values */ - const int multiplier = 24; - float val_f; - uint val_ui; - uint sign; - int exponent; - val_ui = as_uint(x); - sign = val_ui & 0x80000000; - val_ui = val_ui & 0x7fffffff;/* remove the sign bit */ - int val_x = val_ui; - - exponent = val_ui >> 23; /* get the exponent */ - int dexp = exponent; - - /* denormal support */ - int fbh = 127 - (as_uint((float)(as_float(val_ui | 0x3f800000) - 1.0f)) >> 23); - int dexponent = 25 - fbh; - uint dval_ui = (( (val_ui << fbh) & 0x007fffff) | (dexponent << 23)); - int ex = dexponent + n - multiplier; - dexponent = ex; - uint val = sign | (ex << 23) | (dval_ui & 0x007fffff); - int ex1 = dexponent + multiplier; - ex1 = -ex1 +25; - dval_ui = (((dval_ui & 0x007fffff )| 0x800000) >> ex1); - dval_ui = dexponent > 0 ? val :dval_ui; - dval_ui = dexponent > 254 ? 0x7f800000 :dval_ui; /*overflow*/ - dval_ui = dexponent < -multiplier ? 0 : dval_ui; /*underflow*/ - dval_ui = dval_ui | sign; - val_f = as_float(dval_ui); - - exponent += n; - - val = sign | (exponent << 23) | (val_ui & 0x007fffff); - ex1 = exponent + multiplier; - ex1 = -ex1 +25; - val_ui = (((val_ui & 0x007fffff )| 0x800000) >> ex1); - val_ui = exponent > 0 ? val :val_ui; - val_ui = exponent > 254 ? 0x7f800000 :val_ui; /*overflow*/ - val_ui = exponent < -multiplier ? 0 : val_ui; /*underflow*/ - val_ui = val_ui | sign; - - val_ui = dexp == 0? dval_ui : val_ui; - val_f = as_float(val_ui); - - val_f = isnan(x) | isinf(x) | val_x == 0 ? x : val_f; - return val_f; -} - -// This defines all the ldexp(floatN, intN) variants. -_CLC_BINARY_VECTORIZE(_CLC_OVERLOAD _CLC_DEF, float, ldexp, float, int) +_CLC_DEFINE_BINARY_BUILTIN(float, ldexp, __clc_ldexp, float, int) #ifdef cl_khr_fp64 #pragma OPENCL EXTENSION cl_khr_fp64 : enable -_CLC_DEF _CLC_OVERLOAD double ldexp(double x, int n) { - long l = as_ulong(x); - int e = (l >> 52) & 0x7ff; - long s = l & 0x8000000000000000; - - ulong ux = as_ulong(x * 0x1.0p+53); - int de = ((int)(ux >> 52) & 0x7ff) - 53; - int c = e == 0; - e = c ? de: e; - - ux = c ? ux : l; - - int v = e + n; - v = clamp(v, -0x7ff, 0x7ff); - - ux &= ~EXPBITS_DP64; - - double mr = as_double(ux | ((ulong)(v+53) << 52)); - mr = mr * 0x1.0p-53; - - mr = v > 0 ? as_double(ux | ((ulong)v << 52)) : mr; - - mr = v == 0x7ff ? as_double(s | PINFBITPATT_DP64) : mr; - mr = v < -53 ? as_double(s) : mr; - - mr = ((n == 0) | isinf(x) | (x == 0) ) ? x : mr; - return mr; -} - -// This defines all the ldexp(doubleN, intN) variants. -_CLC_BINARY_VECTORIZE(_CLC_OVERLOAD _CLC_DEF, double, ldexp, double, int) +_CLC_DEFINE_BINARY_BUILTIN(float, ldexp, __clc_ldexp, float, int) #endif _______________________________________________ cfe-commits mailing list [email protected] http://lists.cs.uiuc.edu/mailman/listinfo/cfe-commits
