summaryrefslogtreecommitdiff
path: root/external/clBLAS/src/client/clfunc_xher2.hpp
diff options
context:
space:
mode:
Diffstat (limited to 'external/clBLAS/src/client/clfunc_xher2.hpp')
-rw-r--r--external/clBLAS/src/client/clfunc_xher2.hpp335
1 files changed, 0 insertions, 335 deletions
diff --git a/external/clBLAS/src/client/clfunc_xher2.hpp b/external/clBLAS/src/client/clfunc_xher2.hpp
deleted file mode 100644
index aec7cc83..00000000
--- a/external/clBLAS/src/client/clfunc_xher2.hpp
+++ /dev/null
@@ -1,335 +0,0 @@
-/* ************************************************************************
- * Copyright 2013 Advanced Micro Devices, Inc.
- *
- * Licensed under the Apache License, Version 2.0 (the "License");
- * you may not use this file except in compliance with the License.
- * You may obtain a copy of the License at
- *
- * http://www.apache.org/licenses/LICENSE-2.0
- *
- * Unless required by applicable law or agreed to in writing, software
- * distributed under the License is distributed on an "AS IS" BASIS,
- * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
- * See the License for the specific language governing permissions and
- * limitations under the License.
- * ************************************************************************/
-
-
-// $Id
-
-#ifndef CLBLAS_BENCHMARK_XHER2_HXX__
-#define CLBLAS_BENCHMARK_XHER2_HXX__
-
-#include "clfunc_common.hpp"
-
-template <typename T>
-struct xHer2Buffer
-{
- clblasOrder order;
- clblasUplo uplo;
- size_t N;
- T alpha;
- T* cpuX;
- cl_mem X;
- size_t offx;
- int incx;
- T* cpuY;
- cl_mem Y;
- size_t offy;
- int incy;
- T* cpuA;
- cl_mem A;
- size_t offa;
- size_t lda;
-}; // struct buffer
-
-template <typename T>
-class xHer2 : public clblasFunc
-{
-public:
- xHer2(StatisticalTimer& timer, cl_device_type devType) : clblasFunc(timer, devType)
- {
- timer.getUniqueID("clHer2", 0);
- }
-
- ~xHer2()
- {
- delete buffer.cpuA;
- delete buffer.cpuX;
- OPENCL_V_THROW( clReleaseMemObject(buffer.A), "releasing buffer A");
- OPENCL_V_THROW( clReleaseMemObject(buffer.X), "releasing buffer C");
- }
-
- double gflops()
- {
- return static_cast<double>((2 * buffer.N * buffer.N)/time_in_ns());
- }
-
- std::string gflops_formula()
- {
- return "2*N*N/time";
- }
-
- void setup_buffer(int order_option, int side_option, int
- uplo_option, int diag_option, int transA_option, int
- transB_option, size_t M, size_t N, size_t K,
- size_t lda, size_t ldb, size_t ldc,size_t offA,
- size_t offB, size_t offC, double alpha,
- double beta);
- void initialize_cpu_buffer();
- void initialize_gpu_buffer();
- void reset_gpu_write_buffer();
- void call_func();
- void read_gpu_buffer()
- {
- //cl_int err;
- //to-do need to fill up
- }
- void roundtrip_func()
- {//to-do need to fill up
- }
- void roundtrip_setup_buffer(int order_option, int side_option, int uplo_option,
- int diag_option, int transA_option, int transB_option,
- size_t M, size_t N, size_t K, size_t lda, size_t ldb,
- size_t ldc, size_t offA, size_t offBX, size_t offCY,
- double alpha, double beta)
- {}
- void releaseGPUBuffer_deleteCPUBuffer()
- {
- //this is necessary since we are running a iteration of tests and calculate the average time. (in client.cpp)
- //need to do this before we eventually hit the destructor
- //to do
- }
-protected:
-protected:
- void initialize_scalars(double alpha, double beta)
- {
- buffer.alpha = makeScalar<T>(alpha);
- //buffer.beta = makeScalar<T>(beta);
- }
-
-private:
- xHer2Buffer<T> buffer;
-};
-
-template <typename T>
-void xHer2<T>::setup_buffer(int order_option, int side_option, int
- uplo_option, int diag_option, int transA_option, int
- transB_option, size_t M, size_t N, size_t K,
- size_t lda, size_t ldb, size_t ldc,size_t offA,
- size_t offB, size_t offC, double alpha,
- double beta)
-{
- initialize_scalars(alpha, beta);
- buffer.offa = offA;
- buffer.offx = offB;
- buffer.incx = 1;//If this changes, remember to adjust size of Y in rest of the file
- buffer.offy = offC;
- buffer.incy = 1;//If this changes, remember to adjust size of Y in rest of the file
- buffer.N = M;
- if (order_option == 0)
- {
- buffer.order = clblasRowMajor;
- }
- else
- {
- buffer.order = clblasColumnMajor;
- }
- if (uplo_option == 0)
- {
- buffer.uplo = clblasUpper;
- }
- else
- {
- buffer.uplo = clblasLower;
- }
- if (lda == 0)
- {
- buffer.lda = buffer.N;
- }
- else if (lda < buffer.N)
- {
- std::cerr << "lda:wrong size\n";
- exit(1);
- }
- else
- {
- buffer.lda = lda;
- }
- buffer.cpuX = new T[buffer.N];
- buffer.cpuY = new T[buffer.N];
- buffer.cpuA = new T[buffer.N * buffer.lda];
- cl_int err;
- buffer.A = clCreateBuffer(ctx_, CL_MEM_READ_ONLY,
- buffer.N * buffer.lda*sizeof(T),
- NULL, &err);
-
- buffer.X = clCreateBuffer(ctx_, CL_MEM_READ_WRITE,
- buffer.N*sizeof(T),
- NULL, &err);
- buffer.Y = clCreateBuffer(ctx_, CL_MEM_READ_WRITE,
- buffer.N*sizeof(T),
- NULL, &err);
-}
-
-template <typename T>
-void xHer2<T>::initialize_cpu_buffer()
-{
- srand(10);
- for (size_t i = 0; i < buffer.N; ++i)
- {
- for (size_t j = 0; j < buffer.lda; ++j)
- {
- buffer.cpuA[i*buffer.lda+j] = random<T>(UPPER_BOUND<T>()) /
- randomScale<T>();
- }
- }
-
- for (size_t i = 0; i < buffer.N; ++i)
- {
- buffer.cpuX[i] = random<T>(UPPER_BOUND<T>()) /
- randomScale<T>();
- buffer.cpuY[i] = random<T>(UPPER_BOUND<T>()) /
- randomScale<T>();
- }
-}
-
-//template <>
-//void xHer2<cl_float2>::initialize_cpu_buffer()
-//{
-// srand(10);
-// for (size_t i = 0; i < buffer.N; ++i)
-// {
-// for (size_t j = 0; j < buffer.lda; ++j)
-// {
-// buffer.cpuA[i*buffer.lda+j].s[0] = static_cast<cl_float>(rand())/
-// static_cast<cl_float>(RAND_MAX);
-// buffer.cpuA[i*buffer.lda+j].s[1] = static_cast<cl_float>(rand())/
-// static_cast<cl_float>(RAND_MAX);
-// }
-// }
-//
-// for (size_t i = 0; i < buffer.N; ++i)
-// {
-// buffer.cpuX[i].s[0] = static_cast<cl_float>(rand())/
-// static_cast<cl_float>(RAND_MAX);
-// buffer.cpuX[i].s[1] = static_cast<cl_float>(rand())/
-// static_cast<cl_float>(RAND_MAX);
-// buffer.cpuY[i].s[0] = static_cast<cl_float>(rand())/
-// static_cast<cl_float>(RAND_MAX);
-// buffer.cpuY[i].s[1] = static_cast<cl_float>(rand())/
-// static_cast<cl_float>(RAND_MAX);
-// }
-//}
-//template <>
-//void xHer2<cl_double2>::initialize_cpu_buffer()
-//{
-// srand(10);
-// for (size_t i = 0; i < buffer.N; ++i)
-// {
-// for (size_t j = 0; j < buffer.lda; ++j)
-// {
-// buffer.cpuA[i*buffer.lda+j].s[0] = static_cast<cl_double>(rand())/
-// static_cast<cl_double>(RAND_MAX);
-// buffer.cpuA[i*buffer.lda+j].s[1] = static_cast<cl_double>(rand())/
-// static_cast<cl_double>(RAND_MAX);
-// }
-// }
-//
-// for (size_t i = 0; i < buffer.N; ++i)
-// {
-// buffer.cpuX[i].s[0] = static_cast<cl_double>(rand())/
-// static_cast<cl_double>(RAND_MAX);
-// buffer.cpuX[i].s[1] = static_cast<cl_double>(rand())/
-// static_cast<cl_double>(RAND_MAX);
-// buffer.cpuY[i].s[0] = static_cast<cl_double>(rand())/
-// static_cast<cl_double>(RAND_MAX);
-// buffer.cpuY[i].s[1] = static_cast<cl_double>(rand())/
-// static_cast<cl_double>(RAND_MAX);
-// }
-//}
-//
-
-template <typename T>
-void xHer2<T>::initialize_gpu_buffer()
-{
- cl_int err;
-
- err = clEnqueueWriteBuffer(queue_, buffer.A, CL_TRUE,
- buffer.offa * sizeof(T),
- buffer.N * buffer.lda*sizeof(T),
- buffer.cpuA, 0, NULL, NULL);
-
- err = clEnqueueWriteBuffer(queue_, buffer.X, CL_TRUE, 0,
- buffer.N*sizeof(T),
- buffer.cpuX, 0, NULL, NULL);
- err = clEnqueueWriteBuffer(queue_, buffer.Y, CL_TRUE, 0,
- buffer.N*sizeof(T),
- buffer.cpuY, 0, NULL, NULL);
-}
-
-template <typename T>
-void xHer2<T>::reset_gpu_write_buffer()
-{
- cl_int err;
- err = clEnqueueWriteBuffer(queue_, buffer.A, CL_TRUE,
- buffer.offa * sizeof(T),
- buffer.N * buffer.lda*sizeof(T),
- buffer.cpuA, 0, NULL, NULL);;
-}
-
-template <>
-void xHer2<cl_float2>::call_func()
-{
- timer.Start(timer_id);
- clblasCher2(buffer.order, buffer.uplo, buffer.N, buffer.alpha, buffer.X, buffer.offx,
- buffer.incx, buffer.Y, buffer.offy, buffer.incy, buffer.A, buffer.offa,
- buffer.lda, 1, &queue_, 0, NULL,&event_);
- clWaitForEvents(1, &event_);
- timer.Stop(timer_id);
-}
-
-template <>
-void xHer2<cl_double2>::call_func()
-{
- timer.Start(timer_id);
- clblasZher2(buffer.order, buffer.uplo, buffer.N, buffer.alpha, buffer.X, buffer.offx,
- buffer.incx, buffer.Y, buffer.offy, buffer.incy, buffer.A, buffer.offa,
- buffer.lda, 1, &queue_, 0, NULL,&event_);
- clWaitForEvents(1, &event_);
- timer.Stop(timer_id);
-}
-
-template<>
-double
-xHer2<cl_float2>::
-gflops()
-{
- return static_cast<double>((8 * buffer.N * buffer.N)/time_in_ns());
-}
-
-template<>
-double
-xHer2<cl_double2>::
-gflops()
-{
- return static_cast<double>((8 * buffer.N * buffer.N)/time_in_ns());
-}
-
-template<>
-std::string
-xHer2<cl_float2>::
-gflops_formula()
-{
- return "8*N*N/time";
-}
-
-template<>
-std::string
-xHer2<cl_double2>::
-gflops_formula()
-{
- return "8*N*N/time";
-}
-
-#endif // ifndef CLBLAS_BENCHMARK_XSYR_HXX__ \ No newline at end of file