R Under development (unstable) (2023-12-20 r85711 ucrt) -- "Unsuffered Consequences" Copyright (C) 2023 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > source("incl/start.R") Loading required package: future [13:13:49.032] plan(): Setting new future strategy stack: [13:13:49.033] List of future strategies: [13:13:49.033] 1. sequential: [13:13:49.033] - args: function (..., envir = parent.frame(), workers = "") [13:13:49.033] - tweaked: FALSE [13:13:49.033] - call: future::plan("sequential") [13:13:49.048] plan(): nbrOfWorkers() = 1 > library("tools") ## toTitleCase() > > message("*** future_mapply() - globals ...") *** future_mapply() - globals ... > > #plan(cluster, workers = "localhost") > plan(sequential) [13:13:49.056] plan(): Setting new future strategy stack: [13:13:49.056] List of future strategies: [13:13:49.056] 1. sequential: [13:13:49.056] - args: function (..., envir = parent.frame(), workers = "") [13:13:49.056] - tweaked: FALSE [13:13:49.056] - call: plan(sequential) [13:13:49.070] plan(): nbrOfWorkers() = 1 > > options(future.debug = FALSE) > a <- 1 > b <- 2 > > globals_set <- list( + A = FALSE, + B = TRUE, + C = c("a", "b"), + D = list(a = 1, b = 2) + ) > > x <- list(1) > > for (name in names(globals_set)) { + globals <- globals_set[[name]] + message("Globals set ", sQuote(name)) + str(globals) + + y_truth <- tryCatch({ + mapply(function(x) { + median(c(x, a, b)) + }, x) + }, error = identity) + + y1 <- tryCatch({ + future_mapply(function(x) { + median(c(x, a, b)) + }, x) + }, error = identity) + print(y1) + stopifnot(identical(y1, y_truth)) + + y2 <- tryCatch({ + future_mapply(function(x) { + median(c(x, a, b)) + }, x, future.globals = globals, future.packages = "utils") + }, error = identity) + print(y2) + stopifnot(identical(y2, y_truth)) + } Globals set 'A' logi FALSE [1] 1 [1] 1 Globals set 'B' logi TRUE [1] 1 [1] 1 Globals set 'C' chr [1:2] "a" "b" [1] 1 [1] 1 Globals set 'D' List of 2 $ a: num 1 $ b: num 2 [1] 1 [1] 1 > > message("*** future_mapply() - globals ... DONE") *** future_mapply() - globals ... DONE > > > > ## Test adopted from http://stackoverflow.com/questions/42561088/nested-do-call-within-a-foreach-dopar-environment-cant-find-function-passed-w > > message("*** future_mapply() - tricky globals ...") *** future_mapply() - tricky globals ... > > my_add <- function(a, b) a + b > > call_my_add <- function(a, b) { + do.call(my_add, args = list(a = a, b = b)) + } > > call_my_add_caller <- function(a, b, FUN = call_my_add) { + do.call(FUN, args = list(a = a, b = b)) + } > > main <- function(x = 1:2, caller = call_my_add_caller, + args = list(FUN = call_my_add)) { + results <- future_mapply(function(i) { + do.call(caller, args = c(list(a = i, b = i + 1L), args)) + }, x) + results + } > > x <- list(list(1:2)) > z_length <- mapply(do.call, args = x, MoreArgs = list(what = length)) > fun <- function(...) sum(...) > z_fun <- mapply(do.call, args = x, MoreArgs = list(what = fun)) > > y0 <- NULL > for (strategy in supportedStrategies()) { + plan(strategy) + + y <- main(1:3) + if (is.null(y0)) y0 <- y + stopifnot(identical(y, y0)) + + message("- future_mapply(do.call, x, ...) ...") + z <- future_mapply(do.call, args = x, MoreArgs = list(what = length)) + stopifnot(identical(z, z_length)) + z <- future_mapply(do.call, args = x, MoreArgs = list(what = fun)) + stopifnot(identical(z, z_fun)) + + message("- future_mapply(FUN, x, ...) - passing arguments via '...' ...") + ## typeof() == "list" + obj <- data.frame(a = 1:2) + stopifnot(typeof(obj) == "list") + y <- future_mapply(function(a, b) typeof(b), 1L, MoreArgs = list(b = obj)) + stopifnot(identical(y[[1]], typeof(obj))) + + ## typeof() == "environment" + obj <- new.env() + stopifnot(typeof(obj) == "environment") + y <- future_mapply(function(a, b) typeof(b), 1L, MoreArgs = list(b = obj)) + stopifnot(identical(y[[1]], typeof(obj))) + + ## typeof() == "S4" + if (requireNamespace("methods")) { + obj <- methods::getClass("MethodDefinition") + stopifnot(typeof(obj) == "S4") + y <- future_mapply(function(a, b) typeof(b), 1L, MoreArgs = list(b = obj)) + stopifnot(identical(y[[1]], typeof(obj))) + } + + message("- future_mapply(FUN, X, ...) - 'X' containing globals ...") + ## From https://github.com/HenrikBengtsson/future.apply/issues/12 + a <- 42 + b <- 21 + X <- list( + function(b) 2 * a, + function() b / 2, + function() a + b, + function() nchar(toTitleCase("hello world")) + ) + z0 <- mapply(function(s, f) f() + s, s = seq_along(X), X) + str(z0) + z1 <- future_mapply(function(s, f) f() + s, s = seq_along(X), X) + str(z1) + stopifnot(identical(z1, z0)) + } - future_mapply(do.call, x, ...) ... - future_mapply(FUN, x, ...) - passing arguments via '...' ... - future_mapply(FUN, X, ...) - 'X' containing globals ... num [1:4] 85 12.5 66 15 num [1:4] 85 12.5 66 15 - future_mapply(do.call, x, ...) ... - future_mapply(FUN, x, ...) - passing arguments via '...' ... - future_mapply(FUN, X, ...) - 'X' containing globals ... num [1:4] 85 12.5 66 15 num [1:4] 85 12.5 66 15 > > message("*** future_mapply() - tricky globals ... DONE") *** future_mapply() - tricky globals ... DONE > > > message("*** future_mapply() - missing arguments ...") *** future_mapply() - missing arguments ... > > ## Here 'abc' becomes missing, i.e. missing(abc) is TRUE > foo <- function(x, abc) mapply(function(y) y, x) > y0 <- foo(1:2) > foo <- function(x, abc) future_mapply(function(y) y, x) > y <- foo(1:2) > stopifnot(identical(y, y0)) > > message("*** future_mapply() - missing arguments ... DONE") *** future_mapply() - missing arguments ... DONE > > > message("*** future_mapply() - false positives ...") *** future_mapply() - false positives ... > > ## Here 'abc' becomes a promise, which fails to resolve > ## iff 'xyz' does not exist. (Issue #161) > suppressWarnings(rm(list = "xyz")) > foo <- function(x, abc) mapply(function(y) y, x) > y0 <- foo(1:2, abc = (xyz >= 3.14)) > foo <- function(x, abc) future_mapply(function(y) y, x) > y <- foo(1:2, abc = (xyz >= 3.14)) > stopifnot(identical(y, y0)) > > message("*** future_mapply() - false positives ... DONE") *** future_mapply() - false positives ... DONE > > > message("*** future_mapply() - too large ...") *** future_mapply() - too large ... > > X <- replicate(10L, 1:100, simplify = FALSE) > FUN <- function(x) { + getOption("future.globals.maxSize") + } > > y0 <- mapply(FUN = FUN, X) > > sizes <- unclass(c(FUN = object.size(FUN), X = object.size(X))) > cat(sprintf("Baseline size of globals: %.2f KiB\n", sizes[["FUN"]] / 1024)) Baseline size of globals: 1.80 KiB > > message("- true positive ...") - true positive ... > oMaxSize <- getOption("future.globals.maxSize") > options(future.globals.maxSize = 1L) > res <- tryCatch({ + y <- future_mapply(FUN = FUN, X) + }, error = identity) > stopifnot(inherits(res, "error")) > res <- NULL > options(future.globals.maxSize = oMaxSize) > > maxSize <- getOption("future.globals.maxSize") > y <- future_mapply(FUN = FUN, X) > str(y) List of 10 $ : NULL $ : NULL $ : NULL $ : NULL $ : NULL $ : NULL $ : NULL $ : NULL $ : NULL $ : NULL > stopifnot(all(sapply(y, FUN = identical, oMaxSize))) > > message("- approximately invariant to chunk size ...") - approximately invariant to chunk size ... > maxSize <- sizes[["FUN"]] + sizes[["X"]] / length(X) > maxSize <- 3.0 * maxSize ## Add a bit of leeway > options(future.globals.maxSize = maxSize) > > for (chunk.size in c(1L, 2L, 5L, structure(10L, ordering = "random"))) { + y <- future_mapply(FUN = FUN, X, future.chunk.size = chunk.size) + str(y) + stopifnot(all(unlist(y) == maxSize)) + cat(sprintf("maxSize = %g bytes\nfuture.globals.maxSize = %g bytes\n", + maxSize, getOption("future.globals.maxSize"))) + stopifnot(getOption("future.globals.maxSize") == maxSize) + } num [1:10] 6917 6917 6917 6917 6917 ... maxSize = 6916.8 bytes future.globals.maxSize = 6916.8 bytes num [1:10] 6917 6917 6917 6917 6917 ... maxSize = 6916.8 bytes future.globals.maxSize = 6916.8 bytes num [1:10] 6917 6917 6917 6917 6917 ... maxSize = 6916.8 bytes future.globals.maxSize = 6916.8 bytes num [1:10] 6917 6917 6917 6917 6917 ... maxSize = 6916.8 bytes future.globals.maxSize = 6916.8 bytes > y <- NULL > options(future.globals.maxSize = oMaxSize) > > message("*** future_mapply() - too large ... DONE") *** future_mapply() - too large ... DONE > > > message("*** future_mapply() - globals exceptions ...") *** future_mapply() - globals exceptions ... > > res <- tryCatch({ + y <- future_mapply(function(x) x, 1, future.globals = 42) + }, error = identity) > stopifnot(inherits(res, "error")) > > res <- tryCatch({ + y <- future_mapply(function(x) x, 1, future.globals = list(1)) + }, error = identity) > stopifnot(inherits(res, "error")) > > res <- tryCatch({ + y <- future_mapply(function(x) x, 1, future.globals = "...future.FUN") + }, error = identity) > stopifnot(inherits(res, "error")) > > ...future.elements_ii <- 42L > X <- list(function() 2 * ...future.elements_ii) > res <- tryCatch({ + y <- future_mapply(FUN = function(f) f(), X) + }, error = identity) > stopifnot(inherits(res, "error")) > > message("*** future_mapply() - globals exceptions ... DONE") *** future_mapply() - globals exceptions ... DONE > > source("incl/end.R") > > proc.time() user system elapsed 1.23 0.12 2.54