I’ve implemented a version of a model to use map_rect() and it appears as if it works just like the original version. But I’m not sure how to run it in rstan across all of the available processors. I’ll probably be running the model as an sbatch job on a MPI cluster, or locally on a 20-core machine. I’m using rstan, but I’ve looked through the manual, on CRAN, and on mc-stan.org, and I can’t find out how to start using all of the available resources. Also, if I run 4 chains on a 20 core machine, does rstan simply work correctly to distribute the resources as necessary?
After getting this working on my personal computer, I’ve been trying to get it on a high performance computer cluster. However, unlike on my personal computer, it doesn’t appear to be using any more cores than there are chains. So, right now, I have 4 chains on a 20-core system, and htop is showing me just 4 cores at 100% and the rest at 0%. This is on a CentOS machine. My .R/Makevars is: