great, thank you, I included a tolerance on the x at convergence and now it returns a solution, so just to be clear you did the following:
1- transposed the constraint gradient
2- you used ```inequality_constraint!(opt, (z,x,g) → ps_con(z,x,g,[1,1]), [1e-8,1e-8])
question: if the solver does not need gradients does it matter how I define the gradient? I guess no