+ int result;
+ process* proc;
+ try {
+ proc = opt::loba_algorithms.new_instance(opt::loba_algo, argc, argv);
+ result = proc->run();
+ loads.push_back(proc->get_load());
+ delete proc;
+ }
+ catch (std::invalid_argument& e) {
+ THROW1(arg_error, 0, "%s", e.what());
+ }
+ return result;
+}
+
+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 compute_load_imbalance()
+{
+ using std::tr1::bind;
+ using std::tr1::placeholders::_1;
+
+ unsigned n = loads.size();
+ load_avg = std::accumulate(loads.begin(), loads.end(), 0.0) / n;
+
+ std::vector<double> diff(loads);
+ std::transform(diff.begin(), diff.end(), diff.begin(),
+ bind(std::minus<double>(), _1, load_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;
+ load_stddev = sqrt(variance);