demo/cslavePI.R

In Rmpi: Interface (Wrapper) to MPI (Message-Passing Interface)

cslavePI <- function (n)
{
	if (mpi.comm.size(1) > 1)
		stop ("It seems some slaves running on comm 1.")
        mpi.comm.spawn("cslavePI")
	mpi.intercomm.merge(2,0,1)
        mpi.bcast(as.integer(n),1)
        out <-mpi.reduce(0)
        mpi.comm.free()
        out
}