Pm.sample gets stuck after init with cores > 1

Have you tried specifiying your priors a bit more? Maybe do pm.Normal(0,5) for predictors/ intercepts?

Yup. I’ve tried more informative priors like pm.Normal.dist(0., 1.) and couldn’t get the sampler to start.