 /* read block of B into shared memory */
     bs[ FIXME ][ FIXME ] = b[ boff ];
 		
+    __syncthreads();
 
 /* increment A and B offsets  for next round of data reads */
     boff += BLOCK_K;


/* call CUBLAS dgemm */

  checkCUBLAS( 
  cublasDgemm( handle, CUBLAS_OP_N, CUBLAS_OP_N,
               size, size, size,
               &alpha, 
               d_a, size,
               d_b, size,
               &beta,
               d_c, size )
             );

/* stop timers */


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

typedef float floatType_t;

/* macro for index calculations */

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