 
   if( myrow < m && mycol < m )
   {
-    register double temp = 0.0;
+    register floatType_t temp = 0.0;
 
     for( int k = 0; k < m; k++ ) 
       temp += a[INDX( myrow, k, m )] * b[INDX( k, mycol, m )];


  for( int i = 0; i < size * size; i++ )
  {
    temp += ( h_c[i] - h_c1[i] ) * ( h_c[i] - h_c1[i] );
  } /* end for */

  printf("error is %f\n",temp);
  if( temp > 10 ) printf("FAIL\n");
  else printf("PASS\n");

/* cleanup */

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

typedef float floatType_t;

/* macro for index calculations */

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