MATLAB: Handling bound constraints by the Levenberg-Marquardt algorithm

bound constraintslevenberg-marquardt algorithmlsqnonlinOptimization Toolbox

Hello,
Could anyone please tell me why the Levenberg-Marquardt algorithm embedded in "lsqnonlin" does not handle bound constraints, while the other one ("trust-region-reflective") does?
There are implementations of the Levenberg-Marquardt algorithm that do accept bound constraints, so what is the principle limitation explaining why this has not been implemented in "lsqnonlin"?
Thank you!
Igor.

Best Answer

The theory of Levenberg-Marquardt does not define a way to handle bound constraints. If, as you claim, there are modifications of classical LM that support bounds, I surmise that they involve manipulations similar in spirit to what is done in the trust-region algorithm.