Stan Math Library  2.20.0
reverse mode automatic differentiation
pack.hpp
Go to the documentation of this file.
1 #ifndef STAN_MATH_OPENCL_KERNELS_PACK_HPP
2 #define STAN_MATH_OPENCL_KERNELS_PACK_HPP
3 #ifdef STAN_OPENCL
4 
7 
8 namespace stan {
9 namespace math {
10 namespace opencl_kernels {
11 // \cond
12 static const char* pack_kernel_code = STRINGIFY(
13  // \endcond
30  __kernel void pack(__global double* A, __global double* B,
31  unsigned int rows, unsigned int cols,
32  unsigned int part) {
33  int i = get_global_id(0);
34  int j = get_global_id(1);
35  if (i < rows && j < cols) {
36  // the packed matrices are stored in row major
37  if (part == LOWER) {
38  const int column_offset = j * rows - (j * (j - 1)) / 2;
39  const int row_offset = (i - j);
40  if (j <= i) {
41  A[column_offset + row_offset] = B(i, j);
42  }
43  } else {
44  const int column_offset = j * (j + 1) / 2;
45  if (j >= i) {
46  A[column_offset + i] = B(i, j);
47  }
48  }
49  }
50  }
51  // \cond
52 );
53 // \endcond
54 
59  "pack", {indexing_helpers, pack_kernel_code});
60 
61 } // namespace opencl_kernels
62 } // namespace math
63 } // namespace stan
64 #endif
65 #endif
int rows(const Eigen::Matrix< T, R, C > &m)
Return the number of rows in the specified matrix, vector, or row vector.
Definition: rows.hpp:20
const kernel_cl< out_buffer, in_buffer, int, int, TriangularViewCL > pack("pack", {indexing_helpers, pack_kernel_code})
See the docs for pack() .
#define STRINGIFY(src)
Definition: kernel_cl.hpp:22
static const char * indexing_helpers
Definition: helpers.hpp:14
int cols(const Eigen::Matrix< T, R, C > &m)
Return the number of columns in the specified matrix, vector, or row vector.
Definition: cols.hpp:20
Creates functor for kernels.
Definition: kernel_cl.hpp:201

     [ Stan Home Page ] © 2011–2018, Stan Development Team.