Optimize gather_tracks

Author: beomki-yeo

device/common/include/traccc/ambiguity_resolution/device/gather_tracks.hpp

@@ -41,6 +41,11 @@ struct gather_tracks_payload {
      */
     vecmem::data::vector_view<unsigned int> sorted_ids_view;
 
+    /**
+     * @brief View object to the updated track
+     */
+    vecmem::data::vector_view<unsigned int> updated_tracks_view;
+
     /**
      * @brief View object to the whether track id is updated
      */

device/cuda/src/ambiguity_resolution/greedy_ambiguity_resolution_algorithm.cu

@@ -412,9 +412,6 @@ greedy_ambiguity_resolution_algorithm::operator()(
     unsigned int nBlocks_warp =
         (n_accepted + nThreads_warp - 1) / nThreads_warp;
 
-    unsigned int nThreads_full = 1024;
-    unsigned int nBlocks_full = (n_tracks + 1023) / 1024;
-
     // Compute the threadblock dimension for scanning kernels
     auto compute_scan_config = [&](unsigned int n_accepted) {
         unsigned int nThreads_scan = m_warp_size * 4;
@@ -580,13 +577,14 @@ greedy_ambiguity_resolution_algorithm::operator()(
             .temp_sorted_ids_view = temp_sorted_ids_buffer,
         });
 
-        kernels::gather_tracks<<<nBlocks_full, nThreads_full, 0, stream>>>(
+        kernels::gather_tracks<<<nBlocks_adaptive, nThreads_adaptive, 0, stream>>>(
             device::gather_tracks_payload{
                 .terminate = terminate_device.get(),
                 .n_accepted = n_accepted_device.get(),
                 .n_updated_tracks = n_updated_tracks_device.get(),
                 .temp_sorted_ids_view = temp_sorted_ids_buffer,
                 .sorted_ids_view = sorted_ids_buffer,
+                .updated_tracks_view = updated_tracks_buffer,
                 .is_updated_view = is_updated_buffer});
 
         cudaStreamEndCapture(stream, &graph);

device/cuda/src/ambiguity_resolution/kernels/gather_tracks.cu

@@ -23,14 +23,16 @@ __global__ void gather_tracks(device::gather_tracks_payload payload) {
     vecmem::device_vector<const unsigned int> temp_sorted_ids(
         payload.temp_sorted_ids_view);
     vecmem::device_vector<unsigned int> sorted_ids(payload.sorted_ids_view);
+    vecmem::device_vector<unsigned int> updated_tracks(
+        payload.updated_tracks_view);
     vecmem::device_vector<int> is_updated(payload.is_updated_view);
 
     auto globalIndex = threadIdx.x + blockIdx.x * blockDim.x;
     const unsigned int n_accepted = *(payload.n_accepted);
 
     // Reset is_updated vector
-    if (globalIndex < is_updated.size()) {
-        is_updated[globalIndex] = 0;
+    if (globalIndex < *(payload.n_updated_tracks)) {
+        is_updated[updated_tracks[globalIndex]] = 0;
     }
 
     if (globalIndex >= n_accepted) {