CarpetLib: Avoid nans in timer output

Avoid nans in timer output; instead of dividing by zero, set the
result to zero.

darcs-hash:20060925214858-dae7b-0c82bb1a3dc1a3080955b86324487e4a7efc044e.gz

1 files changed, 12 insertions, 4 deletions

diff --git a/Carpet/CarpetLib/src/timestat.cc b/Carpet/CarpetLib/src/timestat.cc
index 5e8b0b599..b126471be 100644
--- a/Carpet/CarpetLib/src/timestat.cc
+++ b/Carpet/CarpetLib/src/timestat.cc
@@ -63,10 +63,18 @@ void timestat::stop (double const b)
 
 ostream& operator<< (ostream& os, const timestat& wt)
 {
-  double const avg = wt.wtime / wt.count;
-  double const stddev = sqrt (max (0.0, wt.wtime2 / wt.count - pow (avg, 2)));
-  double const bavg = wt.bytes / wt.count;
-  double const bstddev = sqrt (max (0.0, wt.bytes2 / wt.count - pow (bavg, 2)));
+  double avg, stddev, bavg, bstddev;
+  if (wt.count == 0.0) {
+    avg = 0.0;
+    stddev = 0.0;
+    bavg = 0.0;
+    bstddev = 0.0;
+  } else {
+    avg = wt.wtime / wt.count;
+    stddev = sqrt (max (0.0, wt.wtime2 / wt.count - pow (avg, 2)));
+    bavg = wt.bytes / wt.count;
+    bstddev = sqrt (max (0.0, wt.bytes2 / wt.count - pow (bavg, 2)));
+  }
   os << "timestat[seconds]:"
      << " cnt: " << wt.count
      << "   time: sum: " << wt.wtime