54   for (CeedInt k = i; k < BASIS_Q_1D * BASIS_P_1D; k += blockDim.x) {  in InterpAtPoints()  local
55     s_chebyshev_interp_1d[k] = chebyshev_interp_1d[k];  in InterpAtPoints()
84         for (CeedInt k = i; k < writeLen; k += blockDim.x) {  in InterpAtPoints()  local
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()
91           out[k] = v_k;  in InterpAtPoints()
137   for (CeedInt k = i; k < BASIS_Q_1D * BASIS_P_1D; k += blockDim.x) {  in InterpTransposeAtPoints()  local
138     s_chebyshev_interp_1d[k] = chebyshev_interp_1d[k];  in InterpTransposeAtPoints()
158       for (CeedInt k = i; k < BASIS_NUM_QPTS; k += blockDim.x) {  in InterpTransposeAtPoints()  local
159         s_chebyshev_coeffs[k] = 0.0;  in InterpTransposeAtPoints()
203         for (CeedInt k = i; k < writeLen; k += blockDim.x) {  in InterpTransposeAtPoints()  local
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()
210           if (d == BASIS_DIM - 1) out[k] += v_k;  in InterpTransposeAtPoints()
211           else out[k] = v_k;  in InterpTransposeAtPoints()
233   for (CeedInt k = i; k < BASIS_Q_1D * BASIS_P_1D; k += blockDim.x) {  in GradAtPoints()  local
234     s_chebyshev_interp_1d[k] = chebyshev_interp_1d[k];  in GradAtPoints()
264         for (CeedInt k = i; k < writeLen; k += blockDim.x) {  in GradAtPoints()  local
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()
271           out[k] = v_k;  in GradAtPoints()
322   for (CeedInt k = i; k < BASIS_Q_1D * BASIS_P_1D; k += blockDim.x) {  in GradTransposeAtPoints()  local
323     s_chebyshev_interp_1d[k] = chebyshev_interp_1d[k];  in GradTransposeAtPoints()
344       for (CeedInt k = i; k < BASIS_NUM_QPTS; k += blockDim.x) {  in GradTransposeAtPoints()  local
345         s_chebyshev_coeffs[k] = 0.0;  in GradTransposeAtPoints()
394         for (CeedInt k = i; k < writeLen; k += blockDim.x) {  in GradTransposeAtPoints()  local
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()
401           if (d == BASIS_DIM - 1) out[k] += v_k;  in GradTransposeAtPoints()
402           else out[k] = v_k;  in GradTransposeAtPoints()