 	std::cout << "trying serial" << std::endl;
 
 	//Warmup
-	contractFieldFieldScalarSerial(out2_c_l_r, in_l_c_p, in_r_c_p);
+	contractFieldFieldScalarSerial(out2_c_l_r, in_c_l_p, in_c_r_p);
 
 	timespec tic;
 	clock_gettime(CLOCK_MONOTONIC, &tic);
 
 	//repeat the calculation 5 times so we can average out some randomness
 	for(int i = 0; i < 5; ++i){
-		contractFieldFieldScalarSerial(out2_c_l_r, in_l_c_p, in_r_c_p);
+		contractFieldFieldScalarSerial(out2_c_l_r, in_c_l_p, in_c_r_p);
 	}
 
 	timespec toc;

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

	}
}
*/
/*
void
cudaDocontractFieldFieldScalarUnrolled(double * h_out,
double * h_inLeft,
