 #include "cublas_v2.h"
 #include "../debug.h"
 
+typedef float floatType_t;
+
 /* macro for index calculations */
 
 #define INDX( row, col, ld ) ( ( (col) * (ld) ) + (row) )

  cublasHandle_t handle;
  checkCUBLAS( cublasCreate( &handle ) );

  double alpha = 1.0;
  double beta  = 0.0;

/* start timers */



/* call CUBLAS dgemm */

  if( sizeof( floatType_t ) == 4 ) 
  {
  checkCUBLAS( 
  cublasSgemm( handle, CUBLAS_OP_N, CUBLAS_OP_N,
               size, size, size,
               (float *)&alpha, 
               (float *)d_a, size,
               (float *)d_b, size,
               (float *)&beta,
               (float *)d_c, size )
              );
  } /* end if */
  else
  {
  checkCUBLAS( 
  cublasDgemm( handle, CUBLAS_OP_N, CUBLAS_OP_N,
               size, size, size,
               (double *)&alpha, 
               (double *)d_a, size,
               (double *)d_b, size,
               (double *)&beta,
               (double *)d_c, size )
              );
  } /* end else */

/* stop timers */

