42#define coef_normal_area_idx(i, j, k, l, lx, nf) \
43 (((i) + (lx) * (((j) - 1) + (lx) * (((k) - 1) + (nf) * (((l) - 1))))) - 1)
50 const int idx2 = idx -1;
51 index[3] =
idx2/(lx * lx * lx) ;
52 index[2] = (
idx2 - (lx*lx*lx)*index[3])/(lx * lx);
53 index[1] = (
idx2 - (lx*lx*lx)*index[3] - (lx*lx) * index[2]) / lx;
54 index[0] = (
idx2 - (lx*lx*lx)*index[3] - (lx*lx) * index[2]) - lx*index[1];
__inline__ __device__ void nonlinear_index(const int idx, const int lx, int *index)
__global__ void dirichlet_apply_scalar_kernel(const int *__restrict__ msk, T *__restrict__ x, const T g, const int m)