22   data.slice[data.t_id_x] = *U;  in ContractX1d()
27       *V += B[i + data.t_id_x * P_1D] * data.slice[i];  // Contract x direction  in ContractX1d()
38   data.slice[data.t_id_x] = *U;  in ContractTransposeX1d()
43       *V += B[data.t_id_x + i * P_1D] * data.slice[i];  // Contract x direction  in ContractTransposeX1d()
131   data.slice[data.t_id_x + data.t_id_y * T_1D] = *U;  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()
147   data.slice[data.t_id_x + data.t_id_y * T_1D] = *U;  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()
163   data.slice[data.t_id_x + data.t_id_y * T_1D] = *U;  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()
179   data.slice[data.t_id_x + data.t_id_y * T_1D] = *U;  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()
195   data.slice[data.t_id_x + data.t_id_y * T_1D] = *U;  in ContractTransposeAddX2d()
199       *V += B[data.t_id_x + i * P_1D] * data.slice[i + data.t_id_y * T_1D];  // Contract x direction  in ContractTransposeAddX2d()
330     data.slice[data.t_id_x + data.t_id_y * T_1D] = U[k];  in ContractX3d()
335         V[k] += r_B[i] * data.slice[i + data.t_id_y * T_1D];  // Contract x direction  in ContractX3d()
353     data.slice[data.t_id_x + data.t_id_y * T_1D] = U[k];  in ContractY3d()
358         V[k] += r_B[i] * data.slice[data.t_id_x + i * T_1D];  // Contract y direction  in ContractY3d()
406     data.slice[data.t_id_x + data.t_id_y * T_1D] = U[k];  in ContractTransposeY3d()
411         V[k] += r_B[i] * data.slice[data.t_id_x + i * T_1D];  // Contract y direction  in ContractTransposeY3d()
429     data.slice[data.t_id_x + data.t_id_y * T_1D] = U[k];  in ContractTransposeAddY3d()
433         V[k] += r_B[i] * data.slice[data.t_id_x + i * T_1D];  // Contract y direction  in ContractTransposeAddY3d()
451     data.slice[data.t_id_x + data.t_id_y * T_1D] = U[k];  in ContractTransposeX3d()
456         V[k] += r_B[i] * data.slice[i + data.t_id_y * T_1D];  // Contract x direction  in ContractTransposeX3d()
474     data.slice[data.t_id_x + data.t_id_y * T_1D] = U[k];  in ContractTransposeAddX3d()
478         V[k] += r_B[i] * data.slice[i + data.t_id_y * T_1D];  // Contract x direction  in ContractTransposeAddX3d()