Switched from using pop to get in _initiateSolverAndObjFunc.  It was causing the fit to fail on large datasets.