pack.hpp

Go to the documentation of this file.

#ifndef STAN_MATH_OPENCL_KERNELS_PACK_HPP
#define STAN_MATH_OPENCL_KERNELS_PACK_HPP
#ifdef STAN_OPENCL
 
#include <stan/math/opencl/kernel_cl.hpp>
#include <stan/math/opencl/buffer_types.hpp>
#include <string>
 
namespace stan {
namespace math {
namespace opencl_kernels {
// \cond
static constexpr const char* pack_kernel_code = STRINGIFY(
    // \endcond
    __kernel void pack(__global double* A, __global double* B,
                       unsigned int rows, unsigned int cols,
                       unsigned int view) {
      int i = get_global_id(0);
      int j = get_global_id(1);
      if (i < rows && j < cols) {
        // the packed matrices are stored in row major
        if (view == LOWER) {
          const int column_offset = j * rows - (j * (j - 1)) / 2;
          const int row_offset = (i - j);
          if (j <= i) {
            A[column_offset + row_offset] = B(i, j);
          }
        } else {
          const int column_offset = j * (j + 1) / 2;
          if (j >= i) {
            A[column_offset + i] = B(i, j);
          }
        }
      }
    }
    // \cond
);
// \endcond
 
const kernel_cl<out_buffer, in_buffer, int, int, matrix_cl_view> pack(
    "pack", {indexing_helpers, pack_kernel_code});
 
}  // namespace opencl_kernels
}  // namespace math
}  // namespace stan
#endif
#endif