Simple sequential thrust program

gistfile1.cpp

//Define the system to be OpenMP to remove all CUDA references
#define THRUST_DEVICE_SYSTEM THRUST_DEVICE_SYSTEM_OMP
#include <thrust/host_vector.h>
#include <thrust/fill.h>
#include <thrust/reduce.h>
#include <iostream>

template<typename T>
struct my_reduction {
    typedef T result_type;
    T operator()(const T& a, const T& b) const {
        return a + b + 1;
    }
};

int main() {
    int length = 10;
    thrust::host_vector<int> x(length);
    thrust::fill(x.begin(), x.end(), 1);
    int reduction_value = thrust::reduce(x.begin(), x.end(), 0, my_reduction<int>());
    std::cout << "Thrust result: " << reduction_value << std::endl;
    std::cout << "Result should be: " << 2 * length << std::endl;
}

Thank's I'll see what I can do with this.