 	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];
 			}
 		}

	dim3 gridSize((numCells * numLeftFields * numRightFields / 1024) + 1);

	cudaDocontractFieldFieldScalar_kernel<<<gridSize, blockSize>>>(d_left,
			d_right, d_out, numCells, numLeftCells, numRightCells, numPoints);

	cudaMemcpy(h_out, d_out, sizeof(double) * numCells * numLeftCells,
									numRightCells, cudaMemcpyDeviceToHost);

}
/*

				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_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);
