Bug in dual methods solved

bslv_algs.c

@@ -1267,7 +1267,6 @@ void phase1_dual(soltype *const sol, const vlptype *vlp, const opttype *opt, lpt
 	
 	init_P1(sol, vlp, HOMOGENEOUS, lpstr);	// assume that appropriate eta is stored in sol->eta
 	lp_set_options(opt, PHASE1_DUAL, lpstr);
-
 	/* PHASE 1 -- DUAL -- PART 1
 	 *
 	 * determine first outer approximation of upper image
@@ -1334,7 +1333,7 @@ void phase1_dual(soltype *const sol, const vlptype *vlp, const opttype *opt, lpt
 			assert(lp_status == LP_OPTIMAL);
 			double opt_val=lower_image.val[lower_image.dim-1]; // store value before lower_image.val is used for input
 			lp_primal_solution_cols(lower_image.val,vlp->n+1,vlp->q, 1, lpstr);
-			if (opt_val-lp_obj_val(0)>opt->eps_benson_phase1)
+			if (opt_val-lp_obj_val(lpstr)>opt->eps_benson_phase1)
 			{
 				if (opt->message_level >= 3)
 					printf("add primal vertex\n");
@@ -1495,7 +1494,7 @@ void phase2_dual(soltype *const sol, const vlptype *vlp, const opttype *opt, lpt
 			}
 			double opt_val=lower_image->val[lower_image->dim-1]; // store this value before lower_image->val is used for input
 			lp_primal_solution_cols(lower_image->val,vlp->n+1,vlp->q, 1, lpstr);
-			if (opt_val-lp_obj_val(0)>opt->eps_benson_phase2)
+			if (opt_val-lp_obj_val(lpstr)>opt->eps_benson_phase2)
 			{
 				if (opt->message_level >= 3)
 					printf("add primal vertex\n");