         }
       }
     }
-
+    /*
     Kokkos::deep_copy(kokkosLeft, hostLeft);
     Kokkos::deep_copy(kokkosRight, hostRight);
     Kokkos::deep_copy(kokkosOut, hostOut);
-
+    
     ContractFieldFieldTensorFunctor<input_view_t, input_view_t, output_view_t>
       kokkosFunctor(kokkosLeft, kokkosRight, kokkosOut,
 		    numLeftFields, numRightFields, numPoints,
 		    dim1Tensor, dim2Tensor);
-    //Kokkos::parallel_for(numCells, kokkosFunctor);
-
+    Kokkos::parallel_for(numCells, kokkosFunctor);
+    
     Kokkos::fence();
 
     Kokkos::deep_copy(hostOut, kokkosOut);

   */
   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_BLAS; "
		<< " diff-1norm = " << rst::vectorNorm(&out1_c_l_r[0], out1_c_l_r.size(), Intrepid::NORM_ONE) << "\n\n";
   }
   else {
     std::cout << "Cpp and blas get same result" << std::endl;
   }


  contractFieldFieldTensor(out1_c_l_r, in_c_l_p_d_d, in_c_r_p_d_d, 0);
  contractFieldFieldTensor(out2_c_l_r, in_c_l_p_d_d, in_c_r_p_d_d, 1);

  /*
  Kokkos::initialize();



 contractFieldFieldTensor(out1_c_l_r, in_c_l_p_d_d, in_c_r_p_d_d, 0);
 contractFieldFieldTensor(out2_c_l_r, in_c_l_p_d_d, in_c_r_p_d_d, 2);
 

  /* I got rid of the typedefs - for now
   * 0 -> manual computation
