Thank you for your reply, Steven.
As far as I know, NLopt implements Newton with trust region for unconstrained or bound-constrained problems, but not for general non-linear equality and inequality constraints. Am I right?
I am looking for packages implementing trust-region constrained optimization method.