Run simple R/MPI example on CSC Taito

rmpi_script.R

# Tell all slaves to return a message identifying themselves. 
mpi.remote.exec(paste("I am",mpi.comm.rank(),"of",mpi.comm.size())) 
mpi.remote.exec(paste(mpi.comm.get.parent()))

# Send execution commands to the slaves
x <- 5
x <- mpi.remote.exec(rnorm,x) 
length(x)
x

mpi.close.Rslaves()
mpi.quit()

run_r_parallel.sh

#!/bin/bash
#SBATCH -J rmpi_job
#SBATCH -o output.txt
#SBATCH -e errors.txt
#SBATCH -t hh:mm:ss
#SBATCH --nodes=1
#SBATCH --ntasks=12
#SBATCH --mem-per-cpu=200
#
module load R.latest
srun -n 12 Rmpi --no-save < rmpi_script.R