1#ifndef __COMMON_GRADIENT_JUMP_PENALTY_KERNEL_H__
2#define __COMMON_GRADIENT_JUMP_PENALTY_KERNEL_H__
50 const int jk =
ijk / nx;
51 const int i =
ijk -
jk * nx;
52 const int k =
jk / nx;
53 const int j =
jk -
k * nx;
89 const int jk =
ijk / nx;
90 const int i =
ijk -
jk * nx;
91 const int k =
jk / nx;
92 const int j =
jk -
k * nx;
94 * dphidxi_d[0+
i*nx] + \
96 * dphidxi_d[nx-1+
i*nx] + \
98 * dphidxi_d[0+
j*nx] + \
100 * dphidxi_d[nx-1+
j*nx] + \
102 * dphidxi_d[0+
k*nx] + \
104 * dphidxi_d[nx-1+
k*nx];
__global__ void dirichlet_apply_scalar_kernel(const int *__restrict__ msk, T *__restrict__ x, const T g, const int m)
__global__ void gradient_jump_penalty_finalize_kernel(T *__restrict__ penalty_d, T *__restrict__ penalty_facet_d, T *__restrict__ dphidxi_d, const int nx)
__global__ void pick_facet_value_hex_kernel(T *__restrict__ b, T *__restrict__ a, const int nx)