         gold<<<n_blocks, block_size>>>(thrust::raw_pointer_cast(s.data()), thrust::raw_pointer_cast(g.data()));
         correct = thrust::equal(r.begin(), r.end(), g.begin());
     }
+    
     if (correct)
         std::cout << "Results passed";
     else

    int block_size = 256;
    int n = n_blocks * block_size;
    thrust::device_vector<T> s(n);
    thrust::device_vector<T> r(n);
    int iterations = 10;
    cudaEvent_t start, stop;

    int block_size = 256;
    int n = n_blocks * block_size;
    thrust::device_vector<T> s(n);
    fill_test(s);
    thrust::device_vector<T> r(n);
    int iterations = 10;
    cudaEvent_t start, stop;
