 
   hysteresis(d_image,size,t1,t2);
 
-//  gpuOutput(d_image,image,size.z);
   cudaMemcpy(image,d_image,size.z*sizeof(float),cudaMemcpyDeviceToHost);
 
   /// free memory used for image magnitude

  cudaMalloc((void**) &d_magnitude, size.z*sizeof(float));
  CUT_CHECK_ERROR("Memory temporary image creation failed");

  cudaSobel(d_magnitude,d_direction,d_image,size);

  gradientMaximumDetector(d_image,d_magnitude,d_direction,size,maxKernelWidth);


  size.z = width*height;
  int imageSize = size.z*sizeof(float);

  /// Warmup
  unsigned int WarmupTimer = 0;
  cutCreateTimer( &WarmupTimer );
  cutStartTimer( WarmupTimer );
