summaryrefslogtreecommitdiff
path: root/external/clBLAS/src/tests/performance/perf-hpr.cpp
diff options
context:
space:
mode:
Diffstat (limited to 'external/clBLAS/src/tests/performance/perf-hpr.cpp')
-rw-r--r--external/clBLAS/src/tests/performance/perf-hpr.cpp319
1 files changed, 0 insertions, 319 deletions
diff --git a/external/clBLAS/src/tests/performance/perf-hpr.cpp b/external/clBLAS/src/tests/performance/perf-hpr.cpp
deleted file mode 100644
index ea990e23..00000000
--- a/external/clBLAS/src/tests/performance/perf-hpr.cpp
+++ /dev/null
@@ -1,319 +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.
- * ************************************************************************/
-
-#include <stdlib.h> // srand()
-#include <string.h> // memcpy()
-#include <gtest/gtest.h>
-#include <clBLAS.h>
-
-#include <common.h>
-#include <clBLAS-wrapper.h>
-#include <BlasBase.h>
-#include <hpr.h>
-#include <blas-random.h>
-
-#ifdef PERF_TEST_WITH_ACML
-#include <blas-internal.h>
-#include <blas-wrapper.h>
-#endif
-
-#include "PerformanceTest.h"
-
-/*
- * NOTE: operation factor means overall number
- * of multiply and add per each operation involving
- * 2 matrix elements
- */
-
-using namespace std;
-using namespace clMath;
-
-#define CHECK_RESULT(ret) \
-do { \
- ASSERT_GE(ret, 0) << "Fatal error: can not allocate resources or " \
- "perform an OpenCL request!" << endl; \
- EXPECT_EQ(0, ret) << "The OpenCL version is slower in the case" << \
- endl; \
-} while (0)
-
-namespace clMath {
-
-template <typename ElemType> class HprPerformanceTest : public PerformanceTest
-{
-public:
- virtual ~HprPerformanceTest();
-
- virtual int prepare(void);
- virtual nano_time_t etalonPerfSingle(void);
- virtual nano_time_t clblasPerfSingle(void);
-
- static void runInstance(BlasFunction fn, TestParams *params)
- {
- HprPerformanceTest<ElemType> perfCase(fn, params);
- int ret = 0;
- int opFactor;
- BlasBase *base;
-
- base = clMath::BlasBase::getInstance();
-
- opFactor = 1;
-
- if ((fn == FN_ZHPR) &&
- !base->isDevSupportDoublePrecision()) {
-
- std::cerr << ">> WARNING: The target device doesn't support native "
- "double precision floating point arithmetic" <<
- std::endl << ">> Test skipped" << std::endl;
- return;
- }
-
- if (!perfCase.areResourcesSufficient(params)) {
- std::cerr << ">> RESOURCE CHECK: Skip due to unsufficient resources" <<
- std::endl;
- }
- else {
- ret = perfCase.run(opFactor);
- }
-
- ASSERT_GE(ret, 0) << "Fatal error: can not allocate resources or "
- "perform an OpenCL request!" << endl;
- EXPECT_EQ(0, ret) << "The OpenCL version is slower in the case" << endl;
- }
-
-private:
- HprPerformanceTest(BlasFunction fn, TestParams *params);
-
- bool areResourcesSufficient(TestParams *params);
-
- TestParams params_;
- ElemType alpha_;
- ElemType *AP_;
- ElemType *X_;
- ElemType *backAP_;
- cl_mem mobjAP_;
- cl_mem mobjX_;
- ::clMath::BlasBase *base_;
-};
-
-template <typename ElemType>
-HprPerformanceTest<ElemType>::HprPerformanceTest(
- BlasFunction fn,
- TestParams *params) : PerformanceTest(fn, (problem_size_t)((((params->N * params->N) + params->N) * 2 ) * sizeof(ElemType))),
- params_(*params), mobjAP_(NULL), mobjX_(NULL)
-{
- AP_ = new ElemType[( ( params_.N*( params_.N + 1 ) )/2 )+ params_.offa];
- X_ = new ElemType[ 1 + (params_.N-1) * abs(params_.incx) + params_.offBX];
- backAP_ = new ElemType[( ( params_.N*( params_.N + 1 ) )/2 )+ params_.offa];
-
- base_ = ::clMath::BlasBase::getInstance();
-}
-
-template <typename ElemType>
-HprPerformanceTest<ElemType>::~HprPerformanceTest()
-{
- if(AP_ != NULL)
- {
- delete[] AP_;
- }
- if(backAP_ != NULL)
- {
- delete[] backAP_;
- }
- if(X_ != NULL)
- {
- delete[] X_;
- }
-
- if(mobjX_ != NULL) {
- clReleaseMemObject(mobjX_);
- }
- if(mobjAP_ != NULL) {
- clReleaseMemObject(mobjAP_);
- }
-}
-
-/*
- * Check if available OpenCL resources are sufficient to
- * run the test case
- */
-template <typename ElemType> bool
-HprPerformanceTest<ElemType>::areResourcesSufficient(TestParams *params)
-{
- clMath::BlasBase *base;
- size_t gmemSize, allocSize;
- size_t n = params->N;
-
- if((AP_ == NULL) || (backAP_ == NULL) || (X_ == NULL))
- {
- return 0;
- }
-
- base = clMath::BlasBase::getInstance();
- gmemSize = (size_t)base->availGlobalMemSize(0);
- allocSize = (size_t)base->maxMemAllocSize();
-
- bool suff = ( sizeof(ElemType) *(( n*( n + 1 ) )/2 )< allocSize ) && ((1 + (n-1)*abs(params->incx))*sizeof(ElemType) < allocSize); //for individual allocations
- suff = suff && (((( (n*( n + 1 ) )/2 ) + (1 + (n-1)*abs(params->incx))*2)*sizeof(ElemType)) < gmemSize) ; //for total global allocations
-
- return suff ;
-}
-
-template <typename ElemType> int
-HprPerformanceTest<ElemType>::prepare(void)
-{
- size_t lenX = 1 + (params_.N-1) * abs(params_.incx);
- alpha_ = convertMultiplier<ElemType>(params_.alpha);
- randomHerMatrices( params_.order, params_.uplo, params_.N, &alpha_, (AP_ + params_.offa), 0, (X_ + params_.offBX), params_.incx );
- memcpy(backAP_, AP_, ((( ( params_.N*( params_.N + 1 ) )/2 )+ params_.offa)* sizeof(ElemType)));
-
- mobjAP_ = base_->createEnqueueBuffer(AP_, (( ( params_.N*( params_.N + 1 ) )/2 )+ params_.offa)* sizeof(*AP_), 0, CL_MEM_READ_WRITE);
- mobjX_ = base_->createEnqueueBuffer(X_, (lenX + params_.offBX )* sizeof(*X_), 0, CL_MEM_READ_ONLY);
-
- return ( (mobjAP_ != NULL) && (mobjX_ != NULL) ) ? 0 : -1;
-}
-
-template <typename ElemType> nano_time_t
-HprPerformanceTest<ElemType>::etalonPerfSingle(void)
-{
- nano_time_t time = 0;
- clblasOrder order;
-// size_t lda;
-
-#ifndef PERF_TEST_WITH_ROW_MAJOR
- if (params_.order == clblasRowMajor) {
- cerr << "Row major order is not allowed" << endl;
- return NANOTIME_ERR;
- }
-#endif
-
- order = params_.order;
-
-
-#ifdef PERF_TEST_WITH_ACML
-
- clblasOrder fOrder;
- clblasUplo fUplo;
- fOrder = params_.order;
- fUplo = params_.uplo;
-
- if (order != clblasColumnMajor)
- {
- doConjugate( (X_ + params_.offBX), (1 + (params_.N-1) * abs(params_.incx)), 1, 1 );
- fOrder = clblasColumnMajor;
- fUplo = (fUplo == clblasLower)? clblasUpper : clblasLower;
- }
-
- time = getCurrentTime();
- clMath::blas::hpr(fOrder, fUplo, params_.N, CREAL(alpha_), X_, params_.offBX, params_.incx, AP_, params_.offa);
- time = getCurrentTime() - time;
-
-#endif // PERF_TEST_WITH_ACML
-
- return time;
-}
-
-
-template <typename ElemType> nano_time_t
-HprPerformanceTest<ElemType>::clblasPerfSingle(void)
-{
- nano_time_t time;
- cl_event event;
- cl_int status;
- cl_command_queue queue = base_->commandQueues()[0];
-
- status = clEnqueueWriteBuffer(queue, mobjAP_, CL_TRUE, 0,
- ((( params_.N*( params_.N + 1 ) )/2 ) + params_.offa) *
- sizeof(ElemType), backAP_, 0, NULL, &event);
- if (status != CL_SUCCESS) {
- cerr << "Matrix A buffer object enqueuing error, status = " <<
- status << endl;
-
- return NANOTIME_ERR;
- }
-
- status = clWaitForEvents(1, &event);
- if (status != CL_SUCCESS) {
- cout << "Wait on event failed, status = " <<
- status << endl;
-
- return NANOTIME_ERR;
- }
-
- event = NULL;
-
-#define TIMING
-#ifdef TIMING
- clFinish( queue);
- time = getCurrentTime();
-
- int iter = 20;
- for ( int i = 1; i <= iter; i++)
- {
-#endif
- status = (cl_int)clMath::clblas::hpr(params_.order, params_.uplo, params_.N, CREAL(alpha_), mobjX_, params_.offBX, params_.incx,
- mobjAP_, params_.offa, 1, &queue, 0, NULL, &event);
-
- if (status != CL_SUCCESS) {
- cerr << "The CLBLAS HPR function failed, status = " <<
- status << endl;
-
- return NANOTIME_ERR;
- }
-
-#ifdef TIMING
- } // iter loop
- clFinish( queue);
- time = getCurrentTime() - time;
- time /= iter;
-#else
- status = flushAll(1, &queue);
- if (status != CL_SUCCESS) {
- cerr << "clFlush() failed, status = " << status << endl;
- return NANOTIME_ERR;
- }
-
- time = getCurrentTime();
- status = waitForSuccessfulFinish(1, &queue, &event);
- if (status == CL_SUCCESS) {
- time = getCurrentTime() - time;
- }
- else {
- cerr << "Waiting for completion of commands to the queue failed, "
- "status = " << status << endl;
- time = NANOTIME_ERR;
- }
-#endif
-
- return time;
-}
-
-} // namespace clMath
-
-TEST_P(HPR, chpr)
-{
- TestParams params;
-
- getParams(&params);
- HprPerformanceTest<FloatComplex>::runInstance(FN_CHPR, &params);
-}
-
-TEST_P(HPR, zhpr)
-{
- TestParams params;
-
- getParams(&params);
- HprPerformanceTest<DoubleComplex>::runInstance(FN_ZHPR, &params);
-}