Utils: Fix NAN error in Statistics

[senf.git] / senf / Utils / Statistics.cc

diff --git a/senf/Utils/Statistics.cc b/senf/Utils/Statistics.cc
index 6e6eb18..445e1c4 100644 (file)
--- a/senf/Utils/Statistics.cc
+++ b/senf/Utils/Statistics.cc
@@ -30,6 +30,7 @@
 #include <cmath>
 #include <sstream>
 #include <senf/Utils/Console/Console.hh>
+#include <senf/Utils/Format.hh>
 #include "StatisticsTargets.hh"
 
 //#include "Statistics.mpp"
@@ -82,63 +83,23 @@ senf::StatisticsBase::output(unsigned n)
     return OutputProxy<StatisticsBase>(this, &(i->second));
 }
 
-//
-// generate an engineering style notation 
-//
-std::string format_eng( float f)
-{
-    char buf[16];
-    if (f > 0) {
-        int n = 0;
-        while( f >= 1000.0f) {
-                f /= 1000.0f;
-                n+=3;
-        }
-
-        if( n >=3)
-            sprintf( buf, " %3.2fe%+03d", f, n);
-        else
-            sprintf( buf, "      %3.2f", f);
-    }
-    else if (f < 0) {
-        int n = 0;
-        while( f <= -1000.0f) {
-                f *= 1000.0f;
-                n+=3;
-        }
-        if( n >=3)
-            sprintf( buf, " %3.2fe%+03d", f, n);
-        else
-            sprintf( buf, "      %3.2f", f);
-    }
-    else{
-        sprintf( buf, "        0.00");
-    }
-
-    return buf;
-}
-
-
 prefix_ void senf::StatisticsBase::consoleList(unsigned level, std::ostream & os)
     const
 {
-    os << boost::format("%s%-5d%|15t|  %12s  %12s+/-%12s  %12s\n") 
-        % std::string(2*level,' ') 
-        % rank() 
-        % format_eng(min()) 
-        % format_eng(avg()) 
-        % format_eng(dev())
-        % format_eng(max());
+    namespace fmt = senf::format;
+
+    os << boost::format("%s%-5d%|15t|  %12.5g  %19.5g  %12.5g\n") 
+        % std::string(2*level,' ') % rank() 
+        % fmt::eng(min()).setw() % fmt::eng(avg(),dev()).setw() % fmt::eng(max()).setw();
     {
         OutputMap::const_iterator i (outputs_.begin());
         OutputMap::const_iterator i_end (outputs_.end());
         for (; i != i_end; ++i)
-            os << boost::format("            %3d  %12s  %12s+/-%12s  %12s\n")
+            os << boost::format("            %3d  %12.5g  %19.5g  %12.5g\n")
                 % i->second.n 
-                % format_eng(i->second.min/i->second.n) 
-                % format_eng(i->second.avg/i->second.n)
-                % format_eng(i->second.dev/i->second.n)
-                % format_eng(i->second.max/i->second.n);
+                % fmt::eng(i->second.min).setw() 
+                % fmt::eng(i->second.avg, i->second.dev).setw() 
+                % fmt::eng(i->second.max).setw();
     }
     {
         Children::const_iterator i (children_.begin());
@@ -167,7 +128,11 @@ prefix_ void senf::StatisticsBase::generateOutput()
             i->second.max += j->max;
             i->second.dev += j->dev;
         }
-        i->second.signal(i->second.min/n, i->second.avg/n, i->second.max/n, i->second.dev/n);
+        i->second.min /= n;
+        i->second.avg /= n;
+        i->second.max /= n;
+        i->second.dev /= n;
+        i->second.signal(i->second.min, i->second.avg, i->second.max, i->second.dev);
     }
 }
 
@@ -245,7 +210,7 @@ prefix_ senf::Statistics::Statistics()
 
 prefix_ void senf::Statistics::consoleList(std::ostream & os)
 {
-    os << "RANK        WIN           MIN           AVG            DEV           MAX\n";
+    os << "RANK        WIN       MIN          AVG                   MAX\n";
     StatisticsBase::consoleList(0, os);
 }
 
@@ -305,11 +270,9 @@ prefix_ void senf::Collector::enter(float min, float avg, float max, float dev)
     accSumSq_ += avg*avg + dev*dev;
     if (min < accMin_) accMin_ = min;
     if (max > accMax_) accMax_ = max;
-//    std::cerr << "! " << i_ << ' ' << avg << ' ' << dev << ' ' << accSum_ << ' ' << accSumSq_ 
-//              << std::endl;
     if (++i_ >= rank_) {
         float accAvg (accSum_ / i_);
-        float accDev (std::sqrt(accSumSq_ / i_ - accAvg*accAvg));
+        float accDev (std::max(0.0f,std::sqrt(accSumSq_ / i_ - accAvg*accAvg)));
         StatisticsBase::enter(accMin_, accAvg, accMax_, accDev);
         i_ = 0;
         accMin_ = FLT_MAX;