Sign glitch in `detect_simple_bounds`?

[gregreich]

Sign glitch in `detect_simple_bounds`?

[ernestds]

Even simpler example:

import casadi as ca
x = ca.MX.sym("x",1)
p = ca.MX.sym("p",1)
g = -x
print(ca.detect_simple_bounds(x,p,g,ca.MX([-2]),ca.MX([-2])))

g = x
print(ca.detect_simple_bounds(x,p,g,ca.MX([-2]),ca.MX([-2])))

Half assed hacky fix?

x = ca.MX.sym("x",2)
p = ca.MX.sym("p",2)
lbg = ca.MX([-2, -2])
ubg = ca.MX([2, -1])

g = p[0] * x + p[1]
[gi,lbx,ubx,lam_f,lam_b] = ca.detect_simple_bounds(x,p,g,lbg,ubg)
dx = 1.
simple_bounds_idx = list(set(range(g.numel())) - set(gi))
f_g = ca.Function('x',[x,p], [g[simple_bounds_idx]],)
f_lbg = ca.Function('x',[x,p], [lbg[simple_bounds_idx]],)
f_ubg = ca.Function('x',[x,p], [ubg[simple_bounds_idx]],)
constant = f_g(0,p)
coeff = (f_g(dx,p) - constant)/dx

lbg = f_lbg(0,p)
ubg = f_ubg(0,p)

lbx_f = ca.Function('lbx',[p],[ca.if_else(coeff >= 0, (lbg - constant)/coeff,(ubg - constant)/coeff)])
ubx_f = ca.Function('ubx',[p],[ca.if_else(coeff >= 0, (ubg - constant)/coeff,(lbg - constant)/coeff)])
p0 = [-1,2]
print(lbx_f(p0))
print(ubx_f(p0))

p0 = [1,2]
print(lbx_f(p0))
print(ubx_f(p0))