     float log2 = log2_series(data_size);
     double end_time = getclock();
     
-    printf("CPU result: %f\n", log2);
-    printf(" log(2)=%f\n", log(2.0));
-    printf(" time=%fs\n", end_time - start_time);
+//    printf("CPU result: %f\n", log2);
+//    printf(" log(2)=%f\n", log(2.0));
+//    printf(" time=%fs\n", end_time - start_time);
     
     // Parameter definition
     int threads_per_block = 4 * 32;

    float time_per_iter = total_time / (float)data_size;
    float bandwidth = sizeof(float) / time_per_iter; // B/s
    
    printf(" Total time 1: %g s,\n Per iteration: %g ns\n Throughput: %g GB/s\n",
    	total_time,
    	time_per_iter * 1.e9,
    	bandwidth / 1.e9);

    CUDA_SAFE_CALL(cudaEventElapsedTime(&elapsedTime, start_2, stop_2));	// In ms

    total_time = elapsedTime / 1000.;	// s
    time_per_iter = total_time / (float)data_size;
    bandwidth = sizeof(float) / time_per_iter; // B/s
    
    printf(" Total time 2: %g s,\n Per iteration: %g ns\n Throughput: %g GB/s\n",
    	total_time,
    	time_per_iter * 1.e9,
    	bandwidth / 1.e9);
  

    CUDA_SAFE_CALL(cudaEventDestroy(start_2));


    // TODO
    
//    printf("GPU results:\n");
//    printf(" Sum1: %f\n", sum);
//    printf(" Sum2: %f\n", sum_2);
    
    float elapsedTime;
    CUDA_SAFE_CALL(cudaEventElapsedTime(&elapsedTime, start, stop));	// In ms
