Protect GPU kernels with appropriate nemespaces.

resolve/cuda/CMakeLists.txt

@@ -8,13 +8,12 @@
 
 set(ReSolve_CUDA_SRC
     cudaKernels.cu
-    cudaVectorKernels.cu
+    VectorKernels.cu
     MemoryUtils.cu
 )
 
 set(ReSolve_CUDA_HEADER_INSTALL
     cudaKernels.h
-    cudaVectorKernels.h
     CudaMemory.hpp
     cuda_check_errors.hpp
 )

resolve/cuda/VectorKernels.cu

@@ -0,0 +1,47 @@
+/**
+ * @file cudaVectorKernels.cu
+ * @author Slaven Peles ([email protected])
+ * @brief Contains implementation of CUDA vector kernels and their wrappers.
+ * @date 2023-12-08
+ * 
+ * @note Kernel wrappers implemented here are intended for use in hardware
+ * agnostic code.
+ */
+
+#include <resolve/vector/VectorKernels.hpp>
+
+
+namespace ReSolve { namespace vector {
+
+  namespace kernels {
+
+    /**
+     * @brief CUDA kernel that sets values of an array to a constant.
+     *
+     * @param[in]  n   - length of the array
+     * @param[in]  val - the value the array is set to
+     * @param[out] arr - a pointer to the array
+     * 
+     * @pre  `arr` is allocated to size `n`
+     * @post `arr` elements are set to `val`
+     */
+    __global__ void set_const(index_type n, real_type val, real_type* arr)
+    {
+      index_type i = blockIdx.x * blockDim.x + threadIdx.x;
+      if(i < n)
+      {
+        arr[i] = val;
+      }
+    }
+
+  } // namespace kernels
+
+  void set_array_const(index_type n, real_type val, real_type* arr)
+  {
+    index_type num_blocks;
+    index_type block_size = 512;
+    num_blocks = (n + block_size - 1) / block_size;
+    kernels::set_const<<<num_blocks, block_size>>>(n, val, arr);
+  }
+
+}} // namespace ReSolve::vector