Hi authors!

I am exploring your package and like the simplicity of the implementation to adaptively sample search spaces. Are there any ideas on adding uncertainty quantification to the sampling strategies like you would have when using Gaussian processes for Bayesian optimization with a space-filling / expected improvement acquisition function? Any response on this would be much appreciated. Thank you very much for maintaining this project!