setup script to cache artificial interim posterior ellipticities

[ismael-mendoza]

We can start with NUTS with low MD and later iterate and optimize. The advatange is that we don't need to tune to get high efficiency and leaves us the option of using GPUs later if want to