Function randomSpmvMatrices

src/tests/include/spmv.h:128–205

template <typename T>
static void
randomSpmvMatrices(
    clblasOrder order,
    clblasUplo uplo,
    size_t N,
    bool useAlpha,
    T *alpha,
    T *A,
    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(*alpha) == 0.0) {
            *alpha = 1.0;
        }
    }

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

    bound = UPPER_BOUND<T>();

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

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

	fAlpha = module(*alpha);
	fBeta  = module(*beta);

    bound = bound / (fAlpha * N);

    bound = sqrt( ((((((fBeta * fBeta)) / fAlpha) / (4.0)) / fAlpha) / (N * N)) + bound) - ((fBeta) / ((2.0) * (fAlpha) * N));