 
   if( myrow < m && mycol < m )
   {
-    register double temp = 0.0;
-
-/* insert correct index code here */
+    register floatType_t temp = 0.0;
+
     for( int k = 0; k < m; k++ ) 
-    {
       temp += a[INDX( FIXME, FIXME, m )] * b[INDX( FIXME, FIXME, m )];
-    } /* end for */
-
-/* insert index code to write the output to the C matrix */
 
     c[INDX( FIXME, FIXME, m )] = temp;
   } /* end if */


  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");

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

typedef float floatType_t;

/* macro for index calculations */

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