From b5130a3b35742148848fdea581bfb0baa06ea95e Mon Sep 17 00:00:00 2001
From: Axel Kohlmeyer <akohlmey@gmail.com>
Date: Sun, 13 Nov 2016 18:46:55 -0500
Subject: [PATCH] avoid NaN for variance from average output
---
src/finish.cpp | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/src/finish.cpp b/src/finish.cpp
index 00d7bf9d6b..0d767b42cd 100644
--- a/src/finish.cpp
+++ b/src/finish.cpp
@@ -865,7 +865,7 @@ void mpi_timings(const char *label, Timer *t, enum Timer::ttype tt,
time_cpu = tmp/nprocs*100.0;
// % variance from the average as measure of load imbalance
- if (time > 1.0e-10)
+ if ((time_sq/time - time) > 1.0e-10)
time_sq = sqrt(time_sq/time - time)*100.0;
else
time_sq = 0.0;
@@ -917,7 +917,7 @@ void omp_times(FixOMP *fix, const char *label, enum Timer::ttype which,
time_std /= nthreads;
time_total /= nthreads;
- if (time_avg > 1.0e-10)
+ if ((time_std/time_avg -time_avg) > 1.0e-10)
time_std = sqrt(time_std/time_avg - time_avg)*100.0;
else
time_std = 0.0;
--
GitLab