     cublasHandle_t handle;
     checkCUBLAS( cublasCreate( &handle ) );
 
-    double alpha = 1.0;
-    double beta  = 0.0;
+    floatType_t alpha = 1.0;
+    floatType_t beta  = 0.0;
 
 	/* start timers */
 

#include "cublas_v2.h"
#include "../debug.h"

/* macro for index calculations */

#define INDX( row, col, ld ) ( ( (col) * (ld) ) + (row) )

    checkCUDA( cudaEventDestroy( start ) );
    checkCUDA( cudaEventDestroy( stop ) );

/* check CUBLAS versus GPU NAIVE numerical results */

    double temp = 0.0;

    for( int i = 0; i < size * size; i++ )
    {
       temp = max( temp, abs( (double)h_c[i] - (double)h_c1[i] )/
                      abs((double)h_c[i]) );
    } /* end for */
    printf("Maximum error is %e percent \n",temp*100.0);
    if( temp > 0.001 ) printf("FAIL\n");
    else printf("PASS\n");

/* cleanup */


    checkCUDA( cudaFree( d_a ) );
    checkCUDA( cudaFree( d_b ) );
