Active set: update numerical zero.

Before that patch, the active set numerical zero
(aka the threshold under which we disregard the value and say it's 0)
was way too high. This patch makes this threshold a bit more sensible.

src/LB.cpp

@@ -2187,22 +2187,21 @@ relax_info ActiveSetL2pL0::mixed_BigM_featuresign_l2pl0_withoutNode(mat A_S1, ve
 				c1 = A_SBar->t()*r;
 				c2 = A_S1.t()*r;
 			}
-
 			e = abs(c1(ind_x_0)) - lambda ;
-			i1 = arg_max_p(ind_x_0, e, &sigma1, &ie1, eps);
+			i1 = arg_max_p(ind_x_0, e, &sigma1, &ie1, eps_);
 
 			e = (-c1(ind_x_M) % sign(x_SBar(ind_x_M)))+ lambda ;
-			i2 = arg_max_p(ind_x_M, e, &sigma2, &ie2, eps);
+			i2 = arg_max_p(ind_x_M, e, &sigma2, &ie2, eps_);
 
 			e = -c2(ind_w_M) % sign(x_S1->rows(ind_w_M));
-			i3 = arg_max_p(ind_w_M, e, &sigma3, &ie3, eps);
+			i3 = arg_max_p(ind_w_M, e, &sigma3, &ie3, eps_);
 
 			if ( (sigma1<0) && (sigma2<0) && (sigma3<0) ){
 
 				grad = -c1(ind_x_in) % sign(x_SBar(ind_x_in)) + lambda;
 
 				if (!grad.is_empty()) {
-					if (max(abs( grad )) < eps) {
+					if (max(abs( grad )) < eps_) {
 						this->ctx->nbr_iter += it;
 
 						//return x_SBar;