 
 	std::cout << "trying cuda" << std::endl;
 	//Now try the cuda version, start with warmup
-	cudaDoContractFieldFieldTensor(cudaOut,cudaLeft,cudaRight, c, p, 1);
+	cudaDoContractDataDataScalar(cudaOut,cudaLeftColMajor,cudaRightColMajor, c, p, 1);
 
 	clock_gettime(CLOCK_MONOTONIC, &tic);
 	for(int i = 0; i < 5; ++i){
-		cudaDoContractFieldFieldTensor(cudaOut,cudaLeft,cudaRight, c, p, 1);
+		cudaDoContractDataDataScalar(cudaOut,cudaLeftColMajor,cudaRightColMajor, c, p, 1);
 	}
 
 	clock_gettime(CLOCK_MONOTONIC, &toc);

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

void
cudaDocontractDataDataScalar(double * h_out,
		double * h_inLeft,
		double * h_inRight,
		int numCells,

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

