[Beignet] [PATCH] utest: add test case for builtin function exp/exp2/exp10/expm1.

[Lu, Guanqun]

I'm suspicious about the printf_c() part, do we really need this in one of these test case? it looks nice, but it clobbers the output when we're trying to do output redirection. :)

> -----Original Message-----
> From: beignet-bounces+guanqun.lu=intel.com at lists.freedesktop.org
> [mailto:beignet-bounces+guanqun.lu=intel.com at lists.freedesktop.org] On
> Behalf Of Yi Sun
> Sent: Thursday, October 24, 2013 3:57 PM
> To: beignet at lists.freedesktop.org
> Cc: Jin, Gordon; Sun, Yi; Shui, YangweiX
> Subject: [Beignet] [PATCH] utest: add test case for builtin function
> exp/exp2/exp10/expm1.
> 
> Signed-off-by: Yi Sun <yi.sun at intel.com>
> Signed-off-by: Yangwei Shui <yangweix.shui at intel.com>
> 
> diff --git a/kernels/builtin_exp.cl b/kernels/builtin_exp.cl
> new file mode 100644
> index 0000000..ecc1a3e
> --- /dev/null
> +++ b/kernels/builtin_exp.cl
> @@ -0,0 +1,10 @@
> +__kernel void builtin_exp(__global float *dst, __global float *src, __global int
> *max_func) {
> +  int i = get_global_id(0);
> +  float x = src[i];
> +
> +  dst[i * (*max_func) + 0] = exp(x);
> +  dst[i * (*max_func) + 1] = exp2(x);
> +  dst[i * (*max_func) + 2] = exp10(x);
> +  dst[i * (*max_func) + 3] = expm1(x);
> +  dst[i * (*max_func) + 4] = x;
> +};
> diff --git a/utests/CMakeLists.txt b/utests/CMakeLists.txt
> index da4cfd7..c7eb8a9 100644
> --- a/utests/CMakeLists.txt
> +++ b/utests/CMakeLists.txt
> @@ -121,6 +121,7 @@ set (utests_sources
>    builtin_local_id.cpp
>    builtin_acos_asin.cpp
>    builtin_pow.cpp
> +  builtin_exp.cpp
>    runtime_createcontext.cpp
>    runtime_null_kernel_arg.cpp
>  	runtime_event.cpp
> diff --git a/utests/builtin_exp.cpp b/utests/builtin_exp.cpp
> new file mode 100644
> index 0000000..d5288c8
> --- /dev/null
> +++ b/utests/builtin_exp.cpp
> @@ -0,0 +1,102 @@
> +#include "utest_helper.hpp"
> +#include <cmath>
> +#include <algorithm>
> +
> +#define udebug 0
> +
> +#define FLT_MAX 0x1.fffffep127f
> +#define FLT_MIN 0x1.0p-126f
> +#define FLT_ULP  (1.0e-6f)
> +
> +#define printf_c(...) \
> +{\
> +  printf("\033[1m\033[40;31m");\
> +  printf( __VA_ARGS__ );\
> +  printf("\033[0m");\
> +}
> +
> +const float input_data[] = {FLT_MAX, -FLT_MAX, FLT_MIN, -FLT_MIN, 80, -80,
> 3.14, -3.14, -0.5, 0.5, 1, -1, 0.0 };
> +const int count_input = sizeof(input_data) / sizeof(input_data[0]);
> +const int max_function = 5;
> +
> +static void cpu_compiler_math(float *dst, const float *src)
> +{
> +  const float x = *src;
> +
> +  dst[0] = exp(x);
> +  dst[1] = exp2(x);
> +  dst[2] = exp10(x);
> +  dst[3] = expm1(x);
> +  dst[4] = x;
> +}
> +
> +static void builtin_exp(void)
> +{
> +  // Setup kernel and buffers
> +  int k, i, index_cur;
> +  float gpu_data[max_function * count_input] = {0}, cpu_data[max_function *
> count_input] = {0};
> +  float diff;
> +  char log[256] = {0};
> +
> +  OCL_CREATE_KERNEL("builtin_exp");
> +
> +  OCL_CREATE_BUFFER(buf[0], CL_MEM_READ_WRITE, count_input *
> max_function * sizeof(float), NULL);
> +  OCL_CREATE_BUFFER(buf[1], CL_MEM_READ_WRITE, count_input *
> sizeof(float), NULL);
> +  OCL_CREATE_BUFFER(buf[2], CL_MEM_READ_WRITE, sizeof(int), NULL);
> +  OCL_SET_ARG(0, sizeof(cl_mem), &buf[0]);
> +  OCL_SET_ARG(1, sizeof(cl_mem), &buf[1]);
> +  OCL_SET_ARG(2, sizeof(cl_mem), &buf[2]);
> +
> +  globals[0] = count_input;
> +  locals[0] = 1;
> +
> +  clEnqueueWriteBuffer( queue, buf[1], CL_TRUE, 0, count_input *
> sizeof(float), input_data, 0, NULL, NULL);
> +  clEnqueueWriteBuffer( queue, buf[2], CL_TRUE, 0, sizeof(int),
> &max_function , 0, NULL, NULL);
> +
> +   // Run the kernel
> +  OCL_NDRANGE( 1 );
> +
> +  clEnqueueReadBuffer( queue, buf[0], CL_TRUE, 0, sizeof(float) *
> max_function * count_input, gpu_data, 0, NULL, NULL);
> +
> +  for (k = 0; (uint)k < count_input; k++)
> +  {
> +    cpu_compiler_math( cpu_data + k * max_function, input_data + k);
> +
> +    for (i = 0; i < max_function; i++)
> +    {
> +      index_cur = k * max_function + i;
> +      diff = fabs(gpu_data[index_cur]-cpu_data[index_cur]);
> +      sprintf(log, "%d/%d: %f -> gpu:%f  cpu:%f diff:%f expect:%f\n", \
> +         k, i, input_data[k], gpu_data[index_cur], cpu_data[index_cur], \
> +         diff/gpu_data[index_cur], 3 * FLT_ULP);
> +
> +#if udebug
> +      if (isinf(cpu_data[index_cur]) && isinf(gpu_data[index_cur])){
> +        printf(log);
> +      }
> +      else if (isnan(cpu_data[index_cur]) && isnan(gpu_data[index_cur])){
> +        printf(log);
> +      }
> +      else if( diff / cpu_data[index_cur] < 3 * FLT_ULP \
> +        && ( gpu_data[index_cur] > FLT_ULP  || cpu_data[index_cur] >
> FLT_ULP )){
> +        printf(log);
> +      }
> +      else if ( gpu_data[index_cur] < FLT_ULP && gpu_data[index_cur] <
> FLT_ULP)
> +        printf(log);
> +      else
> +        printf_c(log);
> +#else
> +      if (isinf(cpu_data[index_cur]))
> +        OCL_ASSERTM(isinf(gpu_data[index_cur]), log);
> +      else if (isnan(cpu_data[index_cur]))
> +        OCL_ASSERTM(isnan(gpu_data[index_cur]), log);
> +      else if ( gpu_data[index_cur] > FLT_ULP || cpu_data[index_cur] >
> FLT_ULP)
> +        OCL_ASSERTM(fabs( diff / cpu_data[index_cur]) < 3 * FLT_ULP, log);
> +      else
> +        OCL_ASSERTM(fabs(diff) < 3 * FLT_ULP, log);
> +#endif
> +    }
> +  }
> +}
> +
> +MAKE_UTEST_FROM_FUNCTION(builtin_exp)
> --
> 1.7.6.4
> 
> _______________________________________________
> Beignet mailing list
> Beignet at lists.freedesktop.org
> http://lists.freedesktop.org/mailman/listinfo/beignet