 
 __global__
 void
-cudaDocontractFieldFieldScalar_kernelRowMajor(double * d_left, double * d_right,
+cudaDocontractFieldFieldScalar_kernel(double * d_left, double * d_right,
 double * d_out,
 int numCells,
 int numLeftFields,

	dim3 gridSize((numCells * numLeftFields * numRightFields / 1024) + 1);

	cudaDocontractFieldFieldScalar_kernel<<<gridSize, blockSize>>>(d_left,
			d_right, d_out, numCells, numLeftCells, numRightCells, numPoints);

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

}
/*

		int matrixRow = matrixIndex / numRightFields;
		int matrixCol = matrixIndex % numRightFields;

		double temp = 0;
		for (int qp = 0; qp < numPoints; qp++) {
			temp += d_left[myMatrix*numPoints*numLeftFields + matrixRow*numPoints + qp] *
