[Beignet] [PATCH] utest: add test case for builtin function exp/exp2/exp10/expm1.
Shui, YangweiX
yangweix.shui at intel.com
Fri Oct 25 05:27:32 CEST 2013
Function printf_c() is not necessary, it would help us to highlight some error messages, if it troubles you, we will modify it to printf when pushed.
> -----Original Message-----
> From: Lu, Guanqun
> Sent: Friday, October 25, 2013 9:07 AM
> To: Sun, Yi; beignet at lists.freedesktop.org
> Cc: Sun, Yi; Shui, YangweiX
> Subject: RE: [Beignet] [PATCH] utest: add test case for builtin function
> exp/exp2/exp10/expm1.
>
> 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
More information about the Beignet
mailing list