     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;




template<int i>
void run_benchmark_contiguous_load_store(const std::string name, void (*test)(array<int, i>*, array<int, i>*),
                                         void (*gold)(array<int, i>*, array<int, i>*)) {

        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
