I am looking for an optimisation routine within scipy/numpy which could solve a non-linear least-squares type problem (e.g., fitting a parametric function to a large dataset
Have a look at: http://lmfit.github.io/lmfit-py/, it should solve your problem.