+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);
+}
+