From e36182cd8c5dd5c6d0aa1f77cf5cdca87e8bb1f0 Mon Sep 17 00:00:00 2001
From: Joseph Redmon <pjreddie@gmail.com>
Date: Fri, 21 Nov 2014 23:35:19 +0000
Subject: [PATCH] cleaned up data parsing a lot. probably nothing broken?
---
src/axpy.c | 123 +++++++++++++++++++++++++++++++++++++++++
1 files changed, 123 insertions(+), 0 deletions(-)
diff --git a/src/axpy.c b/src/axpy.c
index 750f47e..21293b3 100644
--- a/src/axpy.c
+++ b/src/axpy.c
@@ -12,3 +12,126 @@
for(i = 0; i < N; ++i) X[i*INCX] *= ALPHA;
}
+void copy_cpu(int N, float *X, int INCX, float *Y, int INCY)
+{
+ int i;
+ for(i = 0; i < N; ++i) Y[i*INCY] = X[i*INCX];
+}
+
+float dot_cpu(int N, float *X, int INCX, float *Y, int INCY)
+{
+ int i;
+ float dot = 0;
+ for(i = 0; i < N; ++i) dot += X[i*INCX] * Y[i*INCY];
+ return dot;
+}
+
+#ifdef GPU
+#include "opencl.h"
+
+cl_kernel get_axpy_kernel()
+{
+ static int init = 0;
+ static cl_kernel kernel;
+ if(!init){
+ kernel = get_kernel("src/axpy.cl", "axpy", 0);
+ init = 1;
+ }
+ return kernel;
+}
+
+cl_kernel get_copy_kernel()
+{
+ static int init = 0;
+ static cl_kernel kernel;
+ if(!init){
+ kernel = get_kernel("src/axpy.cl", "copy", 0);
+ init = 1;
+ }
+ return kernel;
+}
+
+cl_kernel get_scal_kernel()
+{
+ static int init = 0;
+ static cl_kernel kernel;
+ if(!init){
+ kernel = get_kernel("src/axpy.cl", "scal", 0);
+ init = 1;
+ }
+ return kernel;
+}
+
+
+void axpy_ongpu(int N, float ALPHA, cl_mem X, int INCX, cl_mem Y, int INCY)
+{
+ axpy_ongpu_offset(N,ALPHA,X,0,INCX,Y,0,INCY);
+}
+
+void axpy_ongpu_offset(int N, float ALPHA, cl_mem X, int OFFX, int INCX, cl_mem Y, int OFFY, int INCY)
+{
+ cl_setup();
+ cl_kernel kernel = get_axpy_kernel();
+ cl_command_queue queue = cl.queue;
+
+ cl_uint i = 0;
+ cl.error = clSetKernelArg(kernel, i++, sizeof(N), (void*) &N);
+ cl.error = clSetKernelArg(kernel, i++, sizeof(ALPHA), (void*) &ALPHA);
+ cl.error = clSetKernelArg(kernel, i++, sizeof(X), (void*) &X);
+ cl.error = clSetKernelArg(kernel, i++, sizeof(OFFX), (void*) &OFFX);
+ cl.error = clSetKernelArg(kernel, i++, sizeof(INCX), (void*) &INCX);
+ cl.error = clSetKernelArg(kernel, i++, sizeof(Y), (void*) &Y);
+ cl.error = clSetKernelArg(kernel, i++, sizeof(OFFY), (void*) &OFFY);
+ cl.error = clSetKernelArg(kernel, i++, sizeof(INCY), (void*) &INCY);
+ check_error(cl);
+
+ const size_t global_size[] = {N};
+
+ cl.error = clEnqueueNDRangeKernel(queue, kernel, 1, 0, global_size, 0, 0, 0, 0);
+ check_error(cl);
+
+}
+void copy_ongpu(int N, cl_mem X, int INCX, cl_mem Y, int INCY)
+{
+ copy_ongpu_offset(N,X,0,INCX,Y,0,INCY);
+}
+void copy_ongpu_offset(int N, cl_mem X, int OFFX, int INCX, cl_mem Y, int OFFY, int INCY)
+{
+ cl_setup();
+ cl_kernel kernel = get_copy_kernel();
+ cl_command_queue queue = cl.queue;
+
+ cl_uint i = 0;
+ cl.error = clSetKernelArg(kernel, i++, sizeof(N), (void*) &N);
+ cl.error = clSetKernelArg(kernel, i++, sizeof(X), (void*) &X);
+ cl.error = clSetKernelArg(kernel, i++, sizeof(OFFX), (void*) &OFFX);
+ cl.error = clSetKernelArg(kernel, i++, sizeof(INCX), (void*) &INCX);
+ cl.error = clSetKernelArg(kernel, i++, sizeof(Y), (void*) &Y);
+ cl.error = clSetKernelArg(kernel, i++, sizeof(OFFY), (void*) &OFFY);
+ cl.error = clSetKernelArg(kernel, i++, sizeof(INCY), (void*) &INCY);
+ check_error(cl);
+
+ const size_t global_size[] = {N};
+
+ cl.error = clEnqueueNDRangeKernel(queue, kernel, 1, 0, global_size, 0, 0, 0, 0);
+ check_error(cl);
+}
+void scal_ongpu(int N, float ALPHA, cl_mem X, int INCX)
+{
+ cl_setup();
+ cl_kernel kernel = get_scal_kernel();
+ cl_command_queue queue = cl.queue;
+
+ cl_uint i = 0;
+ cl.error = clSetKernelArg(kernel, i++, sizeof(N), (void*) &N);
+ cl.error = clSetKernelArg(kernel, i++, sizeof(ALPHA), (void*) &ALPHA);
+ cl.error = clSetKernelArg(kernel, i++, sizeof(X), (void*) &X);
+ cl.error = clSetKernelArg(kernel, i++, sizeof(INCX), (void*) &INCX);
+ check_error(cl);
+
+ const size_t global_size[] = {N};
+
+ cl.error = clEnqueueNDRangeKernel(queue, kernel, 1, 0, global_size, 0, 0, 0, 0);
+ check_error(cl);
+}
+#endif
--
Gitblit v1.10.0