Changeset 1670

Timestamp:

12/28/09 22:19:00 (15 years ago)

Author:

Max Voronkov

Message:

increased robustness of the line detection algorithm by fixing a transient problem with not-a-numbers appear in the statistics due to round-off error

File:

• trunk/src/STLineFinder.cpp (modified) (2 diffs)

trunk/src/STLineFinder.cpp

@@ -304 +304 @@
 void LFNoiseEstimator::add(float in)
 {
+   if (isnan(in)) {
+       // normally it shouldn't happen
+       return;
+   }
    itsVariances[itsSampleNumber] = in;
 
@@ -584 +588 @@
                 (meanch2-square(meanch));
       linmean=coeff*(Float(cur_channel)-meanch)+mean;
-      linvariance=sqrt(sumf2/Float(box_chan_cntr)-square(mean)-
-                    square(coeff)*(meanch2-square(meanch)));
+      linvariance=sumf2/Float(box_chan_cntr)-square(mean)-
+                    square(coeff)*(meanch2-square(meanch));
+      if (linvariance<0.) {
+          // this shouldn't happen normally, but could be due to round-off error
+          linvariance = 0; 
+      } else {
+          linvariance = sqrt(linvariance);
+      }
   }
   need2recalculate=False;