+
+ /* d Dual
+ * Compute the value of ------------- (\lambda^k, \mu^k) this portion
+ * d \lambda_i^k
+ * of code depends on function f(x).
+ */
+ cnst_list=&(sys->active_constraint_set);
+ xbt_swag_foreach(cnst1, cnst_list) {
+
+ lambda_partial = 0;
+
+ elem_list = &(cnst1->active_element_set);
+
+ xbt_swag_foreach(elem1, elem_list) {
+ lambda_partial = 0;
+
+ var2 = elem1->variable;
+
+ //for each link with capacity cnsts[i] that uses flow of variable var1 do
+ for(i=0; i<var2->cnsts_number; i++){
+ elem2 = &(var2->cnsts[i]);
+ tmp += (elem2->constraint)->bound + var2->bound;
+ }
+
+ lambda_partial += -1 / tmp;
+ }
+
+ lambda_partial += cnst1->initial_bound;
+ cnst1->bound = cnst1->bound + sigma_step * lambda_partial;
+ }
+