MAINT: for model fitting, let lbfgs use its own fprime approximation instead of passing a separate function that uses finite differences