146 template <int NUM_COMP_OUT, int NUM_COMP_FIELD, int Q_1D>
149   if (data.t_id_x < Q_1D) {  in WriteLVecStandard1d_QFAssembly()
150     const CeedInt ind = data.t_id_x + elem * Q_1D;  in WriteLVecStandard1d_QFAssembly()
153 …   d_v[ind + (input_offset * NUM_COMP_OUT + output_offset + comp) * (Q_1D * num_elem)] = r_v[comp];  in WriteLVecStandard1d_QFAssembly()
274 template <int NUM_COMP_OUT, int NUM_COMP_FIELD, int Q_1D>
277   if (data.t_id_x < Q_1D && data.t_id_y < Q_1D) {  in WriteLVecStandard2d_QFAssembly()
278     const CeedInt ind = (data.t_id_x + data.t_id_y * Q_1D) + elem * Q_1D * Q_1D;  in WriteLVecStandard2d_QFAssembly()
281 …d_v[ind + (input_offset * NUM_COMP_OUT + output_offset + comp) * (Q_1D * Q_1D * num_elem)] = r_v[c…  in WriteLVecStandard2d_QFAssembly()
353 template <int NUM_COMP, int COMP_STRIDE, int Q_1D>
357   if (data.t_id_x < Q_1D && data.t_id_y < Q_1D) {  in ReadEVecSliceStandard3d()
358     const CeedInt node = data.t_id_x + data.t_id_y * Q_1D + q * Q_1D * Q_1D;  in ReadEVecSliceStandard3d()
359     const CeedInt ind  = indices[node + elem * Q_1D * Q_1D * Q_1D];  in ReadEVecSliceStandard3d()
368 template <int NUM_COMP, int Q_1D, int STRIDES_NODE, int STRIDES_COMP, int STRIDES_ELEM>
371   if (data.t_id_x < Q_1D && data.t_id_y < Q_1D) {  in ReadEVecSliceStrided3d()
372     const CeedInt node = data.t_id_x + data.t_id_y * Q_1D + q * Q_1D * Q_1D;  in ReadEVecSliceStrided3d()
442 template <int NUM_COMP_OUT, int NUM_COMP_FIELD, int Q_1D>
445   if (data.t_id_x < Q_1D && data.t_id_y < Q_1D) {  in WriteLVecStandard3d_QFAssembly()
446     for (CeedInt z = 0; z < Q_1D; z++) {  in WriteLVecStandard3d_QFAssembly()
447 …const CeedInt ind = (data.t_id_x + data.t_id_y * Q_1D + z * Q_1D * Q_1D) + elem * Q_1D * Q_1D * Q_…  in WriteLVecStandard3d_QFAssembly()
450 …_offset * NUM_COMP_OUT + output_offset + comp) * (Q_1D * Q_1D * Q_1D * num_elem)] = r_v[z + comp *…  in WriteLVecStandard3d_QFAssembly()
475 template <int NUM_COMP, int Q_1D, int T_1D>
478   if (data.t_id_x < Q_1D && data.t_id_y < Q_1D) {  in GradColloSlice3d()
481       data.slice[data.t_id_x + data.t_id_y * T_1D] = r_U[q + comp * Q_1D];  in GradColloSlice3d()
485       for (CeedInt i = 0; i < Q_1D; i++) {  in GradColloSlice3d()
486 …      r_V[comp + 0 * NUM_COMP] += c_G[i + data.t_id_x * Q_1D] * data.slice[i + data.t_id_y * T_1D];  in GradColloSlice3d()
490       for (CeedInt i = 0; i < Q_1D; i++) {  in GradColloSlice3d()
491 …      r_V[comp + 1 * NUM_COMP] += c_G[i + data.t_id_y * Q_1D] * data.slice[data.t_id_x + i * T_1D];  in GradColloSlice3d()
495       for (CeedInt i = 0; i < Q_1D; i++) {  in GradColloSlice3d()
496         r_V[comp + 2 * NUM_COMP] += c_G[i + q * Q_1D] * r_U[i + comp * Q_1D];  in GradColloSlice3d()
505 template <int NUM_COMP, int Q_1D, int T_1D>
508   if (data.t_id_x < Q_1D && data.t_id_y < Q_1D) {  in GradColloSliceTranspose3d()
514       for (CeedInt i = 0; i < Q_1D; i++) {  in GradColloSliceTranspose3d()
515         r_V[q + comp * Q_1D] += c_G[data.t_id_x + i * Q_1D] * data.slice[i + data.t_id_y * T_1D];  in GradColloSliceTranspose3d()
521       for (CeedInt i = 0; i < Q_1D; i++) {  in GradColloSliceTranspose3d()
522         r_V[q + comp * Q_1D] += c_G[data.t_id_y + i * Q_1D] * data.slice[data.t_id_x + i * T_1D];  in GradColloSliceTranspose3d()
525       for (CeedInt i = 0; i < Q_1D; i++) {  in GradColloSliceTranspose3d()
526         r_V[i + comp * Q_1D] += c_G[i + q * Q_1D] * r_U[comp + 2 * NUM_COMP];  in GradColloSliceTranspose3d()