Nothing
### SHELL> mpiexec -np 2 Rscript --vanilla [...].r
library(pbdMPI, quiet = TRUE)
library(pbdNCDF4, quiet = TRUE)
### default of pbdMPI
rank <- comm.rank()
size <- comm.size()
### parallel open
if(! file.exists("test_par.nc")){
stop("test_par.nc does not exist. Run ncwrite_par to create.")
}
nc <- nc_open_par("test_par.nc")
nc_var_par_access(nc, "testMatrix")
### get data dimension.
n <- length(nc$dim$rows$vals)
p <- length(nc$dim$columns$vals)
### set up slab for this process
st <- c(rank*n/size + 1, 1)
co <- c(n/size , p)
### read variable values to file
x <- ncvar_get(nc, "testMatrix", start = st, count = co)
### Print results.
comm.cat("n = ", n, " p = ", p, "\n", sep = "")
comm.print(x, all.rank = TRUE)
### close file
nc_close(nc)
ncdump("test_par.nc")
finalize()
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.