diff options
Diffstat (limited to 'libclc/generic/lib')
| -rw-r--r-- | libclc/generic/lib/SOURCES | 1 | ||||
| -rw-r--r-- | libclc/generic/lib/math/logb.cl | 31 |
2 files changed, 32 insertions, 0 deletions
diff --git a/libclc/generic/lib/SOURCES b/libclc/generic/lib/SOURCES index 74bea25bbc2..517daba0ba3 100644 --- a/libclc/generic/lib/SOURCES +++ b/libclc/generic/lib/SOURCES @@ -104,6 +104,7 @@ math/log.cl math/log10.cl math/log1p.cl math/log2.cl +math/logb.cl math/mad.cl math/modf.cl math/native_log.cl diff --git a/libclc/generic/lib/math/logb.cl b/libclc/generic/lib/math/logb.cl new file mode 100644 index 00000000000..31e51616534 --- /dev/null +++ b/libclc/generic/lib/math/logb.cl @@ -0,0 +1,31 @@ +#include <clc/clc.h> +#include "math.h" +#include "../clcmacro.h" + +_CLC_OVERLOAD _CLC_DEF float logb(float x) { + int ax = as_int(x) & EXSIGNBIT_SP32; + float s = -118 - clz(ax); + float r = (ax >> EXPSHIFTBITS_SP32) - EXPBIAS_SP32; + r = ax >= PINFBITPATT_SP32 ? as_float(ax) : r; + r = ax < 0x00800000 ? s : r; + r = ax == 0 ? as_float(NINFBITPATT_SP32) : r; + return r; +} + +_CLC_UNARY_VECTORIZE(_CLC_OVERLOAD _CLC_DEF, float, logb, float); + +#ifdef cl_khr_fp64 +#pragma OPENCL EXTENSION cl_khr_fp64 : enable + +_CLC_OVERLOAD _CLC_DEF double logb(double x) { + long ax = as_long(x) & EXSIGNBIT_DP64; + double s = -1011L - clz(ax); + double r = (int) (ax >> EXPSHIFTBITS_DP64) - EXPBIAS_DP64; + r = ax >= PINFBITPATT_DP64 ? as_double(ax) : r; + r = ax < 0x0010000000000000L ? s : r; + r = ax == 0L ? as_double(NINFBITPATT_DP64) : r; + return r; +} + +_CLC_UNARY_VECTORIZE(_CLC_OVERLOAD _CLC_DEF, double, logb, double) +#endif |
