Experiment: SMC sampler
https://docs.pymc.io/notebooks/SMC2_gaussians.html
Sadly, with Cyclus this algorithm didn't make any progress. I haven't bothered looking further into optimizing it yet, because Slice already gives us good results even for multimodal distributions (with random start values), but SMC is generally harder to adjust as it is less transparent (auto-tuning etc.)