1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
|
iRNGStream <- function(seed) {
# Convert a single number into the appropriate vector for "L'Ecuyer-CMRG"
if (length(seed) == 1) {
seed <- convseed(seed)
}
# Error checking: this will throw an error right away if the seed is bad
nextRNGStream(seed)
# Define the "Next Element" function for the iterator
nextEl <- function() (seed <<- nextRNGStream(seed))
obj <- list(nextElem=nextEl)
class(obj) <- c('abstractiter', 'iter')
obj
}
iRNGSubStream <- function(seed) {
# Convert a single number into the appropriate vector for "L'Ecuyer-CMRG"
if (length(seed) == 1) {
seed <- convseed(seed)
}
# Error checking: this will throw an error right away if the seed is bad
nextRNGSubStream(seed)
# Define the "Next Element" function for the iterator
nextEl <- function() (seed <<- nextRNGSubStream(seed))
obj <- list(nextElem=nextEl)
class(obj) <- c('abstractiter', 'iter')
obj
}
convseed <- function(iseed) {
saveseed <- if (exists('.Random.seed', where=.GlobalEnv, inherits=FALSE))
get('.Random.seed', pos=.GlobalEnv, inherits=FALSE)
saverng <- RNGkind("L'Ecuyer-CMRG")
tryCatch({
set.seed(iseed)
get('.Random.seed', pos=.GlobalEnv, inherits=FALSE)
},
finally={
RNGkind(saverng[1], saverng[2])
if (is.null(saveseed))
rm('.Random.seed', pos=.GlobalEnv)
else
assign('.Random.seed', saveseed, pos=.GlobalEnv)
})
}
|
