Combine reduction kernels

Author: sami-hatna66

external/Stats.cpp

@@ -44,27 +44,16 @@ void updateStats(sycl::queue myQueue, sycl::buffer<Actor> actorBuf,
 
             auto forceSumReduction =
                 sycl::reduction(forceSumBuf, cgh, sycl::plus<float>());
-
-            cgh.parallel_for(sycl::range<1>{actorAcc.size()}, forceSumReduction,
-                             [=](sycl::id<1> index, auto &sum) {
-                                 if (!actorAcc[index].getAtDestination()) {
-                                     sum += actorAcc[index].getForce();
-                                 }
-                             });
-        });
-        myQueue.throw_asynchronous();
-
-        myQueue.submit([&](sycl::handler &cgh) {
-            auto actorAcc = actorBuf.get_access<sycl::access::mode::read>(cgh);
-
+            
             auto activeActorsReduction =
                 sycl::reduction(activeActorsBuf, cgh, sycl::plus<int>());
 
-            cgh.parallel_for(sycl::range<1>{actorAcc.size()},
-                             activeActorsReduction,
-                             [=](sycl::id<1> index, auto &sum) {
-                                 if (!actorAcc[index].getAtDestination()) {
-                                     sum += 1;
+            int globalSize = (std::floor(actorAcc.size()/16) + 1) * 16;
+            cgh.parallel_for(sycl::nd_range<1>{globalSize, 16}, forceSumReduction, activeActorsReduction,
+                             [=](sycl::nd_item<1> index, auto &forceSum, auto& actorSum) {
+                                 if (!actorAcc[index.get_global_id()].getAtDestination()) {
+                                     forceSum += actorAcc[index.get_global_id()].getForce();
+                                     actorSum += 1;
                                  }
                              });
         });