Implement the interval constraint propagation scheme presented in Vu 2008. For improved convergences.
A parametric bisection routine will be updated that can divide the (X,P) space into a series of boxes that all contain unique branches of the implicit function p->y(p).
Provide a better interface the nonconvex semi-infinite programs solvers (JuMPeR extension?).
Add additional McCormick relaxations.
Add handling for domain reduction of special expression forms.