 #include "custom_cutil_math.h"
 #include "SPHCudaUtils.h"
-#include "omp.h"
+#include <sys/time.h>
 #include <thrust/sort.h>
 #include <thrust/scan.h>
 #include <thrust/reduce.h>

		//************
		myGpuTimer.Stop();
		real_ time2 = (real_)myGpuTimer.Elapsed();
		real_ ompTimerEnd = (real_)omp_get_wtime();


		if (tStep % 2 == 0) {
			printf("step: %d, step Time (CUDA): %f, step Time (OMP): %f\n ", tStep, time2, 1000 * (ompTimerEnd - ompTimerStart));
			//printf("a \n");
		}
		fflush(stdout);


		GpuTimer myGpuTimer;
		myGpuTimer.Start();

		struct timeval cpuT_start, cpuT_end;
		struct timezone cpuT_timezone;
		gettimeofday(&cpuT_start, &cpuT_timezone);

//		if (tStep < 1000) delT = 0.25 * delTOrig; else delT = delTOrig;

