Function setDiagonalUnityOrNonUnity

src/tests/include/blas-random.h:803–860

template <typename T>
static void
setDiagonalUnityOrNonUnity(int unity, T* data, size_t rows, size_t cols, size_t lda, int vectorLength, int creationFlags, cl_double bound)
{

	if (creationFlags & PACKED_MATRIX)
	{

		// Rows = Cols for PACKED Matrix
		for(size_t i=0;i< rows;i++)
		{
			if (creationFlags & UPPER_HALF_ONLY)
    			{
				(unity==1)? setElementWithUnity( ((creationFlags & ROW_MAJOR_ORDER))?RMUPacked(i,i):RMLPacked(i,i), vectorLength):
                      			(unity == 0)? setElementWithZero( ((creationFlags & ROW_MAJOR_ORDER))?RMUPacked(i,i):RMLPacked(i,i), vectorLength):
                      				setElementWithRandomData( ((creationFlags & ROW_MAJOR_ORDER))?RMUPacked(i,i):RMLPacked(i,i), vectorLength, bound);
    			}
    			else
    			{
          			(unity==1)? setElementWithUnity( (creationFlags & ROW_MAJOR_ORDER)?RMLPacked(i,i):RMUPacked(i,i), vectorLength):
                      			(unity==0)? setElementWithZero( (creationFlags & ROW_MAJOR_ORDER)?RMLPacked(i,i):RMUPacked(i,i), vectorLength):
                      				setElementWithRandomData( (creationFlags & ROW_MAJOR_ORDER)?RMLPacked(i,i):RMUPacked(i,i) , vectorLength, bound);
    			}
		}
	}
	else
	{
		// Row Major - rows x lda
		// Col major - lda x cols
		size_t firstdimension;
		T *p;

		if (creationFlags & ROW_MAJOR_ORDER)
		{
			firstdimension = rows;
		} else {
			firstdimension = cols;
		}

		for(size_t i=0; i<firstdimension; i++)
		{
			p = (T *)data + (i*lda)*vectorLength;
			p += i*vectorLength;

			if (unity == 0)
			{
				setElementWithZero(p, vectorLength);
			}
			else if (unity == 1)
			{
				setElementWithUnity(p, vectorLength);
			}
			else
			{
				setElementWithRandomData(p, vectorLength, bound);
			}
		}
	}
}