std::accumulate(loads.begin(), loads.end(), 0.0);
std::transform(loads.begin(), loads.end(), loads.begin(),
std::bind(std::multiplies<double>(), _1, factor));
- for (unsigned i = 0 ; i < hosts.size() ; ++i)
- set_load(i, loads[i]);
+ if (opt::integer_transfer) {
+ // round the loads
+ std::vector<double> iloads(hosts.size());
+ std::transform(loads.begin(), loads.end(), iloads.begin(), round);
+ // compute the differences between each load and its rounded value
+ std::vector<double> diffs(hosts.size());
+ std::transform(loads.begin(), loads.end(), iloads.begin(),
+ diffs.begin(), std::minus<double>());
+ // compute the absolute values of the diffs
+ std::vector<double> adiffs(hosts.size());
+ std::transform(diffs.begin(), diffs.end(), adiffs.begin(), fabs);
+ // find i, index of the element farthest from its rounded value
+ unsigned i;
+ i = std::max_element(adiffs.begin(), adiffs.end()) - adiffs.begin();
+ // remove element i from diffs, and compute the residual part...
+ diffs[i] = diffs.back();
+ diffs.pop_back();
+ double residue = std::accumulate(diffs.begin(), diffs.end(), 0.0);
+ // ... and compute element i (rounded to avoid numerical errors)
+ iloads[i] = fabs(round(loads[i] + residue));
+ // final sanity check
+ xbt_assert(opt::auto_depl::load ==
+ std::accumulate(iloads.begin(), iloads.end(), 0.0));
+ for (unsigned i = 0 ; i < hosts.size() ; ++i)
+ set_load(i, iloads[i]);
+ } else {
+ for (unsigned i = 0 ; i < hosts.size() ; ++i)
+ set_load(i, loads[i]);
+ }
}
void deployment_generator::deploy()