utests: add utest for sqrt-div optimization

Signed-off-by: rander.wang <rander.wang@intel.com>
Reviewed-by: Yang Rong <rong.r.yang@intel.com>

1 files changed, 61 insertions, 0 deletions

diff --git a/utests/compiler_sqrt_div.cpp b/utests/compiler_sqrt_div.cpp
new file mode 100644
index 00000000..c22c5a39
--- /dev/null
+++ b/utests/compiler_sqrt_div.cpp
@@ -0,0 +1,61 @@
+#include "utest_helper.hpp"
+#include <cmath>
+
+void compiler_sqrt_div(void) {
+  const int n = 1024;
+  float src[n];
+
+  // Setup kernel and buffers
+  OCL_CREATE_KERNEL("compiler_sqrt_div");
+  OCL_CREATE_BUFFER(buf[0], 0, n * sizeof(float), NULL);
+  OCL_CREATE_BUFFER(buf[1], 0, n * 4 * sizeof(float), NULL);
+  OCL_SET_ARG(0, sizeof(cl_mem), &buf[0]);
+  OCL_SET_ARG(1, sizeof(cl_mem), &buf[1]);
+  globals[0] = n;
+  locals[0] = 16;
+
+  for (int j = 0; j < 1024; j++) {
+    OCL_MAP_BUFFER(0);
+    for (int i = 0; i < n; ++i) {
+      src[i] = ((float *)buf_data[0])[i] = (j * n + i + 1) * 0.001f;
+    }
+    OCL_UNMAP_BUFFER(0);
+
+    OCL_NDRANGE(1);
+
+    OCL_MAP_BUFFER(1);
+    float *dst = (float *)buf_data[1];
+    for (int i = 0; i < n; ++i) {
+      float cpu = 1.0f / sqrt(src[i]);
+      float gpu = dst[4 * i];
+      if (fabsf(cpu - gpu) >= 1e-3) {
+        printf("%f %f %f", src[i], cpu, gpu);
+        OCL_ASSERT(0);
+      }
+
+      cpu = i / sqrt(src[i]);
+      gpu = dst[4 * i + 1];
+      if (fabsf(cpu - gpu) >= 1e-3) {
+        printf("%f %f %f", src[i], cpu, gpu);
+        OCL_ASSERT(0);
+      }
+
+      cpu = 2.0f / sqrt(src[i]);
+      gpu = dst[4 * i + 2];
+      if (fabsf(cpu - gpu) >= 1e-3) {
+        printf("%f %f %f", src[i], cpu, gpu);
+        OCL_ASSERT(0);
+      }
+
+      cpu = 1.0f / sqrt(src[i]) + sqrt(src[i]);
+      gpu = dst[4 * i + 3];
+      if (fabsf(cpu - gpu) >= 1e-3) {
+        printf("%f %f %f", src[i], cpu, gpu);
+        OCL_ASSERT(0);
+      }
+    }
+    OCL_UNMAP_BUFFER(1);
+  }
+}
+
+MAKE_UTEST_FROM_FUNCTION(compiler_sqrt_div);