search for: cmrg

...brary/base/man/Random.Rd says that L'Ecuyer requires a seed of length 6. ./src/library/parallel/man/RngStream.Rd also mentions this, but only in the text part; In the "Arguments"-part, it says that "seed" has to be of length 7 Also: ,---- | > RNGkind("L'Ecuyer-CMRG") | > length(.Random.seed) | [1] 7 `---- Is the docu old? Some clarification on the two .Rd files would be great. Cheers, Marius

...## create cluster object cl <- snow::makeCluster(parallel::detectCores(), type=type) on.exit(snow::stopCluster(cl)) ## shut down cluster and terminate execution environment registerDoSNOW(cl) ## register the cluster object with foreach ## seed if(seed=="L'Ecuyer-CMRG") { if(!exists(".Random.seed")) stop(".Random.seed does not exist - in l'Ecuyer setting") .t <- snow::clusterSetupRNG(cl, seed=.Random.seed[2:7]) # => fails! } ## actual work foreach(i=seq_len(n)) %dopar% { runif(1) } } ##...

Hi, This is most likely only a minor technicality, but I saw the following: On page 6 of the 'parallel' vignette (http://stat.ethz.ch/R-manual/R-devel/library/parallel/doc/parallel.pdf), the random-number generator "L'Ecuyer-CMRG" is said to have seed "(x_n, x_{n-1}, x_{n-2}, y_n, y_{n-1}, y_{n-2})". However, in L'Ecuyer et al. (2002), the seed is given with 'increasing' indices, so should rather be "(x_{n-2}, x_{n-1}, x_n, y_{n-2}, y_{n-1}, y_n)" (or, even more intuitively, "(x_{n-...

...am not yet aware of, so I'm asking you what can go wrong? I see danger in conflicts between my effort to manage seeds and the work of parallel functions that try to manage seeds for me. That'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:nR...

Dear All, Is the following expected behaviour? set.seed(1) library(parallel) cl = makeCluster(5) clusterSetRNGStream(cl, iseed = NULL) parSapply(cl, 1:5, function(i) sample(1:10, 1)) # 7 4 2 10 10 clusterSetRNGStream(cl, iseed = NULL) # 7 4 2 10 10 parSapply(cl, 1:5, function(i) sample(1:10, 1)) stopCluster(cl) The documentation could be read either way, e.g. * iseed: An integer to be

...at uwaterloo.ca> wrote: > Hi, > > This is most likely only a minor technicality, but I saw the > following: On page 6 of the 'parallel' vignette > (http://stat.ethz.ch/R-manual/R-devel/library/parallel/doc/parallel.pdf), > the random-number generator "L'Ecuyer-CMRG" is said to have seed > "(x_n, x_{n-1}, x_{n-2}, y_n, y_{n-1}, y_{n-2})". However, in L'Ecuyer > et al. (2002), the seed is given with 'increasing' indices, so should > rather be "(x_{n-2}, x_{n-1}, x_n, y_{n-2}, y_{n-1}, y_n)" (or, even > more intui...

In a vanilla R 2.14.0 GUI session (on Windows XP SP3): > library(parallel) > cl<-makePSOCKcluster(2) > RNGkind() [1] "Mersenne-Twister" "Inversion" > clusterSetRNGStream(cl) > RNGkind() [1] "L'Ecuyer-CMRG" "Inversion" > stopCluster(cl) Is it intentional that clusterSetRNGStream() changes the RNG kind in the master process? The code of clusterSetRNGStream() suggests that it is not: the old random number seed is saved in 'oldseed' and then assigned to '.Random.seed'...

...pposed to run on the RMPI cluster or just in a single workstation. The code for each run of the model begins by loading the required libraries and loading the seed file, if it exists, or generating a new "projSeed" object if it is not found. library(parallel) RNGkind("L'Ecuyer-CMRG") set.seed(234234) if (file.exists("projSeeds.rda")) { load("projSeeds.rda") } else { source("seedCreator.R") } ## Suppose the "run" number is: run <- 232 initSeedStreams(run) After that, R's random generator functions will draw values from...

...ode ID ddde6c2c45de8e60518aafa077f0f3867ff68e17 # Parent ccbee5bcb31b72706497725381f4e6836b9df657 libxl/xl: implement support for guest iooprt and irq permissions. This is useful for passing legacy ISA devices (e.g. com ports, parallel ports) to guests. Supported syntax is as described in http://cmrg.fifthhorseman.net/wiki/xen#grantingaccesstoserialhardwaretoadomU I tested this using Xen''s ''q'' key handler which prints out the I/O port and IRQ ranges allowed for each domain. e.g.: (XEN) Rangesets belonging to domain 31: (XEN) I/O Ports { 2e8-2ef, 2f8-2ff } (XEN)...

...ction in package tools now has a stylesheet argument, allowing pages to be displayed in alternate formats. o New function requireNamespace() analogous to require(), returning a logical value after attempting to load a namespace. o There is a new type of RNG, "L'Ecuyer-CMRG", implementing L'Ecuyer (1999)'s 'combined multiple-recursive generator' MRG32k3a. See the comments on ?RNG. o help.search() and ?? can now display vignettes and demos as well as help pages. The new option "help.search.types" controls the...

...ction in package tools now has a stylesheet argument, allowing pages to be displayed in alternate formats. o New function requireNamespace() analogous to require(), returning a logical value after attempting to load a namespace. o There is a new type of RNG, "L'Ecuyer-CMRG", implementing L'Ecuyer (1999)'s 'combined multiple-recursive generator' MRG32k3a. See the comments on ?RNG. o help.search() and ?? can now display vignettes and demos as well as help pages. The new option "help.search.types" controls the...