+ DEBUG1("====> returning %e", (min+max)/2.0);
+ return ((min+max)/2.0);
+}
+
+/*
+ *
+ */
+double partial_diff_mu(double mu, void *param_var){
+ double mu_partial=0.0;
+ double sigma_mu=0.0;
+ lmm_variable_t var = (lmm_variable_t)param_var;
+ int i;
+
+ //compute sigma_i
+ for(i=0; i<var->cnsts_number; i++)
+ sigma_mu += (var->cnsts[i].constraint)->lambda;
+
+ //compute sigma_i + mu_i
+ sigma_mu += var->mu;
+
+ //use auxiliar function passing (sigma_i + mu_i)
+ mu_partial = diff_aux(var, sigma_mu) ;
+
+ //add the RTT limit
+ mu_partial += var->bound;
+
+ return mu_partial;
+}
+
+/*
+ *
+ */
+double partial_diff_lambda(double lambda, void *param_cnst){
+
+ int i;
+ xbt_swag_t elem_list = NULL;
+ lmm_element_t elem = NULL;
+ lmm_variable_t var = NULL;
+ lmm_constraint_t cnst= (lmm_constraint_t) param_cnst;
+ double lambda_partial=0.0;
+ double sigma_mu=0.0;
+
+ elem_list = &(cnst->element_set);
+
+ DEBUG2("Computting diff of cnst (%p) %s", cnst, (char *)cnst->id);
+
+ xbt_swag_foreach(elem, elem_list) {
+ var = elem->variable;
+ if(var->weight<=0) continue;