OptimizationMOI Ipopt violating inequality constraint

That’s a fair question. If you look at the assumptions underlying IPOPT, they say that all problem functions should be defined and differentiable over the entire space (\mathbb{R}^n). IPOPT belongs to the category of infeasible interior-point methods, first developed for linear optimization.

Feasible interior-point methods also exist in the literature, but are not typically implemented because they first require finding a strictly feasible initial guess (g(x) < 0), which can be as difficult as solving the problem in the first place. But of course, as is the case in your problem, if the objective is only defined when g(x) < 0, feasible methods are the only option.

All the details for implementing an efficient feasible interior-point method are in the literature. Alas, that’s a fairly nontrivial and time-consuming task.

You could get a trial version of KNITRO for free, at least to see if you manage to solve your problem to your satisfaction with it.

EDIT: since your problem is so small, you may be satisfied with a “prototype” implementation of a feasible interior-point method, one that doesn’t care about sparse linear algebra. It wouldn’t be too hard to get a preliminary version going.

I also notice in your original post that you have a strict inequality constraint. Optimization mostly considers closed feasible sets, so you’re unlikely to find a method that guarantees that the inequality will be satisfied strictly. However, feasible interior-point methods are probably your best bet.