1 files changed, 46 insertions, 0 deletions
diff --git a/benchmarks/CUDA/WP/util4.cu b/benchmarks/CUDA/WP/util4.cu
new file mode 100644
index 0000000..4b629d4
--- /dev/null
+++ b/benchmarks/CUDA/WP/util4.cu
@@ -0,0 +1,46 @@
+#ifndef PREPASS
+#include <stdio.h>
+#include <stdlib.h>
+#include <math.h>
+#include "cublas.h"
+#endif
+
+__device__ float4 max4 ( const float4 a , const float4 b )
+{
+    float4 c  ;
+    c.x = (a.x>b.x)?a.x:b.x;
+    c.y = (a.y>b.y)?a.y:b.y;
+    c.z = (a.z>b.z)?a.z:b.z;
+    c.w = (a.w>b.w)?a.w:b.w;
+    return(c) ;
+}
+__device__ float4 min4 ( const float4 a , const float4 b )
+{
+    float4 c  ;
+    c.x = (a.x<b.x)?a.x:b.x;
+    c.y = (a.y<b.y)?a.y:b.y;
+    c.z = (a.z<b.z)?a.z:b.z;
+    c.w = (a.w<b.w)?a.w:b.w;
+    return(c) ;
+}
+
+__device__ float4 log4 ( const float4 a )
+{
+    float4 c  ;
+    c.x = log(a.x) ; c.y = log(a.y) ; c.z = log(a.z) ; c.w = log(a.w) ;
+    return(c) ;
+}
+
+__device__ float4 exp4 ( const float4 a )
+{
+    float4 c  ;
+    c.x = exp(a.x) ; c.y = exp(a.y) ; c.z = exp(a.z) ; c.w = exp(a.w) ;
+    return(c) ;
+}
+
+__device__ float4 sqrt4 ( const float4 a )
+{
+    float4 c  ;
+    c.x = sqrt(a.x) ; c.y = sqrt(a.y) ; c.z = sqrt(a.z) ; c.w = sqrt(a.w) ;
+    return(c) ;
+}