1 /* Copyright (c) 2007-2011, 2013-2017. The SimGrid Team.
2 * All rights reserved. */
4 /* This program is free software; you can redistribute it and/or modify it
5 * under the terms of the license (GNU LGPL) which comes with this package. */
7 #include "src/kernel/lmm/maxmin.hpp"
9 #include "xbt/sysdep.h"
15 XBT_LOG_EXTERNAL_DEFAULT_CATEGORY(surf_maxmin);
16 #define SHOW_EXPR_G(expr) XBT_DEBUG(#expr " = %g", expr);
17 #define SHOW_EXPR_D(expr) XBT_DEBUG(#expr " = %d", expr);
18 #define SHOW_EXPR_P(expr) XBT_DEBUG(#expr " = %p", expr);
20 void simgrid::kernel::lmm::bottleneck_solve(lmm_system_t sys)
23 lmm_element_t elem = nullptr;
24 xbt_swag_t elem_list = nullptr;
26 if (not sys->modified)
29 XBT_DEBUG("Variable set : %zu", sys->variable_set.size());
30 for (s_lmm_variable_t& var : sys->variable_set) {
32 XBT_DEBUG("Handling variable %p", &var);
33 if (var.sharing_weight > 0.0 && std::find_if(begin(var.cnsts), end(var.cnsts), [](s_lmm_element_t const& x) {
34 return x.consumption_weight != 0.0;
35 }) != end(var.cnsts)) {
36 sys->saturated_variable_set.push_back(var);
38 XBT_DEBUG("Err, finally, there is no need to take care of variable %p", &var);
39 if (var.sharing_weight > 0.0)
44 XBT_DEBUG("Active constraints : %zu", sys->active_constraint_set.size());
45 for (s_lmm_constraint_t& cnst : sys->active_constraint_set) {
46 sys->saturated_constraint_set.push_back(cnst);
48 for (s_lmm_constraint_t& cnst : sys->saturated_constraint_set) {
49 cnst.remaining = cnst.bound;
53 XBT_DEBUG("Fair bottleneck Initialized");
56 * Compute Usage and store the variables that reach the maximum.
58 auto& var_list = sys->saturated_variable_set;
59 auto& cnst_list = sys->saturated_constraint_set;
61 if (XBT_LOG_ISENABLED(surf_maxmin, xbt_log_priority_debug)) {
62 XBT_DEBUG("Fair bottleneck done");
65 XBT_DEBUG("******* Constraints to process: %zu *******", cnst_list.size());
66 for (auto iter = std::begin(cnst_list); iter != std::end(cnst_list);) {
67 s_lmm_constraint_t& cnst = *iter;
69 XBT_DEBUG("Processing cnst %p ", &cnst);
70 elem_list = &cnst.enabled_element_set;
72 xbt_swag_foreach(_elem, elem_list)
74 elem = static_cast<lmm_element_t>(_elem);
75 xbt_assert(elem->variable->sharing_weight > 0);
76 if (elem->consumption_weight > 0 && elem->variable->saturated_variable_set_hook.is_linked())
79 XBT_DEBUG("\tThere are %d variables", nb);
80 if (nb > 0 && not cnst.sharing_policy)
85 iter = cnst_list.erase(iter);
87 cnst.usage = cnst.remaining / nb;
88 XBT_DEBUG("\tConstraint Usage %p : %f with %d variables", &cnst, cnst.usage, nb);
93 for (auto iter = std::begin(var_list); iter != std::end(var_list);) {
94 s_lmm_variable_t& var = *iter;
95 double min_inc = DBL_MAX;
96 for (s_lmm_element_t const& elm : var.cnsts) {
97 if (elm.consumption_weight > 0)
98 min_inc = std::min(min_inc, elm.constraint->usage / elm.consumption_weight);
101 min_inc = std::min(min_inc, var.bound - var.value);
103 XBT_DEBUG("Updating variable %p maximum increment: %g", &var, var.mu);
105 if (var.value == var.bound)
106 iter = var_list.erase(iter);
111 for (auto iter = std::begin(cnst_list); iter != std::end(cnst_list);) {
112 s_lmm_constraint_t& cnst = *iter;
113 XBT_DEBUG("Updating cnst %p ", &cnst);
114 elem_list = &cnst.enabled_element_set;
115 xbt_swag_foreach(_elem, elem_list)
117 elem = static_cast<lmm_element_t>(_elem);
118 xbt_assert(elem->variable->sharing_weight > 0);
119 if (cnst.sharing_policy) {
120 XBT_DEBUG("\tUpdate constraint %p (%g) with variable %p by %g", &cnst, cnst.remaining, elem->variable,
122 double_update(&cnst.remaining, elem->consumption_weight * elem->variable->mu, sg_maxmin_precision);
124 XBT_DEBUG("\tNon-Shared variable. Update constraint usage of %p (%g) with variable %p by %g", &cnst,
125 cnst.usage, elem->variable, elem->variable->mu);
126 cnst.usage = std::min(cnst.usage, elem->consumption_weight * elem->variable->mu);
129 if (not cnst.sharing_policy) {
130 XBT_DEBUG("\tUpdate constraint %p (%g) by %g", &cnst, cnst.remaining, cnst.usage);
132 double_update(&cnst.remaining, cnst.usage, sg_maxmin_precision);
135 XBT_DEBUG("\tRemaining for %p : %g", &cnst, cnst.remaining);
136 if (cnst.remaining <= 0.0) {
137 XBT_DEBUG("\tGet rid of constraint %p", &cnst);
139 iter = cnst_list.erase(iter);
140 xbt_swag_foreach(_elem, elem_list)
142 elem = static_cast<lmm_element_t>(_elem);
143 if (elem->variable->sharing_weight <= 0)
145 if (elem->consumption_weight > 0 && elem->variable->saturated_variable_set_hook.is_linked()) {
146 XBT_DEBUG("\t\tGet rid of variable %p", elem->variable);
147 var_list.erase(var_list.iterator_to(*elem->variable));
154 } while (not var_list.empty());
157 sys->modified = true;
158 if (XBT_LOG_ISENABLED(surf_maxmin, xbt_log_priority_debug)) {
159 XBT_DEBUG("Fair bottleneck done");