fix

alpaka-group · Nov 29, 2021 · 2bcd53d · 2bcd53d
1 parent 8942dfa
commit 2bcd53d
Showing 1 changed file with 27 additions and 15 deletions.
diff --git a/examples/thrust/thrust.cu b/examples/thrust/thrust.cu
@@ -321,6 +321,29 @@ auto thrustDeviceAlloc = [](auto alignment, std::size_t size)
     return p;
 };
 
+template <typename View>
+struct ViewIteratorAt
+{
+    View view;
+
+    LLAMA_FN_HOST_ACC_INLINE auto operator()(std::size_t i)
+    {
+        return *(view.begin() + i);
+    }
+};
+
+template<typename View>
+auto viewIteratorAt(View& view, std::size_t index)
+{
+    ViewIteratorAt<View> t{view};
+    using ViewTransformIterator = thrust::transform_iterator<
+        decltype(t),
+        thrust::counting_iterator<std::size_t>,
+        typename View::iterator::reference,
+        typename View::iterator::value_type>;
+    return ViewTransformIterator{thrust::counting_iterator<std::size_t>{index}, t};
+}
+
 template<int Mapping>
 void run(std::ostream& plotFile)
 {
@@ -374,15 +397,8 @@ void run(std::ostream& plotFile)
     std::cout << mappingName << '\n';
 
     auto view = llama::allocView(mapping, thrustDeviceAlloc);
-
-    auto makeViewIteratorFromIndexCreator = [](decltype(view) view)
-    { return [view] __host__ __device__(std::size_t i) mutable { return *(view.begin() + i); }; };
-    auto b = thrust::make_transform_iterator(
-        thrust::counting_iterator<std::size_t>{0},
-        makeViewIteratorFromIndexCreator(view));
-    auto e = thrust::make_transform_iterator(
-        thrust::counting_iterator<std::size_t>{N},
-        makeViewIteratorFromIndexCreator(view));
+    auto b = viewIteratorAt(view, 0);
+    auto e = viewIteratorAt(view, N);
     // auto b = view.begin();
     // auto e = view.end();
 
@@ -541,9 +557,7 @@ void run(std::ostream& plotFile)
 
         {
             auto dstView = llama::allocView(mapping, thrustDeviceAlloc);
-            auto db = thrust::make_transform_iterator(
-                thrust::counting_iterator<std::size_t>{0},
-                makeViewIteratorFromIndexCreator(dstView));
+            auto db = viewIteratorAt(dstView, 0);
             Stopwatch stopwatch;
             if constexpr(usePSTL)
                 std::copy(exec, b, e, db);
@@ -559,9 +573,7 @@ void run(std::ostream& plotFile)
 
         {
             auto dstView = llama::allocView(mapping, thrustDeviceAlloc);
-            auto db = thrust::make_transform_iterator(
-                thrust::counting_iterator<std::size_t>{0},
-                makeViewIteratorFromIndexCreator(dstView));
+            auto db = viewIteratorAt(dstView, 0);
             Stopwatch stopwatch;
             if constexpr(usePSTL)
                 std::copy_if(exec, b, e, db, Predicate{});