Displaying 2 results from an estimated 2 matches for "projseeds".
2012 Jan 13
0
Example of "task seeds" with R parallel. Critique?
...#39;s why I wish I could integrate task-based
seeds into parallel itself.
RNGkind("L'Ecuyer-CMRG")
set.seed(23456)
library(parallel)
## nrep = number of repetitions (or tasks)
## streamsPerRep = number of streams needed by each repetition
nReps <- 2000
streamsPerRep <- 2
## projSeeds=list of lists of stream seeds
projSeeds <- vector(mode="list", nReps)
for (i in 1:nReps) projSeeds[[i]] <- vector(mode="list", streamsPerRep)
runif(1) ##establishes .Random.seed
##Grab first seed
s <- .Random.seed
origSeed <- s
x <- rnorm(4) ##will compare later...
2012 Feb 17
3
portable parallel seeds project: request for critiques
......
run1000 seed1000.1 seed1000.2 seed1000.3
This example provides 3 separate streams of random numbers within each
run. Because we will use the L'Ecuyer "many separate streams"
approach, we are confident that there is no correlation or overlap
between any of the runs.
The projSeeds has to have one row per project, but it is not a huge
file. I created seeds for 2000 runs of a project that requires 2 seeds
per run. The saved size of the file 104443kb, which is very small. By
comparison, a 1400x1050 jpg image would usually be twice that size.
If you save 10,000 runs-worth of se...