   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 */
 


  fprintf(stdout, "Matrix size is %d\n",size);

  double *h_a, *h_b, *h_c, *h_c1;
  double *d_a, *d_b, *d_c;
 
  size_t numbytes = (size_t ) size * (size_t ) size * sizeof( double );

  h_a = (double *) malloc( numbytes );
  if( h_a == NULL )
  {
    fprintf(stderr,"Error in host malloc\n");
    return 911;
  }

  h_b = (double *) malloc( numbytes );
  if( h_b == NULL )
  {
    fprintf(stderr,"Error in host malloc\n");
    return 911;
  }

  h_c = (double *) malloc( numbytes );
  if( h_c == NULL )
  {
    fprintf(stderr,"Error in host malloc\n");
    return 911;
  }

  h_c1 = (double *) malloc( numbytes );
  if( h_c1 == NULL )
  {
    fprintf(stderr,"Error in host malloc\n");


/* 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 */

