Fixed LoG offset

src/c/Cuda/CudaUtilities.cu

@@ -193,17 +193,17 @@ Vec<size_t> createLoGKernel(Vec<float> sigma, float** kernelOut, size_t& kernSiz
 		Vec<int> curPos(0, 0, 0);
 		for(curPos.z = 0; curPos.z<kernelDims.z||curPos.z==0; ++curPos.z)
 		{
-			double zMuSub = SQR(curPos.z-center.z+0.5);
+			double zMuSub = SQR(curPos.z-center.z);
 			double zMuSubSigSqr = zMuSub/(2*sigmaSqr.z);
 			double zAdditive = zMuSub/sigma4th.z-1.0/sigmaSqr.z;
 			for(curPos.y = 0; curPos.y<kernelDims.y||curPos.y==0; ++curPos.y)
 			{
-				double yMuSub = SQR(curPos.y-center.y+0.5);
+				double yMuSub = SQR(curPos.y-center.y);
 				double yMuSubSigSqr = yMuSub/(2*sigmaSqr.y);
 				double yAdditive = yMuSub/sigma4th.y-1.0/sigmaSqr.y;
 				for(curPos.x = 0; curPos.x<kernelDims.x||curPos.x==0; ++curPos.x)
 				{
-					double xMuSub = SQR(curPos.x-center.x+0.5);
+					double xMuSub = SQR(curPos.x-center.x);
 					double xMuSubSigSqr = xMuSub/(2*sigmaSqr.x);
 					double xAdditive = xMuSub/sigma4th.x-1.0/sigmaSqr.x;
 					
@@ -233,19 +233,19 @@ Vec<size_t> createLoGKernel(Vec<float> sigma, float** kernelOut, size_t& kernSiz
 		{
 			if(sigma.z!=0)
 			{
-				gaussComp.z = SQR(curPos.z-center.z+0.5)/twoSigmaSqr.z;
+				gaussComp.z = SQR(curPos.z-center.z)/twoSigmaSqr.z;
 			}
 			for(curPos.y = 0; curPos.y<kernelDims.y || curPos.y==0; ++curPos.y)
 			{
 				if(sigma.y!=0)
 				{
-					gaussComp.y = SQR(curPos.y-center.y+0.5)/twoSigmaSqr.y;
+					gaussComp.y = SQR(curPos.y-center.y)/twoSigmaSqr.y;
 				}
 				for(curPos.x =0; curPos.x<kernelDims.x || curPos.x==0; ++curPos.x)
 				{
 					if(sigma.x!=0)
 					{
-						gaussComp.x = SQR(curPos.x-center.x+0.5)/twoSigmaSqr.x;
+						gaussComp.x = SQR(curPos.x-center.x)/twoSigmaSqr.x;
 					}
 					double gauss = gaussComp.sum();
 					kernel[kernelDims.linearAddressAt(curPos)] = ((1.0-gauss)*exp(-gauss))/denominator;