I wouldn't change your initial estimation, if that works much better, but it's still possible to call super().fit() inside the function where currently fmin_bfgs is.