#include #include #include #include #include #include #include #include // this example illustrates how to make repeated access to a range of values // examples: // repeated_range([0, 1, 2, 3], 1) -> [0, 1, 2, 3] // repeated_range([0, 1, 2, 3], 2) -> [0, 0, 1, 1, 2, 2, 3, 3] // repeated_range([0, 1, 2, 3], 3) -> [0, 0, 0, 1, 1, 1, 2, 2, 2, 3, 3, 3] // ... template class repeated_range { public: typedef typename thrust::iterator_difference::type difference_type; struct repeat_functor : public thrust::unary_function { difference_type repeats; repeat_functor(difference_type repeats) : repeats(repeats) {} __host__ __device__ difference_type operator()(const difference_type& i) const { return i / repeats; } }; typedef typename thrust::counting_iterator CountingIterator; typedef typename thrust::transform_iterator TransformIterator; typedef typename thrust::permutation_iterator PermutationIterator; // type of the repeated_range iterator typedef PermutationIterator iterator; // construct repeated_range for the range [first,last) repeated_range(Iterator first, Iterator last, difference_type repeats) : first(first), last(last), repeats(repeats) {} iterator begin(void) const { return PermutationIterator(first, TransformIterator(CountingIterator(0), repeat_functor(repeats))); } iterator end(void) const { return begin() + repeats * (last - first); } protected: Iterator first; Iterator last; difference_type repeats; }; int main(void) { thrust::device_vector data(4); data[0] = 10; data[1] = 20; data[2] = 30; data[3] = 40; // print the initial data std::cout << "range "; thrust::copy(data.begin(), data.end(), std::ostream_iterator(std::cout, " ")); std::cout << std::endl; typedef thrust::device_vector::iterator Iterator; // create repeated_range with elements repeated twice repeated_range twice(data.begin(), data.end(), 2); std::cout << "repeated x2: "; thrust::copy(twice.begin(), twice.end(), std::ostream_iterator(std::cout, " ")); std::cout << std::endl; // create repeated_range with elements repeated x3 repeated_range thrice(data.begin(), data.end(), 3); std::cout << "repeated x3: "; thrust::copy(thrice.begin(), thrice.end(), std::ostream_iterator(std::cout, " ")); std::cout << std::endl; return 0; }