xref JuliaNLSolvers/Optim.jl#977, specially point 1 of the TODO list:
The passing of array type for the Jacobian of the constraints is ugly. I made it an argument to IPNewton, but it would be better as a field of TwiceDifferentiableConstraints, just like H is a field of TwiceDifferentiable. That requires changing NLSolversBase