add.hpp

Go to the documentation of this file.

#ifndef STAN_MATH_OPENCL_KERNELS_ADD_HPP
#define STAN_MATH_OPENCL_KERNELS_ADD_HPP
#ifdef STAN_OPENCL
 
#include <stan/math/opencl/kernel_cl.hpp>
#include <stan/math/opencl/buffer_types.hpp>
#include <stan/math/opencl/matrix_cl_view.hpp>
#include <string>
 
namespace stan {
namespace math {
namespace opencl_kernels {
 
// \cond
static constexpr const char *add_batch_kernel_code = STRINGIFY(
    // \endcond
    __kernel void add_batch(__global double *B, __global double *A,
                            unsigned int rows, unsigned int cols,
                            unsigned int batch_size) {
      const int i = get_global_id(0);
      const int j = get_global_id(1);
      if (i < rows && j < cols) {
        double temp = 0.0;
        for (int k = 0; k < batch_size; k++) {
          temp += A_batch(i, j, k);
        }
        B(i, j) = temp;
      }
    }
    // \cond
);
// \endcond
 
const kernel_cl<out_buffer, in_buffer, int, int, int> add_batch(
    "add_batch", {indexing_helpers, add_batch_kernel_code});
}  // namespace opencl_kernels
}  // namespace math
}  // namespace stan
#endif
#endif