#include <math.h>
XBT_LOG_EXTERNAL_DEFAULT_CATEGORY(surf_maxmin);
+#define SHOW_EXPR_G(expr) DEBUG1(#expr " = %g",expr);
+#define SHOW_EXPR_D(expr) DEBUG1(#expr " = %d",expr);
+#define SHOW_EXPR_P(expr) DEBUG1(#expr " = %p",expr);
void bottleneck_solve(lmm_system_t sys)
{
cnst->usage = 0.0;
}
- if (XBT_LOG_ISENABLED(surf_maxmin, xbt_log_priority_debug)) {
- DEBUG0("Fair bottleneck Init");
- lmm_print(sys);
- }
+ DEBUG0("Fair bottleneck Initialized");
/*
* Compute Usage and store the variables that reach the maximum.
*/
while (1) {
+/* if (XBT_LOG_ISENABLED(surf_maxmin, xbt_log_priority_debug)) { */
+/* DEBUG0("Fair bottleneck done"); */
+/* lmm_print(sys); */
+/* } */
DEBUG1("******* Constraints to process: %d *******", xbt_swag_size(cnst_list));
+ min_usage = -1;
xbt_swag_foreach_safe(cnst, cnst_next, cnst_list) {
int nb = 0;
double max_elem = 0.0;
xbt_swag_remove(cnst, cnst_list);
continue;
}
+
/* Saturated constraints update */
if (min_usage < 0 || min_usage > cnst->remaining / cnst->usage) {
DEBUG3("Update min_usage (%g) with cnst %p -> %g",min_usage, cnst,
DEBUG2("Keep min_usage (%g) with cnst %p",min_usage, cnst);
xbt_swag_remove(cnst, cnst_list);
xbt_swag_insert(cnst, &(cnst_to_update));
+ } else {
+ DEBUG1("\tmin_usage: %f. No update",min_usage);
}
}
while ((cnst_next = xbt_swag_extract(&cnst_to_update))) {
int nb = 0;
+ double remaining = cnst_next->remaining;
elem_list = &(cnst_next->element_set);
xbt_swag_foreach(elem, elem_list) {
if (elem->variable->weight <= 0)
DEBUG1("Updating for cnst %p",cnst_next);
xbt_swag_foreach(elem, elem_list) {
+ if (elem->value <= 0) continue;
+
var = elem->variable;
+
if (var->weight <= 0)
break;
if (var->value == 0.0) {
- var->value = var->weight * cnst_next->remaining / (nb * elem->value);
+ DEBUG2("\tUpdating var %p (%g)",var,var->value);
+ var->value = var->weight * remaining / (nb * elem->value);
/* Update usage */
- DEBUG1("\tUpdating var %p",var);
for (i = 0; i < var->cnsts_number; i++) {
lmm_element_t elm = &var->cnsts[i];
cnst = elm->constraint;
sys->modified = 0;
if (XBT_LOG_ISENABLED(surf_maxmin, xbt_log_priority_debug)) {
+ DEBUG0("Fair bottleneck done");
lmm_print(sys);
}
}