 
 
 int main(int argc, char* argv[]) {
-	int c=5, p=9, l=3, r=7;
+	int c=550000, p=9, l=3, r=7;
 
 	FieldContainer<double> in_c_l_p(c, l, p);
 	FieldContainer<double> in_c_r_p(c, r, p);
 	FieldContainer<double> out1_c_l_r(c, l, r);
-	FieldContainer<double> out2_c_l_r_l_r_l_r(c, l, r);
-	double zero = INTREPID_TOL*10000.0;
+	FieldContainer<double> out2_c_l_r(c, l, r);
+	double zero = Intrepid::INTREPID_TOL*10000.0;
 
 	// fill with random numbers
 	for (int i=0; i<in_c_l_p.size(); i++) {

		double *                                  time = 0) {

	// get sizes
	int numCells        = leftFields.dimension(0);
	int numLeftFields   = leftFields.dimension(1);
	int numRightFields  = rightFields.dimension(1);
	int numPoints       = leftFields.dimension(2);

	// Deep copy Kokkos host views into device views
	Kokkos::deep_copy(leftDevice, leftHost);

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

}
*/
template<class DeviceType, class LeftViewType, class RightViewType, class OutputViewType>
struct contractFieldFieldScalarFunctor {
	typedef DeviceType device_type;
