I’ve recently been experimenting more with pymc3’s library for probabilistic programming, specifically Bayesian networks. It was insanely slow and I was getting around 20/samples a second from a very benign y=b1X1 + b2X1 + alpha
model.
Thinking this was just my slow machine, I kept going, until I …