   //smultiply<<<1,1>>>(g_a, g_b, g_c);
   //cudaEventRecord(stop, 0);  
 
-  cudaEventSynchronize(stop);
-  cudaEventElapsedTime(&time, start, stop);
+  //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)

  cudaFree(g_d);
  //free up all unused user allocated memory on Cuda device

  //printf("\n");  

  /*The next 2 for loops print out the values of both
  answer matrices. This can be used to ensure that both
