Function spr2CorrectnessTest

src/tests/correctness/corr-spr2.cpp:63–191

template <typename T>
void
spr2CorrectnessTest(TestParams *params)
{
    cl_int err;
    T *blasAP, *clblasAP, *X, *Y;
    cl_mem bufAP, bufX, bufY;
    clMath::BlasBase *base;
    cl_event *events;
	bool useAlpha;
	T alpha;

    base = clMath::BlasBase::getInstance();

    if ((typeid(T) == typeid(cl_double)) &&
        !base->isDevSupportDoublePrecision()) {

        std::cerr << ">> WARNING: The target device doesn't support native "
                     "double precision floating point arithmetic" <<
                     std::endl << ">> Test skipped" << std::endl;
        SUCCEED();
        return;
    }

	printf("number of command queues : %d\n\n", params->numCommandQueues);

    events = new cl_event[params->numCommandQueues];
    memset(events, 0, params->numCommandQueues * sizeof(cl_event));

    size_t lengthAP = (params->N *( params->N + 1 ))/2 ;
    size_t lengthX = (1 + ((params->N - 1) * abs(params->incx)));
	size_t lengthY = (1 + ((params->N - 1) * abs(params->incy)));

    blasAP 		= new T[lengthAP + params->offa ];
    clblasAP 	= new T[lengthAP + params->offa ];
    X		 	= new T[lengthX + params->offBX ];
	Y			= new T[lengthY + params->offCY ];

    srand(params->seed);

	if((blasAP == NULL) || (clblasAP == NULL) || (X == NULL) || (Y == NULL))
    {
        deleteBuffers<T>(blasAP, clblasAP, X, Y);
        ::std::cerr << "Cannot allocate memory on host side\n" << "!!!!!!!!!!!!Test skipped.!!!!!!!!!!!!" << ::std::endl;
        delete[] events;
        SUCCEED();
        return;
    }

	alpha =  convertMultiplier<T>(params->alpha);
	useAlpha = true;

    randomSyr2Matrices<T>(params->order, params->uplo, params->N, useAlpha, &alpha, (blasAP + params->offa), params->lda,
							(X + params->offBX), params->incx, (Y + params->offCY), params->incy);

	// Copy blasAP to clblasAP
    memcpy(clblasAP, blasAP, (lengthAP + params->offa)* sizeof(*blasAP));

	// Allocate buffers