I think, eventually (?) we need to use a more general setup that allows for the choice of optimizer and options.
lbfgsb might be also a good choice.
It would currently be possible by calling super().fit but needs a loglike method to optimize, Skipper started to outsource our optimizer wrapper.