 	int *V, *out;				//host copies
 	int *d_V, *d_out;			//device copies
 	int *D;						
-
-//int *h_D;
-
 	FILE *fp;
 	if(argc != 2)
 	{

		float time;
		// get event elapsed time
		cudaEventElapsedTime(&time, start, stop);

		// copy result back to host
		cudaMemcpy(out, d_out, m*k*sizeof(int), cudaMemcpyDeviceToHost);

//		// record event and synchronize
//		cudaEventRecord(stop);
//		cudaEventSynchronize(stop);
//		float time;
//		// get event elapsed time
//		cudaEventElapsedTime(&time, start, stop);

		showResult(m, k, out);
		if(m == 1024) {

#define TILE_WIDTH 32
#define TILE_DEPTH 128
#define MAX_BLOCK_SIZE 256
#define MAX_PTRNUM_IN_SMEM 1024 

void showResult(int m, int k, int *out);

// compute the square of distance of the ith point and jth point
__global__ void computeDist(int m, int n, int *V, int *D)
