<div dir="ltr">Hi, 

 I have the following code that uses Pyomo and Couenne to optimize the energy consumption of a pumping system. The code works; however, something curious happens. I know the optimal result in advance, which is 3810. Most of the time, I get that result, but in some cases, the constraint is violated:

<div><br></div><div>  # Restricciones de flujo total por sistema<br>            def flujo_restriccion_rule(model):<br>                yield (model.flujo_bb_1_fase_total + model.flujo_bb_2_fase_total) - requerimientos_hidraulicos['BPC S1'][0] == 0<br>                yield (model.flujo_bb_1_fase_total + model.flujo_bb_2_fase_total) - model.flujo_bpc_s1_total == 0<br>                yield model.flujo_bpc_s1_total - requerimientos_hidraulicos['BPC S1'][0] == 0<br><br>            model.del_component('flujo_restricciones')<br>            model.flujo_restricciones = ConstraintList()<br><br>            for constraint in flujo_restriccion_rule(model):<br>                model.flujo_restricciones.add(constraint)<br></div></div>