I don't know if this is in PyMC3 scope but I was thinking since we have variational inference perhaps we could consider adding [INLA](https://inla.r-inla-download.org/r-inla.org/papers/inla-rss.pdf) [You can find a brief intro is here](https://www.precision-analytics.ca/blog-1/inla)