Function randomHemvMatrices

src/tests/include/blas-random.h:642–742

template <typename T>
static void
randomHemvMatrices(
    clblasOrder order,
    clblasUplo uplo,
    size_t N,
    bool useAlpha,
    T *alpha,
    T *A,
    size_t lda,
    T *X,
    int incx,
	bool useBeta,
	T *beta,
    T *Y,
    int incy
    )
{
    size_t i, j;
    size_t lengthX;
    size_t lengthY;
    cl_double bound;
	cl_double fAlpha, fBeta;

    if (!useAlpha) {
        *alpha = random<T>(100);
        if (module(CREAL(*alpha)) == 0.0) {
            CREAL(*alpha) = 1.0;
        }
    }

	if (!useBeta) {
        *beta = random<T>(100);
        if (module(CREAL(*beta)) == 0.0) {
            CREAL(*beta) = 1.0;
        }
    }

    #ifdef DEBUG_HEMV
    printf("ALPHA in randomSyr2Matrices %f.%f\n", CREAL(*alpha), CIMAG(*alpha));
    printf("BETA in randomSyr2Matrices %f.%f\n", CREAL(*beta), CIMAG(*beta));
    #endif

    // bound is calculated by solving the equation (2*alpha*x^2 + x - UPPER_BOUND) < 0

    bound = UPPER_BOUND<T>();

    if((module(CREAL(*alpha)) > bound) || (module(CIMAG(*alpha)) > bound))
        *alpha = random<T>((sqrt(bound) / ((2.0) * N)));
	if (module(CREAL(*alpha)) == 0.0) {
            CREAL(*alpha) = 1.0;
    }

	if((module(CREAL(*beta)) > bound) || (module(CIMAG(*beta)) > bound))
        *beta = random<T>((sqrt(bound)));
	if (module(CREAL(*beta)) == 0.0) {
            CREAL(*beta) = 1.0;
    }