 
 		GpuTimer myGpuTimer;
 		myGpuTimer.Start();
-		real_ ompTimerStart = (real_)omp_get_wtime();
-
+
+		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;
 

		//************
		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);

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