1 comment

Hi Dominique,

m.addConstr(y - radius <= -bd, "ellipsoid_UBy ")

m.addConstr(y - radius >= -bd, "ellipsoid_UBy ")
Moreover, the default lower bound for optimization variables is $$0$$. Once you make your coordinate variables $$x,y$$ free variables by setting
x = m.addVar(lb=-GRB.INFINITY,name="x")y = m.addVar(lb=-GRB.INFINITY,name="y")