I stumbled upon this problem as well. I want to perform local minimization with least_sq for my basinhopping global minimization, instead of using arbitrary methods such as BFGS, Powell, etc. This is because my residuals (objective function) is exactly a nonlinear least-square equation.
I was wondering if you were able to solve it?
Best, Matt