Index: src/histogram.cpp =================================================================== --- src/histogram.cpp (revision 5743) +++ src/histogram.cpp (working copy) @@ -1343,7 +1343,7 @@ for( j = 0; j < len; j++ ) { double a = h1[j] - h2[j]; - double b = h1[j] + h2[j]; + double b = h1[j]; if( fabs(b) > FLT_EPSILON ) result += a*a/b; } @@ -1429,7 +1429,7 @@ else { double a = v1 - v2; - double b = v1 + v2; + double b = v1; if( b > FLT_EPSILON ) result += a*a/b; }