 
 /* shared memory arrays for A and B */
 
-  __shared__ double as[ THREADS_PER_BLOCK_X ][ BLOCK_K + 1 ];
-  __shared__ double bs[ BLOCK_K ][ THREADS_PER_BLOCK_Y + 1 ];
+  __shared__ floatType_t as[ THREADS_PER_BLOCK_X ][ BLOCK_K + 1 ];
+  __shared__ floatType_t bs[ BLOCK_K ][ THREADS_PER_BLOCK_Y + 1 ];
 	
 /* space for C to be held in registers */
 
-  double c_tmp = 0.0 ;
+  floatType_t c_tmp = 0.0 ;
 
 /* calculate my initial offset into A and B */
 

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

/* macro for index calculations */

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


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