histogram.cpp.patch

patch to fix chi squared formula - Nicu Stiurca, 2011-07-28 04:03 am

Download (717 Bytes)

 
src/histogram.cpp (working copy)
1343 1343 
            for( j = 0; j < len; j++ )
1344 1344 
            {
1345 1345 
                double a = h1[j] - h2[j];
1346 
                double b = h1[j] + h2[j];
1346 
                double b = h1[j];
1347 1347 
                if( fabs(b) > FLT_EPSILON )
1348 1348 
                    result += a*a/b;
1349 1349 
            }
......
1429 1429 
            else
1430 1430 
            {
1431 1431 
                double a = v1 - v2;
1432 
                double b = v1 + v2;
1432 
                double b = v1;
1433 1433 
                if( b > FLT_EPSILON )
1434 1434 
                    result += a*a/b;
1435 1435 
            }