网上流传的一个Canny边缘检测算法C语言实现的bug

主程序部分
// 对原图象进行滤波
GaussianSmooth(pUnchImage, nWidth, nHeight, sigma, pUnchSmooth) ;
// 计算方向导数
DirGrad(pUnchSmooth, nWidth, nHeight, pnGradX, pnGradY) ;
// 计算梯度的幅度
GradMagnitude(pnGradX, pnGradY, nWidth, nHeight, pnGradMag) ;
// 应用non-maximum 抑制
NonmaxSuppress(pnGradMag, pnGradX, pnGradY, nWidth, nHeight, pUnchEdge) ;
// 应用Hysteresis，找到所有的边界
Hysteresis(pnGradMag, nWidth, nHeight, dRatioLow, dRatioHigh, pUnchEdge);

这段代码在很多地方都看到过，最近在dsp上用到了这段代码，自己花了点事件优化了一下代码，然后发现其中的一个bug，
void MakeGauss(double sigma, double **pdKernel, int *pnWindowSize)
这个函数中有一句
dValue = exp(-(1/2)*dDis*dDis/(sigma*sigma)) / (sqrt(2 * PI) * sigma );
看出有什么问题了么？如果你在下一句设断点，然后调试一下，会发现每次dValue的值都一样，因为(1/2)就结果是0 。
解决办法是用-0.5代替-(1/2)就行了