Cross-chain warmup adaptation using MPI

Excellent

Good question. I’ll prepare them but before that @stevebronder, @rok_cesnovar or @jonah, how do I run

mpiexec -n 4 --tag-output ./radon sample data file=radon.data.R

Because if I now do

modelc1 = cmdstan_model("bernoulli_logit_glm_rhs.stan", quiet = FALSE)
fit1 = modelc1$sample(data = data)`

it compiles it correctly but doesn’t run it with mpi.