If you really need a faster approximation than MCMC, try ADVI or pathfinder. We just buffed the implementation of pathfinder in pymc-extras, so that’s worth a second look if you’ve tried it before. There’s also a MAP+Laplace approximation in extras that can be good too (if MAP works on your problem in the first place, of course). Example of that one is here.