 
 /* determine my threads's row and col indices in the global C matrix */
 
-	const int myrow = blockDim.x * blockIdx.x + threadIdx.x;
-	const int mycol = blockDim.y * blockIdx.y + threadIdx.y;
+  const int myrow = blockDim.x * blockIdx.x + threadIdx.x;
+  const int mycol = blockDim.y * blockIdx.y + threadIdx.y;
 
 /* if my row and col are in the C matrix, then calculate that value of C */
 
-	if( myrow < m && mycol < m )
-	{
-		register double temp = 0.0;
-		for( int k = 0; k < m; k++ ) temp += a[INDX( myrow, k, m )] * b[INDX( k, mycol, m )];
-		c[INDX( myrow, mycol, m )] = temp;
-	} /* end if */
+  if( myrow < m && mycol < m )
+  {
+    register double temp = 0.0;
+
+    for( int k = 0; k < m; k++ ) 
+      temp += a[INDX( myrow, k, m )] * b[INDX( k, mycol, m )];
+
+    c[INDX( myrow, mycol, m )] = temp;
+  } /* end if */
 
 	return;
 } /* end GPU_naive */

/* matrix size and thread dimensions */

#define SIZE 1024
#define THREAD_X 16
#define THREAD_Y 16

/* naive GPU kernel where each element of C is computed by a single thread */


#define CUDA_CHECK() 
#endif

/* macro for index calculations */

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