   cudaFree(g_d);
   //free up all unused user allocated memory on Cuda device
 
-  printf("\n");  
+  //printf("\n");  
 
   /*The next 2 for loops print out the values of both
   answer matrices. This can be used to ensure that both

  //smultiply<<<1,1>>>(g_a, g_b, g_c);
  //cudaEventRecord(stop, 0);  

  cudaEventSynchronize(stop);
  cudaEventElapsedTime(&time, start, stop);
  //get run time

//  cudaMemcpy(c, g_c, g_size, cudaMemcpyDeviceToHost);

      }
    }
  }
}

//parallel matrix multiplication kernel
__global__ void pmultiply(double* g_a, double* g_b, double* g_d , int dim)
