 
 	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 COMP_CPP vs. COMP_KOKKOS; "
+		std::cout << "\n\nINCORRECT contractFieldFieldTensor (0): check COMP_CPP vs. COMP_KOKKOS; "
 			<< " diff-1norm = " << rst::vectorNorm(&out1_c_l_r[0], out1_c_l_r.size(), Intrepid::NORM_ONE) << "\n\n";
 	}
 

		int numCells,
		int numLeftFields,
		int numRightFields,
		int numPoints) {

	double * d_right;
	double * d_left;


		double temp = 0;
		for (int qp = 0; qp < numPoints; qp++) {
			temp += d_left[myMatrix*numPoints*numLeftFields + numPoints*matrixRow + qp] *
							d_right[myMatrix*numPoints*numRightFields + qp*numRightFields + matrixCol];
		}
		d_out[myID]= temp;
