library(rlecuyer) library(snow) machines <- c("hyland", "dogleg", "crab", "sugar", "strike", "pool", "wagonga", "tamarac") cl <- makeSOCKcluster(machines) clusterCall(cl, function() Sys.info()[c("nodename","machine")]) clusterSetupRNG(cl) ### see rcdd package vignette, pages 8 and 16 clusterEvalQ(cl, library(rcdd)) d <- 5 nsplit <- 200 clusterExport(cl, c("d", "nsplit")) out <- clusterEvalQ(cl, xsplit <- matrix(rnorm(d * nsplit), ncol = d)) out <- clusterEvalQ(cl, xsplit <- qdq(d2q(1e6 * xsplit), rep(d2q(1e6), length(xsplit)))) x <- NULL for (i in 1:length(out)) x <- rbind(x, out[[i]]) dim(x) out <- clusterEvalQ(cl, foo <- cbind(0, cbind(1, xsplit))) out <- clusterEvalQ(cl, redundant(foo, representation = "V")) sapply(out, function(foo) nrow(foo$output)) stopCluster(cl)