Loop over a Foreach Expression using Futures
foreach %dofuture% exprA foreach object created by foreach::foreach()
and foreach::times().
An R expression.
The value of the foreach call.
This is a replacement for %dopar% of the foreach package
that leverages the future framework.
When using %dofuture%:
there is no need to use registerDoFuture()
there is no need to use %dorng% of the doRNG package
(but you need to specify .options.future = list(seed = TRUE)
whenever using random numbers in the expr expression)
global variables and packages are identified automatically by the future framework
errors are relayed as-is (with %dopar% they captured and modified)
When using %dofuture%, the future framework identifies globals and
packages automatically (via static code inspection). However, there
are cases where it fails to find some of the globals or packages. When
this happens, one can specify the future::future() arguments globals
and packages via foreach argument .options.future. For example,
if you specify argument
.options.future = list(globals = structure(TRUE, ignore = "b", add = "a"))
then globals are automatically identified (TRUE), but it ignores b and
always adds a.
An alternative to specifying the globals and the packages options via
.options.future, is to use the %globals%
and the %packages% operators.
See the examples for an illustration.
For further details and instructions, see future::future().
The %dofuture% uses the future ecosystem to generate proper random
numbers in parallel in the same way they are generated in, for instance,
future.apply. For this to work, you need to specify
.options.future = list(seed = TRUE). For example,
y <- foreach(i = 1:3, .options.future = list(seed = TRUE)) %dofuture% {
rnorm(1)
}Unless seed is FALSE or NULL, this guarantees that the exact same
sequence of random numbers are generated given the same initial
seed / RNG state - this regardless of type of future backend, number of
workers, and scheduling ("chunking") strategy.
RNG reproducibility is achieved by pregenerating the random seeds for all iterations by using L'Ecuyer-CMRG RNG streams. In each iteration, these seeds are set before evaluating the foreach expression. Note, for large number of iterations this may introduce a large overhead.
If seed = TRUE, then .Random.seed
is used if it holds a L'Ecuyer-CMRG RNG seed, otherwise one is created
randomly.
If seed = FALSE, it is expected that none of the foreach iterations
use random number generation.
If they do, then an informative warning or error is produces depending
on settings. See future::future for more details.
Using seed = NULL, is like seed = FALSE but without the check
whether random numbers were generated or not.
As input, seed may also take a fixed initial seed (integer),
either as a full L'Ecuyer-CMRG RNG seed (vector of 1+6 integers), or
as a seed generating such a full L'Ecuyer-CMRG seed. This seed will
be used to generated one L'Ecuyer-CMRG RNG stream for each iteration.
An alternative to specifying the seed option via .options.future,
is to use the %seed% operator. See
the examples for an illustration.
For further details and instructions, see
future.apply::future_lapply().
Whether load balancing ("chunking") should take place or not can be
controlled by specifying either argument
.options.future = list(scheduling = <ratio>) or
.options.future = list(chunk.size = <count>) to foreach().
The value chunk.size specifies the average number of elements
processed per future ("chunks").
If +Inf, then all elements are processed in a single future (one worker).
If NULL, then argument future.scheduling is used.
The value scheduling specifies the average number of futures
("chunks") that each worker processes.
If 0.0, then a single future is used to process all iterations;
none of the other workers are not used.
If 1.0 or TRUE, then one future per worker is used.
If 2.0, then each worker will process two futures (if there are
enough iterations).
If +Inf or FALSE, then one future per iteration is used.
The default value is scheduling = 1.0.
For further details and instructions, see
future.apply::future_lapply().
Attribute ordering of chunk.size or scheduling can be used to
control the ordering the elements are iterated over, which only affects
the processing order and not the order values are returned.
This attribute can take the following values:
index vector - an numeric vector of length nX.
function - an function taking one argument which is called as
ordering(nX) and which must return an
index vector of length nX, e.g.
function(n) rev(seq_len(n)) for reverse ordering.
"random" - this will randomize the ordering via random index
vector sample.int(nX).
where nX is the number of foreach iterations to be done.
For example,
.options.future = list(scheduling = structure(2.0, ordering = "random")).
Note, when elements are processed out of order, then captured standard output and conditions are also relayed in that order, that is, out of order.
For further details and instructions, see
future.apply::future_lapply().
# \donttest{
plan(multisession) # parallelize futures on the local machine
y <- foreach(x = 1:10, .combine = rbind) %dofuture% {
y <- sqrt(x)
data.frame(x = x, y = y, pid = Sys.getpid())
}
print(y)
#> x y pid
#> 1 1 1.000000 1019182
#> 2 2 1.414214 1019182
#> 3 3 1.732051 1019181
#> 4 4 2.000000 1019180
#> 5 5 2.236068 1019183
#> 6 6 2.449490 1019185
#> 7 7 2.645751 1019187
#> 8 8 2.828427 1019186
#> 9 9 3.000000 1019184
#> 10 10 3.162278 1019184
## Random number generation
y <- foreach(i = 1:3, .combine = rbind, .options.future = list(seed = TRUE)) %dofuture% {
data.frame(i = i, random = runif(n = 1L))
}
print(y)
#> i random
#> 1 1 0.5805619
#> 2 2 0.3676488
#> 3 3 0.4598859
## Random number generation (alternative specification)
y <- foreach(i = 1:3, .combine = rbind) %dofuture% {
data.frame(i = i, random = runif(n = 1L))
} %seed% TRUE
print(y)
#> i random
#> 1 1 0.68864083
#> 2 2 0.22946006
#> 3 3 0.07882239
## Random number generation with the foreach() %:% nested operator
y <- foreach(i = 1:3, .combine = rbind) %:%
foreach(j = 3:5, .combine = rbind, .options.future = list(seed = TRUE)) %dofuture% {
data.frame(i = i, j = j, random = runif(n = 1L))
}
print(y)
#> i j random
#> 1 1 3 0.13393400
#> 2 1 4 0.35087825
#> 3 1 5 0.54217272
#> 4 2 3 0.56372082
#> 5 2 4 0.10511164
#> 6 2 5 0.58537767
#> 7 3 3 0.05640405
#> 8 3 4 0.03938918
#> 9 3 5 0.56056544
## Random number generation with the nested foreach() calls
y <- foreach(i = 1:3, .combine = rbind, .options.future = list(seed = TRUE)) %dofuture% {
foreach(j = 3:5, .combine = rbind, .options.future = list(seed = TRUE)) %dofuture% {
data.frame(i = i, j = j, random = runif(n = 1L))
}
}
print(y)
#> i j random
#> 1 1 3 0.6587098
#> 2 1 4 0.1306959
#> 3 1 5 0.2513151
#> 4 2 3 0.1306959
#> 5 2 4 0.2513151
#> 6 2 5 0.5332196
#> 7 3 3 0.2513151
#> 8 3 4 0.5332196
#> 9 3 5 0.7158639
# }
# \dontshow{
## R CMD check: make sure any open connections are closed afterward
if (!inherits(plan(), "sequential")) plan(sequential)
# }