Function randomSyrMatrices

src/tests/include/blas-random.h:492–560

template <typename T>
static void
randomSyrMatrices(
    clblasOrder order,
    clblasUplo uplo,
    size_t N,
    bool useAlpha,
    T *alpha,
    T *A,
    size_t lda,
    T *X,
	int incx
    )
{
    size_t i, j;
	size_t lengthX;
    cl_double bound;

    if (!useAlpha) {
        *alpha = random<T>(100);
        if (module(*alpha) == 0.0) {
            *alpha = ONE<T>();
        }
    }
	#ifdef DEBUG_SYR
	printf("ALPHA in randomSyrMatrices %f\n", *alpha);
	#endif

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

	bound = UPPER_BOUND<T>();

	if(module(*alpha) > (sqrt(bound) / (2.0)))
		*alpha = random<T>((sqrt(bound) / (2.0)));

	#ifdef DEBUG_SYR
	printf("ALPHA in randomSyrMatrices after check %f bound for alpha is %f\n", *alpha, (sqrt(bound) / (2.0)));
	#endif

	bound = bound / module(*alpha);

    bound = sqrt( ((((1.0) / module(*alpha)) / (4.0)) / module(*alpha)) + bound) - ((1.0) / ((2.0) * (*alpha)));

	#ifdef DEBUG_SYR
	printf("BOUND : %f alpha %f \n", bound, *alpha);
	#endif

     if( lda )
    {
    for (i = 0; i < N; i++) {
        for (j = 0; j < N; j++) {
            setElement<T>(order, clblasNoTrans, i, j, A, lda, random<T>(bound));
        }
    }
    } else {
        for (i = 0; i < N; i++) {
            for (j = 0; j < N; j++) {
                setElementPacked<T>(order, clblasNoTrans, uplo, i, j, A, N, random<T>(bound));
            }