20 … void ContractX1d(SharedData_Cuda &data, const CeedScalar *U, const CeedScalar *B, CeedScalar *V) {  in ContractX1d()  argument
24   *V = 0.0;  in ContractX1d()
27       *V += B[i + data.t_id_x * P_1D] * data.slice[i];  // Contract x direction  in ContractX1d()
36 …tractTransposeX1d(SharedData_Cuda &data, const CeedScalar *U, const CeedScalar *B, CeedScalar *V) {  in ContractTransposeX1d()  argument
40   *V = 0.0;  in ContractTransposeX1d()
43       *V += B[data.t_id_x + i * P_1D] * data.slice[i];  // Contract x direction  in ContractTransposeX1d()
129 … void ContractX2d(SharedData_Cuda &data, const CeedScalar *U, const CeedScalar *B, CeedScalar *V) {  in ContractX2d()  argument
133   *V = 0.0;  in ContractX2d()
136       *V += B[i + data.t_id_x * P_1D] * data.slice[i + data.t_id_y * T_1D];  // Contract x direction  in ContractX2d()
145 … void ContractY2d(SharedData_Cuda &data, const CeedScalar *U, const CeedScalar *B, CeedScalar *V) {  in ContractY2d()  argument
149   *V = 0.0;  in ContractY2d()
152       *V += B[i + data.t_id_y * P_1D] * data.slice[data.t_id_x + i * T_1D];  // Contract y direction  in ContractY2d()
161 …tractTransposeY2d(SharedData_Cuda &data, const CeedScalar *U, const CeedScalar *B, CeedScalar *V) {  in ContractTransposeY2d()  argument
165   *V = 0.0;  in ContractTransposeY2d()
168       *V += B[data.t_id_y + i * P_1D] * data.slice[data.t_id_x + i * T_1D];  // Contract y direction  in ContractTransposeY2d()
177 …tractTransposeX2d(SharedData_Cuda &data, const CeedScalar *U, const CeedScalar *B, CeedScalar *V) {  in ContractTransposeX2d()  argument
181   *V = 0.0;  in ContractTransposeX2d()
184       *V += B[data.t_id_x + i * P_1D] * data.slice[i + data.t_id_y * T_1D];  // Contract x direction  in ContractTransposeX2d()
193 …ctTransposeAddX2d(SharedData_Cuda &data, const CeedScalar *U, const CeedScalar *B, CeedScalar *V) {  in ContractTransposeAddX2d()  argument
199       *V += B[data.t_id_x + i * P_1D] * data.slice[i + data.t_id_y * T_1D];  // Contract x direction  in ContractTransposeAddX2d()
322 … void ContractX3d(SharedData_Cuda &data, const CeedScalar *U, const CeedScalar *B, CeedScalar *V) {  in ContractX3d()  argument
332     V[k] = 0.0;  in ContractX3d()
335         V[k] += r_B[i] * data.slice[i + data.t_id_y * T_1D];  // Contract x direction  in ContractX3d()
345 … void ContractY3d(SharedData_Cuda &data, const CeedScalar *U, const CeedScalar *B, CeedScalar *V) {  in ContractY3d()  argument
355     V[k] = 0.0;  in ContractY3d()
358         V[k] += r_B[i] * data.slice[data.t_id_x + i * T_1D];  // Contract y direction  in ContractY3d()
368 … void ContractZ3d(SharedData_Cuda &data, const CeedScalar *U, const CeedScalar *B, CeedScalar *V) {  in ContractZ3d()  argument
370     V[k] = 0.0;  in ContractZ3d()
373         V[k] += B[i + k * P_1D] * U[i];  // Contract z direction  in ContractZ3d()
383 …tractTransposeZ3d(SharedData_Cuda &data, const CeedScalar *U, const CeedScalar *B, CeedScalar *V) {  in ContractTransposeZ3d()  argument
385     V[k] = 0.0;  in ContractTransposeZ3d()
388         V[k] += B[k + i * P_1D] * U[i];  // Contract z direction  in ContractTransposeZ3d()
398 …tractTransposeY3d(SharedData_Cuda &data, const CeedScalar *U, const CeedScalar *B, CeedScalar *V) {  in ContractTransposeY3d()  argument
408     V[k] = 0.0;  in ContractTransposeY3d()
411         V[k] += r_B[i] * data.slice[data.t_id_x + i * T_1D];  // Contract y direction  in ContractTransposeY3d()
421 …ctTransposeAddY3d(SharedData_Cuda &data, const CeedScalar *U, const CeedScalar *B, CeedScalar *V) {  in ContractTransposeAddY3d()  argument
433         V[k] += r_B[i] * data.slice[data.t_id_x + i * T_1D];  // Contract y direction  in ContractTransposeAddY3d()
443 …tractTransposeX3d(SharedData_Cuda &data, const CeedScalar *U, const CeedScalar *B, CeedScalar *V) {  in ContractTransposeX3d()  argument
453     V[k] = 0.0;  in ContractTransposeX3d()
456         V[k] += r_B[i] * data.slice[i + data.t_id_y * T_1D];  // Contract x direction  in ContractTransposeX3d()
466 …ctTransposeAddX3d(SharedData_Cuda &data, const CeedScalar *U, const CeedScalar *B, CeedScalar *V) {  in ContractTransposeAddX3d()  argument
478         V[k] += r_B[i] * data.slice[i + data.t_id_y * T_1D];  // Contract x direction  in ContractTransposeAddX3d()