Pm.sample gets stuck after init with cores > 1

I would say that you should be able to sample from a model like that, though it will take time and some restructuring to handle the parallelization properly. Maybe @twiecki can give some advice on the subject. I’ve dealt with parallelism by ditching pythons multiprocessing module and running several sample jobs concurrently with different initialization seeds.