F(x) :
x0 :
Barier function Σ c(x, y) > -k ? -k/c(x, y) : 0, where k :
Constraints c :
0