 
 	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(colMajor)
		cudaDocontractDataDataScalar_kernelColMajor<<<gridSize, blockSize>>>(d_left,
			d_right, d_out, numCells,numPoints,);
	else
		cudaDocontractDataDataScalar_kernelRowMajor<<<gridSize, blockSize>>>(d_left,
		d_right, d_out, numCells,numPoints,);
	cudaMemcpy(h_out, d_out, sizeof(double) * numCells, cudaMemcpyDeviceToHost);

}

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

