+static void check_for_lost_load()
+{
+ double total_init = process::get_total_load_init();
+
+ double total_exit = process::get_total_load_exit();
+ double lost = total_init - total_exit;
+ double lost_ratio = 100.0 * lost / total_init;
+ if (lost_ratio < -opt::load_ratio_threshold)
+ CRITICAL2("Gained load at exit! %g (%g%%) <============",
+ lost, lost_ratio);
+ else if (lost_ratio > opt::load_ratio_threshold)
+ CRITICAL2("Lost load at exit! %g (%g%%) <============",
+ lost, lost_ratio);
+ else
+ DEBUG2("Total load at exit looks good: %g (%g%%)", lost, lost_ratio);
+
+ double total_running = process::get_total_load_running();
+ double running_ratio = 100.0 * total_running / total_init;
+ if (running_ratio < -opt::load_ratio_threshold)
+ CRITICAL2("Negative running load at exit! %g (%g%%) <============",
+ total_running, running_ratio);
+ else if (running_ratio > opt::load_ratio_threshold)
+ CRITICAL2("Remaining running load at exit! %g (%g%%) <============",
+ total_running, running_ratio);
+ else
+ DEBUG2("Running load at exit looks good: %g (%g%%)",
+ total_running, running_ratio);
+}
+
+static void statistics(const std::vector<double>& vec,
+ double* sum, double* avg, double* stddev)
+{
+ using std::tr1::bind;
+ using std::tr1::placeholders::_1;
+
+ unsigned n = vec.size();
+ double vec_sum = std::accumulate(vec.begin(), vec.end(), 0.0);
+ double vec_avg = vec_sum / n;
+
+ if (sum)
+ *sum = vec_sum;
+ if (avg)
+ *avg = vec_avg;
+
+ if (stddev) {
+ std::vector<double> diff(vec);
+ std::transform(diff.begin(), diff.end(), diff.begin(),
+ bind(std::minus<double>(), _1, vec_avg));
+ double epsilon = std::accumulate(diff.begin(), diff.end(), 0.0);
+ double square_sum = std::inner_product(diff.begin(), diff.end(),
+ diff.begin(), 0.0);
+ double variance = (square_sum - (epsilon * epsilon) / n) / n;
+ *stddev = sqrt(variance);
+ }
+}
+
+static void compute_metrics()
+{
+ statistics(comps, &comp_total, &comp_avg, &comp_stddev);
+ statistics(loads, NULL, &load_avg, &load_stddev);
+}
+