diff options
Diffstat (limited to 'external/clBLAS/src/client/clfunc_xherk.hpp')
-rw-r--r-- | external/clBLAS/src/client/clfunc_xherk.hpp | 535 |
1 files changed, 0 insertions, 535 deletions
diff --git a/external/clBLAS/src/client/clfunc_xherk.hpp b/external/clBLAS/src/client/clfunc_xherk.hpp deleted file mode 100644 index 74871a39..00000000 --- a/external/clBLAS/src/client/clfunc_xherk.hpp +++ /dev/null @@ -1,535 +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_XHERK_HXX__ -#define CLBLAS_BENCHMARK_XHERK_HXX__ - -#include "clfunc_common.hpp" - -template <typename T> -struct xHerkBuffer -{ - clblasOrder order_; - clblasUplo uplo_; - clblasTranspose transA_; - size_t N_; - size_t K_; - T alpha_; - cl_mem A_; - size_t offa_; - size_t lda_; - T beta_; - cl_mem C_; - size_t offc_; - size_t ldc_; - size_t a_num_vectors_; - size_t c_num_vectors_; - T* cpuA_; - T* cpuC_; -}; // struct buffer - -template <typename T> -class xHerk : public clblasFunc -{ -public: - xHerk(StatisticalTimer& timer, cl_device_type devType) : clblasFunc(timer, devType) - { - timer.getUniqueID("clHerk", 0); - } - - ~xHerk() - { - } - - double gflops() - { - return static_cast<double>(4*(buffer_.K_ * buffer_.N_ * (buffer_.N_+1))/time_in_ns()); - } - - std::string gflops_formula() - { - return "4*K*N*(N+1)/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) - { - DUMMY_ARGS_USAGE_4(side_option, diag_option, transB_option, M); - DUMMY_ARGS_USAGE_2(ldb, offB); - - initialize_scalars(alpha,beta); - - buffer_.N_ = N; - buffer_.K_ = K; - buffer_.offa_ = offA; - buffer_.offc_ = offC; - - if (uplo_option == 0) - { - buffer_.uplo_ = clblasUpper; - } - else - { - buffer_.uplo_ = clblasLower; - } - - if (ldc == 0) - { - buffer_.ldc_ = N; - } - else if (ldc < N) - { - std::cerr << "ldc:wrong size\n"; - } - else - { - buffer_.ldc_ = ldc; - } - - buffer_.c_num_vectors_ = N; - - if (order_option == 0) - { - order_ = clblasRowMajor; - if (transA_option == 0) - { - buffer_.transA_ = clblasNoTrans; - buffer_.a_num_vectors_ = N; - if (lda == 0) - { - buffer_.lda_ = K; - } - else if (lda < K) - { - std::cerr << "lda:wrong size\n"; - exit(1); - } - else - { - buffer_.lda_ = lda; - } - } - else - { - buffer_.a_num_vectors_ = K; - if (transA_option == 1) - { - buffer_.transA_ = clblasTrans; - } - else if (transA_option == 2) - { - buffer_.transA_ = clblasConjTrans; - } - if (lda == 0) - { - buffer_.lda_ = N; - } - else if (lda < N) - { - std::cerr << "lda:wrong size\n"; - exit(1); - } - else - { - buffer_.lda_ = lda; - } - } - } - else - { - order_ = clblasColumnMajor; - if (transA_option == 0) - { - buffer_.a_num_vectors_ = K; - buffer_.transA_ = clblasNoTrans; - if (lda == 0) - { - buffer_.lda_ = N; - } - else if (lda < N) - { - std::cerr << "lda:wrong size\n"; - exit(1); - } - else - { - buffer_.lda_ = lda; - } - } - else - { - buffer_.a_num_vectors_ = N; - if (transA_option == 1) - { - buffer_.transA_ = clblasTrans; - } - else if (transA_option == 2) - { - buffer_.transA_ = clblasConjTrans; - } - - if (lda == 0) - { - buffer_.lda_ = K; - } - else if (lda < K) - { - std::cerr << "lda:wrong size\n"; - exit(1); - } - else - { - buffer_.lda_ = lda; - } - } - } - - buffer_.cpuA_ = new T[buffer_.lda_*buffer_.a_num_vectors_]; - buffer_.cpuC_ = new T[buffer_.ldc_*buffer_.c_num_vectors_]; - - cl_int err; - buffer_.A_ = clCreateBuffer(ctx_, CL_MEM_READ_ONLY, - (buffer_.lda_ * buffer_.a_num_vectors_ + - buffer_.offa_) * sizeof(T), - NULL, &err); - - buffer_.C_ = clCreateBuffer(ctx_, CL_MEM_READ_WRITE, - (buffer_.ldc_ * buffer_.c_num_vectors_ + - buffer_.offc_) * sizeof(T), - NULL, &err); - } - void initialize_cpu_buffer() - { - srand(10); - for (size_t i = 0; i < buffer_.a_num_vectors_; ++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) - { - for (size_t j = 0; j < buffer_.ldc_; ++j) - { - buffer_.cpuC_[i*buffer_.ldc_+j] = random<T>(UPPER_BOUND<T>()) / - randomScale<T>(); - } - } - } - void initialize_gpu_buffer() - { - cl_int err; - - err = clEnqueueWriteBuffer(queue_, buffer_.A_, CL_TRUE, - buffer_.offa_ * sizeof(T), - buffer_.lda_ * buffer_.a_num_vectors_ * - sizeof(T), - buffer_.cpuA_, 0, NULL, NULL); - - err = clEnqueueWriteBuffer(queue_, buffer_.C_, CL_TRUE, - buffer_.offa_ * sizeof(T), - buffer_.ldc_ * buffer_.c_num_vectors_ * - sizeof(T), - buffer_.cpuC_, 0, NULL, NULL); - } - void reset_gpu_write_buffer() - { - cl_int err; - - err = clEnqueueWriteBuffer(queue_, buffer_.C_, CL_TRUE, - buffer_.offc_ * sizeof(T), - buffer_.ldc_ * buffer_.c_num_vectors_ * - sizeof(T), - buffer_.cpuC_, 0, NULL, NULL); - } - void call_func(); - void read_gpu_buffer() - { - cl_int err; - err = clEnqueueReadBuffer(queue_, buffer_.C_, CL_TRUE, - buffer_.offc_*sizeof(T), buffer_.ldc_*buffer_.c_num_vectors_*sizeof(T), - buffer_.cpuC_, 0, NULL, NULL); - } - void roundtrip_func(); - void zerocopy_roundtrip_func() - { - std::cout << "xTrmm::zerocopy_roundtrip_func\n"; - } - 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) - { - DUMMY_ARGS_USAGE_4(side_option, diag_option, transB_option, M); - DUMMY_ARGS_USAGE_2(ldb, offBX); - - initialize_scalars(alpha,beta); - - buffer_.N_ = N; - buffer_.K_ = K; - buffer_.offa_ = offA; - buffer_.offc_ = offCY; - - if (uplo_option == 0) - { - buffer_.uplo_ = clblasUpper; - } - else - { - buffer_.uplo_ = clblasLower; - } - - if (ldc == 0) - { - buffer_.ldc_ = N; - } - else if (ldc < N) - { - std::cerr << "ldc:wrong size\n"; - } - else - { - buffer_.ldc_ = ldc; - } - - buffer_.c_num_vectors_ = N; - - if (order_option == 0) - { - order_ = clblasRowMajor; - if (transA_option == 0) - { - buffer_.transA_ = clblasNoTrans; - buffer_.a_num_vectors_ = N; - if (lda == 0) - { - buffer_.lda_ = K; - } - else if (lda < K) - { - std::cerr << "lda:wrong size\n"; - exit(1); - } - else - { - buffer_.lda_ = lda; - } - } - else - { - buffer_.a_num_vectors_ = K; - if (transA_option == 1) - { - buffer_.transA_ = clblasTrans; - } - else if (transA_option == 2) - { - buffer_.transA_ = clblasConjTrans; - } - if (lda == 0) - { - buffer_.lda_ = N; - } - else if (lda < N) - { - std::cerr << "lda:wrong size\n"; - exit(1); - } - else - { - buffer_.lda_ = lda; - } - } - } - else - { - order_ = clblasColumnMajor; - if (transA_option == 0) - { - buffer_.a_num_vectors_ = K; - buffer_.transA_ = clblasNoTrans; - if (lda == 0) - { - buffer_.lda_ = N; - } - else if (lda < N) - { - std::cerr << "lda:wrong size\n"; - exit(1); - } - else - { - buffer_.lda_ = lda; - } - } - else - { - buffer_.a_num_vectors_ = N; - if (transA_option == 1) - { - buffer_.transA_ = clblasTrans; - } - else if (transA_option == 2) - { - buffer_.transA_ = clblasConjTrans; - } - - if (lda == 0) - { - buffer_.lda_ = K; - } - else if (lda < K) - { - std::cerr << "lda:wrong size\n"; - exit(1); - } - else - { - buffer_.lda_ = lda; - } - } - } - - buffer_.cpuA_ = new T[buffer_.lda_*buffer_.a_num_vectors_]; - buffer_.cpuC_ = new T[buffer_.ldc_*buffer_.c_num_vectors_]; - } - 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 - delete buffer_.cpuA_; - delete buffer_.cpuC_; - OPENCL_V_THROW( clReleaseMemObject(buffer_.A_), "releasing buffer A"); - OPENCL_V_THROW( clReleaseMemObject(buffer_.C_), "releasing buffer C"); - } -protected: -protected: - void initialize_scalars(double alpha, double beta) - { - buffer_.alpha_ = makeScalar<T>(alpha); - buffer_.beta_ = makeScalar<T>(beta); - } - -private: - xHerkBuffer<T> buffer_; -}; - -template<> -void -xHerk<cl_float2>::call_func() -{ - timer.Start(timer_id); - - clblasCherk(order_, buffer_.uplo_, buffer_.transA_, - buffer_.N_, buffer_.K_, buffer_.alpha_.s[0], - buffer_.A_, buffer_.offa_, buffer_.lda_, - buffer_.beta_.s[0], buffer_.C_, buffer_.offc_, - buffer_.ldc_, 1, &queue_, 0, NULL, &event_); - - clWaitForEvents(1, &event_); - timer.Stop(timer_id); -} - -template<> -void -xHerk<cl_float2>::roundtrip_func() -{ - timer.Start(timer_id); - cl_int err; - buffer_.A_ = clCreateBuffer(ctx_, CL_MEM_READ_ONLY, - (buffer_.lda_ * buffer_.a_num_vectors_ + - buffer_.offa_) * sizeof(cl_float2), - NULL, &err); - - buffer_.C_ = clCreateBuffer(ctx_, CL_MEM_READ_WRITE, - (buffer_.ldc_ * buffer_.c_num_vectors_ + - buffer_.offc_) * sizeof(cl_float2), - NULL, &err); - this->initialize_gpu_buffer(); - - clblasCherk(order_, buffer_.uplo_, buffer_.transA_, - buffer_.N_, buffer_.K_, buffer_.alpha_.s[0], - buffer_.A_, buffer_.offa_, buffer_.lda_, - buffer_.beta_.s[0], buffer_.C_, buffer_.offc_, - buffer_.ldc_, 1, &queue_, 0, NULL, NULL); - - err = clEnqueueWriteBuffer(queue_, buffer_.C_, CL_TRUE, - buffer_.offc_ * sizeof(cl_float2), - buffer_.ldc_ * buffer_.c_num_vectors_ * - sizeof(cl_float2), - buffer_.cpuC_, 0, NULL, &event_); - clWaitForEvents(1, &event_); - timer.Stop(timer_id); -} - -template<> -void -xHerk<cl_double2>::call_func() -{ - timer.Start(timer_id); - - clblasZherk(order_, buffer_.uplo_, buffer_.transA_, - buffer_.N_, buffer_.K_, buffer_.alpha_.s[0], - buffer_.A_, buffer_.offa_, buffer_.lda_, - buffer_.beta_.s[0], buffer_.C_, buffer_.offc_, - buffer_.ldc_, 1, &queue_, 0, NULL, &event_); - - clWaitForEvents(1, &event_); - timer.Stop(timer_id); -} - -template<> -void -xHerk<cl_double2>::roundtrip_func() -{ - timer.Start(timer_id); - cl_int err; - buffer_.A_ = clCreateBuffer(ctx_, CL_MEM_READ_ONLY, - (buffer_.lda_ * buffer_.a_num_vectors_ + - buffer_.offa_) * sizeof(cl_double2), - NULL, &err); - - buffer_.C_ = clCreateBuffer(ctx_, CL_MEM_READ_WRITE, - (buffer_.ldc_ * buffer_.c_num_vectors_ + - buffer_.offc_) * sizeof(cl_double2), - NULL, &err); - this->initialize_gpu_buffer(); - - clblasZherk(order_, buffer_.uplo_, buffer_.transA_, - buffer_.N_, buffer_.K_, buffer_.alpha_.s[0], - buffer_.A_, buffer_.offa_, buffer_.lda_, - buffer_.beta_.s[0], buffer_.C_, buffer_.offc_, - buffer_.ldc_, 1, &queue_, 0, NULL, NULL); - - err = clEnqueueWriteBuffer(queue_, buffer_.C_, CL_TRUE, - buffer_.offc_ * sizeof(cl_double2), - buffer_.ldc_ * buffer_.c_num_vectors_ * - sizeof(cl_double2), - buffer_.cpuC_, 0, NULL, &event_); - clWaitForEvents(1, &event_); - timer.Stop(timer_id); -} -#endif // ifndef CLBLAS_BENCHMARK_XSYR_HXX__
\ No newline at end of file |