Performance speedup for updating posterior with new data

Thank you, it should work.

Another quick question (maybe off-topic), is pm.Normal initialised this way (with shape>1) equal to Multivariate normal with only main diagonal in pymc3?