- if((var->bound >= 0) && (var->weight > 0) ){
- var->new_mu = dichotomy(var->mu, partial_diff_mu, var, dichotomy_min_error);
- if(var->new_mu < 0) var->new_mu = 0;
- DEBUG3("====> var->mu (%p) : %g -> %g", var, var->mu, var->new_mu);
+ if ((var->bound >= 0) && (var->weight > 0)) {
+ DEBUG1("Working on var (%p)", var);
+ var->new_mu =
+ dichotomy(var->mu, partial_diff_mu, var, dichotomy_min_error);
+ dual_updated += (fabs(var->new_mu-var->mu)>dichotomy_min_error);
+ DEBUG2("dual_updated (%d) : %1.20f",dual_updated,fabs(var->new_mu-var->mu));
+ DEBUG3("Updating mu : var->mu (%p) : %1.20f -> %1.20f", var, var->mu, var->new_mu);