147 template <int NUM_COMP_OUT, int NUM_COMP_FIELD, int Q_1D>
150   if (data.t_id_x < Q_1D) {  in WriteLVecStandard1d_QFAssembly()
151     const CeedInt ind = data.t_id_x + elem * Q_1D;  in WriteLVecStandard1d_QFAssembly()
154 …   d_v[ind + (input_offset * NUM_COMP_OUT + output_offset + comp) * (Q_1D * num_elem)] = r_v[comp];  in WriteLVecStandard1d_QFAssembly()
276 template <int NUM_COMP_OUT, int NUM_COMP_FIELD, int Q_1D>
279   if (data.t_id_x < Q_1D && data.t_id_y < Q_1D) {  in WriteLVecStandard2d_QFAssembly()
280     const CeedInt ind = (data.t_id_x + data.t_id_y * Q_1D) + elem * Q_1D * Q_1D;  in WriteLVecStandard2d_QFAssembly()
283 …d_v[ind + (input_offset * NUM_COMP_OUT + output_offset + comp) * (Q_1D * Q_1D * num_elem)] = r_v[c…  in WriteLVecStandard2d_QFAssembly()
356 template <int NUM_COMP, int COMP_STRIDE, int Q_1D>
360   if (data.t_id_x < Q_1D && data.t_id_y < Q_1D) {  in ReadEVecSliceStandard3d()
361     const CeedInt node = data.t_id_x + data.t_id_y * Q_1D + q * Q_1D * Q_1D;  in ReadEVecSliceStandard3d()
362     const CeedInt ind  = indices[node + elem * Q_1D * Q_1D * Q_1D];  in ReadEVecSliceStandard3d()
371 template <int NUM_COMP, int Q_1D, int STRIDES_NODE, int STRIDES_COMP, int STRIDES_ELEM>
374   if (data.t_id_x < Q_1D && data.t_id_y < Q_1D) {  in ReadEVecSliceStrided3d()
375     const CeedInt node = data.t_id_x + data.t_id_y * Q_1D + q * Q_1D * Q_1D;  in ReadEVecSliceStrided3d()
445 template <int NUM_COMP_OUT, int NUM_COMP_FIELD, int Q_1D>
448   if (data.t_id_x < Q_1D && data.t_id_y < Q_1D) {  in WriteLVecStandard3d_QFAssembly()
449     for (CeedInt z = 0; z < Q_1D; z++) {  in WriteLVecStandard3d_QFAssembly()
450 …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()
453 …_offset * NUM_COMP_OUT + output_offset + comp) * (Q_1D * Q_1D * Q_1D * num_elem)] = r_v[z + comp *…  in WriteLVecStandard3d_QFAssembly()
478 template <int NUM_COMP, int Q_1D, int T_1D>
481   if (data.t_id_x < Q_1D && data.t_id_y < Q_1D) {  in GradColloSlice3d()
484       data.slice[data.t_id_x + data.t_id_y * T_1D] = r_U[q + comp * Q_1D];  in GradColloSlice3d()
488       for (CeedInt i = 0; i < Q_1D; i++) {  in GradColloSlice3d()
489 …      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()
493       for (CeedInt i = 0; i < Q_1D; i++) {  in GradColloSlice3d()
494 …      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()
498       for (CeedInt i = 0; i < Q_1D; i++) {  in GradColloSlice3d()
499         r_V[comp + 2 * NUM_COMP] += c_G[i + q * Q_1D] * r_U[i + comp * Q_1D];  in GradColloSlice3d()
508 template <int NUM_COMP, int Q_1D, int T_1D>
511   if (data.t_id_x < Q_1D && data.t_id_y < Q_1D) {  in GradColloSliceTranspose3d()
517       for (CeedInt i = 0; i < Q_1D; i++) {  in GradColloSliceTranspose3d()
518         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()
524       for (CeedInt i = 0; i < Q_1D; i++) {  in GradColloSliceTranspose3d()
525         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()
528       for (CeedInt i = 0; i < Q_1D; i++) {  in GradColloSliceTranspose3d()
529         r_V[i + comp * Q_1D] += c_G[i + q * Q_1D] * r_U[comp + 2 * NUM_COMP];  in GradColloSliceTranspose3d()