try methods suited to optimisation of a noisy function

[grahamgower]

https://noisyopt.readthedocs.io/en/latest/
https://github.com/scikit-optimize/scikit-optimize

[grahamgower]

https://clima.github.io/EnsembleKalmanProcesses.jl/dev/ensemble_kalman_inversion/

[grahamgower]

mackelab.org/sbi/ looks amazing!