 	if(myID < numCells) {
 		double temp = 0;
 		for (int qp = 0; qp < numPoints; qp++) {
-			temp += d_left[myId + qp*numCells];
+			temp += d_left[myID + qp*numCells];
 		}
-		d_out[myID]
+		d_out[myID]= temp;
 	}
 }
 
 void
-cudaDocontractDataDataScalar(double * h_out,
+cudaDoContractDataDataScalar(double * h_out,
 		double * h_inLeft,
 		double * h_inRight,
 		int numCells,

			out1_c(cl) = cudaOut[cl];
	}

	rst::subtract(&out1_c_l_r[0], &out2_c_l_r[0], out2_c_l_r.size());
	if (rst::vectorNorm(&out1_c_l_r[0], out1_c_l_r.size(), Intrepid::NORM_ONE) > zero) {
		std::cout << "\n\nINCORRECT contractFieldFieldTensor (1): check cuda; "
		<< " diff-1norm = " << rst::vectorNorm(&out1_c_l_r[0], out1_c_l_r.size(), Intrepid::NORM_ONE) << "\n\n";
	}

	std::cout << "cuda speedup of " << elapsedTime_serial/elapsedTime_cuda << std::endl;

	if(myID < numCells) {
		double temp = 0;
		for (int qp = 0; qp < numPoints; qp++) {
			temp += d_left[myID + qp*numCells];
		}
		d_out[myID]=temp;
	}
}

