+ elem->variable = var;
+
+ if(var->weight) xbt_swag_insert_at_head(elem,&(elem->constraint->element_set));
+ else xbt_swag_insert_at_tail(elem,&(elem->constraint->element_set));
+
+ make_constraint_active(sys, cnst);
+}
+
+void lmm_expand_add(lmm_system_t sys, lmm_constraint_t cnst,
+ lmm_variable_t var, double value)
+{
+ int i ;
+ sys->modified = 1;
+
+ for(i=0; i< var->cnsts_number ; i++)
+ if(var->cnsts[i].constraint == cnst) break;
+
+ if(i<var->cnsts_number) var->cnsts[i].value +=value;
+ else lmm_expand(sys,cnst,var,value);
+}
+
+void lmm_elem_set_value(lmm_system_t sys, lmm_constraint_t cnst,
+ lmm_variable_t var, double value)
+{
+ int i ;
+
+ for(i=0; i< var->cnsts_number ; i++)
+ if(var->cnsts[i].constraint == cnst) break;
+
+ if(i<var->cnsts_number) {
+ var->cnsts[i].value =value;
+ sys->modified = 1;
+ } else DIE_IMPOSSIBLE;
+}
+
+lmm_constraint_t lmm_get_cnst_from_var(lmm_system_t sys,
+ lmm_variable_t var, int num)
+{
+ if (num < var->cnsts_number)
+ return (var->cnsts[num].constraint);
+ else
+ return NULL;
+}
+
+int lmm_get_number_of_cnst_from_var(lmm_system_t sys, lmm_variable_t var)
+{
+ return (var->cnsts_number);
+}
+
+lmm_variable_t lmm_get_var_from_cnst(lmm_system_t sys,
+ lmm_constraint_t cnst,
+ lmm_variable_t * var)
+{
+ if (!(*var))
+ xbt_swag_getFirst(&(cnst->element_set));
+ else
+ *var = xbt_swag_getNext(*var, cnst->element_set.offset);
+ return *var;
+}
+
+void *lmm_constraint_id(lmm_constraint_t cnst)
+{
+ return cnst->id;
+}
+
+void *lmm_variable_id(lmm_variable_t var)
+{
+ return var->id;
+}
+
+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);
+ }
+}
+
+static 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';
+
+ /* 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->weight);
+ 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);
+ DEBUG1("%s",trace_buf);
+ trace_buf[0]='\000';
+ xbt_assert2((sum<=cnst->bound), "Incorrect value (%f is not smaller than %f)",
+ sum,cnst->bound);
+ }