summaryrefslogtreecommitdiff
path: root/external/clBLAS/src/tests/performance/perf-trmv.cpp
diff options
context:
space:
mode:
Diffstat (limited to 'external/clBLAS/src/tests/performance/perf-trmv.cpp')
-rw-r--r--external/clBLAS/src/tests/performance/perf-trmv.cpp384
1 files changed, 0 insertions, 384 deletions
diff --git a/external/clBLAS/src/tests/performance/perf-trmv.cpp b/external/clBLAS/src/tests/performance/perf-trmv.cpp
deleted file mode 100644
index 7bb94fe9..00000000
--- a/external/clBLAS/src/tests/performance/perf-trmv.cpp
+++ /dev/null
@@ -1,384 +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.
- * ************************************************************************/
-
-
-/*
- * Trmv performance test cases
- */
-
-#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 <trmv.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 TrmvPerformanceTest : public PerformanceTest
-{
-public:
- virtual ~TrmvPerformanceTest();
-
- virtual int prepare(void);
- virtual nano_time_t etalonPerfSingle(void);
- virtual nano_time_t clblasPerfSingle(void);
-
- static void runInstance(BlasFunction fn, TestParams *params)
- {
- TrmvPerformanceTest<ElemType> perfCase(fn, params);
- int ret = 0;
- int opFactor;
- BlasBase *base;
-
- base = clMath::BlasBase::getInstance();
-
- /* *************Important*********************
- if (fn == FN_STRMV || fn == FN_DTRMV) {
- opFactor = 2;
- }
- else {
- opFactor = 8;
- } this is only for blas-3 routines- operations factor
- FOR BLAS-2(bandwidth intensive) ROUTINES MAKE opFactor AS 1 and pass the appropriate size
- that is read and written in the constructor below */
- opFactor = 1; //FIX-ME
-
- if ((fn == FN_DTRMV || fn == FN_ZTRMV) &&
- !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;
- return;
- }
- 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:
- TrmvPerformanceTest(BlasFunction fn, TestParams *params);
-
- bool areResourcesSufficient(TestParams *params);
-
- TestParams params_;
- ElemType *A_;
- ElemType *X_;
- ElemType *backX_;
- cl_mem mobjA_;
- cl_mem mobjX_;
- cl_mem scratchBuff;
- ::clMath::BlasBase *base_;
-};
-
-template <typename ElemType>
-TrmvPerformanceTest<ElemType>::TrmvPerformanceTest(
- BlasFunction fn,
- TestParams *params) : PerformanceTest( fn, (problem_size_t)( ( params->N * (params->N+1) * sizeof(ElemType) ) ) ), //**************Gbps formula here***********
- params_(*params), mobjA_(NULL), mobjX_(NULL)
-{
-
- A_ = new ElemType[params_.N * params_.lda + params_.offa];
- X_ = new ElemType[ 1 + (params_.N-1) * abs(params_.incx) + params_.offBX];
- backX_ = new ElemType[ 1 + (params_.N-1) * abs(params_.incx) + params_.offBX];
-
- base_ = ::clMath::BlasBase::getInstance();
- mobjA_ = NULL;
- mobjX_ = NULL;
- scratchBuff = NULL;
-}
-
-template <typename ElemType>
-TrmvPerformanceTest<ElemType>::~TrmvPerformanceTest()
-{
- if(A_ != NULL)
- {
- delete[] A_;
- }
- if(X_ != NULL)
- {
- delete[] X_;
- }
- if(backX_ != NULL)
- {
- delete[] backX_;
- }
-
- if ( mobjA_ != NULL )
- clReleaseMemObject(mobjA_);
- if ( mobjX_ != NULL )
- clReleaseMemObject(mobjX_);
- if ( scratchBuff != NULL )
- clReleaseMemObject(scratchBuff);
-}
-
-/*
- * Check if available OpenCL resources are sufficient to
- * run the test case
- */
-template <typename ElemType> bool
-TrmvPerformanceTest<ElemType>::areResourcesSufficient(TestParams *params)
-{
- clMath::BlasBase *base;
- size_t gmemSize, allocSize;
- size_t n = params->N;
-
- if((A_ == NULL) || (X_ == NULL) || (backX_ == NULL))
- {
- return 0;
- }
-
- base = clMath::BlasBase::getInstance();
- gmemSize = (size_t)base->availGlobalMemSize( 0 );
- allocSize = (size_t)base->maxMemAllocSize();
-
- bool suff = ( sizeof(ElemType)*n*params->lda < allocSize ) && ((1 + (n-1)*abs(params->incx))*sizeof(ElemType) < allocSize); //for individual allocations
- suff = suff && ((( n*params->lda + (1 + (n-1)*abs(params->incx))*2)*sizeof(ElemType)) < gmemSize) ; //for total global allocations
-
- return suff ;
-}
-
-template <typename ElemType> int
-TrmvPerformanceTest<ElemType>::prepare(void)
-{
- size_t lenX, N;
- N = params_.N;
- lenX = 1 + (N-1) * abs(params_.incx);
-
-
- int creationFlags = 0;
- creationFlags = creationFlags | RANDOM_INIT;
-
- // Default is Column-Major
- creationFlags = ( (this-> params_.order) == clblasRowMajor)? (creationFlags | ROW_MAJOR_ORDER) : (creationFlags);
- creationFlags = ( (this-> params_.uplo) == clblasLower)? (creationFlags | LOWER_HALF_ONLY) : (creationFlags | UPPER_HALF_ONLY);
- BlasRoutineID BlasFn = CLBLAS_TRMV;
-
- // Matrix A
- populate( (A_ + params_.offa), N, N, params_.lda, BlasFn, creationFlags);
- populate( X_ , lenX + params_.offBX, 1, lenX + params_.offBX, BlasFn);
- memcpy(backX_, X_, ((1 + (params_.N-1) * abs(params_.incx))+ params_.offBX )* sizeof(ElemType));
-
-
- mobjA_ = base_->createEnqueueBuffer(A_, (params_.N * params_.lda + params_.offa)* sizeof(*A_), 0, CL_MEM_READ_ONLY);
- mobjX_ = base_->createEnqueueBuffer(X_, (lenX + params_.offBX )* sizeof(*X_), 0, CL_MEM_WRITE_ONLY);
- scratchBuff = base_->createEnqueueBuffer(NULL , lenX * sizeof(*X_), 0, CL_MEM_READ_ONLY);
-
- return ( (mobjA_ != NULL) && (mobjX_ != NULL) && (scratchBuff != NULL) ) ? 0 : -1;
-}
-
-template <typename ElemType> nano_time_t
-TrmvPerformanceTest<ElemType>::etalonPerfSingle(void)
-{
- nano_time_t time = 0;
- clblasOrder order;
- clblasUplo fUplo;
- clblasTranspose fTrans;
- 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;
- fUplo = params_.uplo;
- fTrans = params_.transA;
- lda = params_.lda;
-
-#ifdef PERF_TEST_WITH_ACML
-
- if (order != clblasColumnMajor)
- {
- order = clblasColumnMajor;
- fUplo = (params_.uplo == clblasUpper)? clblasLower : clblasUpper;
- fTrans = (params_.transA == clblasNoTrans)? clblasTrans : clblasNoTrans;
-
- if( params_.transA == clblasConjTrans )
- doConjugate( (A_+params_.offa), params_.N, params_.N, lda );
- }
-
-
- time = getCurrentTime();
- clMath::blas::trmv(order, fUplo,fTrans, params_.diag,
- params_.N, A_, params_.offa, lda, X_, params_.offBX, params_.incx);
- time = getCurrentTime() - time;
-
-#endif // PERF_TEST_WITH_ACML
-
- return time;
-}
-
-
-template <typename ElemType> nano_time_t
-TrmvPerformanceTest<ElemType>::clblasPerfSingle(void)
-{
- nano_time_t time;
- cl_event event;
- cl_int status;
- cl_command_queue queue = base_->commandQueues()[0];
- size_t lenX = 1 + (params_.N-1) * abs(params_.incx);
-
- status = clEnqueueWriteBuffer(queue, mobjX_, CL_TRUE, 0,
- (lenX + params_.offBX )* sizeof(ElemType), backX_, 0, NULL, &event);
- if (status != CL_SUCCESS) {
- cerr << "Vector X 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;
-
- DataType type;
- type = ( typeid(ElemType) == typeid(float))? TYPE_FLOAT:( typeid(ElemType) == typeid(double))? TYPE_DOUBLE:
- ( typeid(ElemType) == typeid(FloatComplex))? TYPE_COMPLEX_FLOAT: TYPE_COMPLEX_DOUBLE;
- time = getCurrentTime();
-#define TIMING
-#ifdef TIMING
- clFinish( queue);
-
- int iter = 20;
- for ( int i = 1; i <= iter; i++)
- {
-#endif
- status = (cl_int)clMath::clblas::trmv(type, params_.order, params_.uplo,
- params_.transA, params_.diag, params_.N, mobjA_, params_.offa, params_.lda,
- mobjX_, params_.offBX, params_.incx, scratchBuff,
- 1, &queue, 0, NULL, &event);
-
-
- if (status != CL_SUCCESS) {
- cerr << "The CLBLAS TRMV 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;
- }
-
- //printf("Time elapsed : %lu\n", time);
-#endif
-
- return time;
-}
-
-} // namespace clMath
-
-// strmv performance test
-TEST_P(TRMV, strmv)
-{
- TestParams params;
-
- getParams(&params);
- TrmvPerformanceTest<float>::runInstance(FN_STRMV, &params);
-}
-
-// dtrmv performance test case
-TEST_P(TRMV, dtrmv)
-{
- TestParams params;
-
- getParams(&params);
- TrmvPerformanceTest<double>::runInstance(FN_DTRMV, &params);
-}
-// ctrmv performance test case
-TEST_P(TRMV, ctrmv)
-{
- TestParams params;
-
- getParams(&params);
- TrmvPerformanceTest<FloatComplex>::runInstance(FN_CTRMV, &params);
-}
-// ztrmv performance test case
-TEST_P(TRMV, ztrmv)
-{
- TestParams params;
-
- getParams(&params);
- TrmvPerformanceTest<DoubleComplex>::runInstance(FN_ZTRMV, &params);
-}
-
generated by cgit v1.2.3 (git 2.39.1) at 2024-08-02 23:31:46 +0000