72       CeedInt           post  = 1;  in InterpAtPoints()  local
81         const CeedInt     writeLen = pre * post * Q;  in InterpAtPoints()
85           const CeedInt c   = k % post;  in InterpAtPoints()
86           const CeedInt j   = (k / post) % Q;  in InterpAtPoints()
87           const CeedInt a   = k / (post * Q);  in InterpAtPoints()
90 …t b = 0; b < P; b++) v_k += s_chebyshev_interp_1d[j * BASIS_P_1D + b] * in[(a * P + b) * post + c];  in InterpAtPoints()
93         post *= Q;  in InterpAtPoints()
100         post = 1;  in InterpAtPoints()
112             for (CeedInt c = 0; c < post; c++) {  in InterpAtPoints()
115               for (CeedInt b = 0; b < Q; b++) v_k += chebyshev_x[b] * in[(a * Q + b) * post + c];  in InterpAtPoints()
116               out[a * post + c] = v_k;  in InterpAtPoints()
119           post *= 1;  in InterpAtPoints()
155       CeedInt           post  = 1;  in InterpTransposeAtPoints()  local
167         post = 1;  in InterpTransposeAtPoints()
179             for (CeedInt c = 0; c < post; c++) {  in InterpTransposeAtPoints()
181 … j < Q; j++) atomicAdd(&out[(a * Q + (j + p) % Q) * post + c], chebyshev_x[(j + p) % Q] * in[a * p…  in InterpTransposeAtPoints()
183 …   for (CeedInt j = 0; j < Q; j++) out[(a * Q + j) * post + c] = chebyshev_x[j] * in[a * post + c];  in InterpTransposeAtPoints()
187           post *= Q;  in InterpTransposeAtPoints()
193       post = 1;  in InterpTransposeAtPoints()
200         const CeedInt     writeLen = pre * post * P;  in InterpTransposeAtPoints()
204           const CeedInt c   = k % post;  in InterpTransposeAtPoints()
205           const CeedInt j   = (k / post) % P;  in InterpTransposeAtPoints()
206           const CeedInt a   = k / (post * P);  in InterpTransposeAtPoints()
209 …t b = 0; b < Q; b++) v_k += s_chebyshev_interp_1d[j + b * BASIS_P_1D] * in[(a * Q + b) * post + c];  in InterpTransposeAtPoints()
213         post *= P;  in InterpTransposeAtPoints()
252       CeedInt           post  = 1;  in GradAtPoints()  local
261         const CeedInt     writeLen = pre * post * Q;  in GradAtPoints()
265           const CeedInt c   = k % post;  in GradAtPoints()
266           const CeedInt j   = (k / post) % Q;  in GradAtPoints()
267           const CeedInt a   = k / (post * Q);  in GradAtPoints()
270 …t b = 0; b < P; b++) v_k += s_chebyshev_interp_1d[j * BASIS_P_1D + b] * in[(a * P + b) * post + c];  in GradAtPoints()
273         post *= Q;  in GradAtPoints()
283           post = 1;  in GradAtPoints()
296               for (CeedInt c = 0; c < post; c++) {  in GradAtPoints()
299                 for (CeedInt b = 0; b < Q; b++) v_k += chebyshev_x[b] * in[(a * Q + b) * post + c];  in GradAtPoints()
300                 out[a * post + c] = v_k;  in GradAtPoints()
303             post *= 1;  in GradAtPoints()
341       CeedInt     post  = 1;  in GradTransposeAtPoints()  local
356           post = 1;  in GradTransposeAtPoints()
369               for (CeedInt c = 0; c < post; c++) {  in GradTransposeAtPoints()
371 … j < Q; j++) atomicAdd(&out[(a * Q + (j + p) % Q) * post + c], chebyshev_x[(j + p) % Q] * in[a * p…  in GradTransposeAtPoints()
373 …   for (CeedInt j = 0; j < Q; j++) out[(a * Q + j) * post + c] = chebyshev_x[j] * in[a * post + c];  in GradTransposeAtPoints()
377             post *= Q;  in GradTransposeAtPoints()
384       post = 1;  in GradTransposeAtPoints()
391         const CeedInt     writeLen = pre * post * P;  in GradTransposeAtPoints()
395           const CeedInt c   = k % post;  in GradTransposeAtPoints()
396           const CeedInt j   = (k / post) % P;  in GradTransposeAtPoints()
397           const CeedInt a   = k / (post * P);  in GradTransposeAtPoints()
400 …t b = 0; b < Q; b++) v_k += s_chebyshev_interp_1d[j + b * BASIS_P_1D] * in[(a * Q + b) * post + c];  in GradTransposeAtPoints()
404         post *= P;  in GradTransposeAtPoints()