Support global bayesian optimization

[HDembinski]

Finding starting values for a fit is often challenging. Bayesian optimization is a global optimization strategy which could be very useful in that regard. Bayesian optimization is implemented, for example, in hyperopt.