If the constraint defines a differentiable manifold, then I would use Manopt.jl . Obviously it does not allow inequality constraints. I find that understanding the manifold structure of the constraint has a huge advantage for efficiency.
1 Like