 				cuda_hostLeft(cl, lbf, qp) = in_c_l_p(cl,lbf, qp);
 				omp_hostLeft(cl,lbf, qp) = in_c_l_p(cl,lbf,qp);
 
-				cudaLeft[cl * p * r + lbf * p + qp] = in_c_ls_p(cl,lbf,qp);
+				cudaLeft[cl * p * r + lbf * p + qp] = in_c_l_p(cl,lbf,qp);
 			}
 			//cudaRightColMajor[cl + c*qp] = in_r_c_p(cl,qp);
 			//cudaLeftColMajor[cl + c*qp] = in_l_c_p(cl,qp);

	for (int cl = 0; cl < c; ++cl) {
		for(int lbf = 0; lbf < l; ++lbf) {
			for(int rbf = 0; rbf < r; ++rbf) {
				cudaOut[cl];
				out1_c_l_r(cl,lbf,rbf) = cudaOut[cl * l * r + lbf * r + rbf];
			}
		}


__global__
void
cudaDocontractFieldFieldScalar_kernel(double * d_left, double * d_right,
double * d_out,
int numCells,
int numLeftFields,
