Python crashes when trying to sample model with covariates

Indeed, that avoids crashing pm.sample, in my setup as well. But is there a better solution, that would allow one to use the multi-cores? Has the issue evolved since 2018?
I am using pymc 3.11.0.
Thanks!