Alan_Weiss <firstname.lastname@example.org> wrote in message <email@example.com>... > > When using the sqp or interior-point algorithms there is no need to set > x(i)>=eps; these algorithms obey bounds at all iterations, and your > extension of the definition to the boundary is fine. =================
I'm not so sure. The function is not differentiable at the the boundary. The directional derivative goes to infinity when approached from the right and has no definition when approached from the left. This violates a lot of textbook assumptions.
A solution might be to make the change of variables
leading to the reformulated problem
J = summmation(i = 1:N) k_i*exp(z_i)*(z_i*+ c_i)
Subject to : summation(i=1:N) k_i*exp(z_i) = 1
This is now finitely differentiable everywhere. It's true that this breaks convexity, but it is a monotonic change of variables, so the problem should still be unimodal (I think).