+static void saturated_constraint_set_update(lmm_system_t sys,
+ lmm_constraint_t cnst,
+ double *min_usage)
+{
+ lmm_constraint_t useless_cnst = NULL;
+
+ XBT_IN3("sys=%p, cnst=%p, min_usage=%f", sys, cnst, *min_usage);
+ if (cnst->usage <= 0) {
+ XBT_OUT;
+ return;
+ }
+ if (cnst->remaining <= 0) {
+ XBT_OUT;
+ return;
+ }
+ if ((*min_usage < 0) || (*min_usage > cnst->remaining / cnst->usage)) {
+ *min_usage = cnst->remaining / cnst->usage;
+ LOG3(xbt_log_priority_trace,
+ "min_usage=%f (cnst->remaining=%f, cnst->usage=%f)", *min_usage,
+ cnst->remaining, cnst->usage);
+ while ((useless_cnst =
+ xbt_swag_getFirst(&(sys->saturated_constraint_set))))
+ xbt_swag_remove(useless_cnst, &(sys->saturated_constraint_set));
+
+ xbt_swag_insert(cnst, &(sys->saturated_constraint_set));
+ } else if (*min_usage == cnst->remaining / cnst->usage) {
+ xbt_swag_insert(cnst, &(sys->saturated_constraint_set));
+ }
+ XBT_OUT;
+}
+
+static void saturated_variable_set_update(lmm_system_t sys)
+{
+ lmm_constraint_t cnst = NULL;
+ xbt_swag_t cnst_list = NULL;
+ lmm_element_t elem = NULL;
+ xbt_swag_t elem_list = NULL;
+
+ cnst_list = &(sys->saturated_constraint_set);
+ while ((cnst = xbt_swag_getFirst(cnst_list))) {
+ elem_list = &(cnst->active_element_set);
+ xbt_swag_foreach(elem, elem_list) {
+ if (elem->variable->weight <= 0)
+ break;
+ if ((elem->value > 0))
+ xbt_swag_insert(elem->variable, &(sys->saturated_variable_set));
+ }
+ xbt_swag_remove(cnst, cnst_list);
+ }
+}
+
+void lmm_print(lmm_system_t sys)
+{
+ lmm_constraint_t cnst = NULL;
+ lmm_element_t elem = NULL;
+ lmm_variable_t var = NULL;
+ xbt_swag_t cnst_list = NULL;
+ xbt_swag_t var_list = NULL;
+ xbt_swag_t elem_list = NULL;
+ char print_buf[1024];
+ char *trace_buf = xbt_malloc0(sizeof(char));
+ double sum = 0.0;
+
+ /* Printing Objective */
+ var_list = &(sys->variable_set);
+ sprintf(print_buf, "MAX-MIN ( ");
+ trace_buf =
+ xbt_realloc(trace_buf, strlen(trace_buf) + strlen(print_buf) + 1);
+ strcat(trace_buf, print_buf);
+ xbt_swag_foreach(var, var_list) {
+ sprintf(print_buf, "'%p'(%f) ", var, var->weight);
+ trace_buf =
+ xbt_realloc(trace_buf, strlen(trace_buf) + strlen(print_buf) + 1);
+ strcat(trace_buf, print_buf);
+ }
+ sprintf(print_buf, ")");
+ trace_buf =
+ xbt_realloc(trace_buf, strlen(trace_buf) + strlen(print_buf) + 1);
+ strcat(trace_buf, print_buf);
+ DEBUG1("%s", trace_buf);
+ trace_buf[0] = '\000';
+
+ DEBUG0("Constraints");
+ /* Printing Constraints */
+ cnst_list = &(sys->active_constraint_set);
+ xbt_swag_foreach(cnst, cnst_list) {
+ sum = 0.0;
+ elem_list = &(cnst->element_set);
+ sprintf(print_buf, "\t");
+ trace_buf =
+ xbt_realloc(trace_buf, strlen(trace_buf) + strlen(print_buf) + 1);
+ strcat(trace_buf, print_buf);
+ xbt_swag_foreach(elem, elem_list) {
+ sprintf(print_buf, "%f.'%p'(%f) + ", elem->value,
+ elem->variable, elem->variable->value);
+ trace_buf =
+ xbt_realloc(trace_buf,
+ strlen(trace_buf) + strlen(print_buf) + 1);
+ strcat(trace_buf, print_buf);
+ sum += elem->value * elem->variable->value;
+ }
+ sprintf(print_buf, "0 <= %f ('%p')", cnst->bound, cnst);
+ trace_buf =
+ xbt_realloc(trace_buf, strlen(trace_buf) + strlen(print_buf) + 1);
+ strcat(trace_buf, print_buf);
+
+ if (!cnst->shared) {
+ sprintf(print_buf, " [MAX-Constraint]");
+ trace_buf =
+ xbt_realloc(trace_buf,
+ strlen(trace_buf) + strlen(print_buf) + 1);
+ strcat(trace_buf, print_buf);
+ }
+ DEBUG1("%s", trace_buf);
+ trace_buf[0] = '\000';
+ xbt_assert3(!double_positive(sum - cnst->bound),
+ "Incorrect value (%f is not smaller than %f): %g",
+ sum, cnst->bound, sum - cnst->bound);
+ }
+
+ DEBUG0("Variables");
+ /* Printing Result */
+ xbt_swag_foreach(var, var_list) {
+ if (var->bound > 0) {
+ DEBUG4("'%p'(%f) : %f (<=%f)", var, var->weight, var->value,
+ var->bound);
+ xbt_assert2(!double_positive(var->value - var->bound),
+ "Incorrect value (%f is not smaller than %f",
+ var->value, var->bound);
+ } else
+ DEBUG3("'%p'(%f) : %f", var, var->weight, var->value);
+ }
+
+ free(trace_buf);