I am comparing a result of optimization using R and Mathematica.
NMaximize) has found a significantly better (global optimum) result then R (using L-BFGS-B).
The function has 18 variables. The R routine stopped at a place, where none of the gradient is zero. So it's not even a local optimum. I am trying to find out why.
I wonder if I can duplicate the same result in Mathematica using the same L-BFGS-B method?