Could someone give some general instructions on how one can parallelize the PyMC MCMC code. I am trying to run LASSO regression following the example given here. I read somewhere that parallel sampling is done by default, but do I still need to use something like Parallel Python to get it to work? Here is some reference code that I would like to be able to parallelize on my machine. x1 = norm.rvs(