在 Foreach 循环中修改变量？

Question

我正在尝试构建一个并行的 for-each 循环来修改两个矩阵。我在开始循环之前创建并初始化矩阵。这是一个虚拟程序，它演示了我的代码应该如何工作。程序执行没有错误，但在 foreach 循环完成后矩阵仍然为空。注意：这是我的代码的简化版本，而不是实际代码本身。

#Assume that I've loaded parallel and doParallel and that my computer has 2+ cores 
cluster <- makeCluster(detectCores())
registerDoParallel(cluster)
a1 <- array(dim=c(9,9))
a2 <- array(dim=c(9,9))
numbers <- 1:18
foreach (i=1:9, .combine='c') %dopar%{
    a1[i,] = numbers[1:9]
    a2[i,] = numbers[10:18]
}
stopCluster(cluster)

为什么这个程序不填充 a1 和 a2 的行？

编辑：我在之前的任何帖子中都找不到我的问题的答案。 C# 和 Perl 也有类似的线程，但没有一个与 R 相关。

Answer 1

您不能使用 foreach 来执行这种副作用。循环中的分配发生在工作进程中，并且故意不发送回主进程。你执行 foreach 是为了它的返回值，所以如果你想生成两个矩阵，foreach 应该返回一个包含两个矩阵的列表。

这是一个示例，它使用组合函数从工作人员生成的行列表中构造多个矩阵：

rcomb <- function(...) {
  args <- list(...)
  lapply(seq_along(args[[1]]), function(i)
         do.call('rbind', lapply(args, function(a) a[[i]])))
}

numbers <- 1:18
m <- foreach (i=1:9, .combine='rcomb', .multicombine=TRUE) %dopar% {
    list(numbers[1:9], numbers[10:18])
}

虽然rcomb 有点复杂，但它可以处理任意数量的矩阵，并且可以调用许多任务结果，这对效率来说很重要。

Answer 2

也许这会有所帮助。无法测试，请见谅。

 require(snowfall)

sfInit(cpus=2,parallel=TRUE)

a1 <- matrix(NA,ncol=9,nrow=0)
a2 <- matrix(NA,ncol=9,nrow=0)
numbers <- 1:18

a1 = rbind(a1,sfSapply(1:9,function(x) return(numbers[1:9]), numbers))
a2 = rbind(a2,sfSapply(1:9,function(x) return(numbers[10:18]), numbers))

sfStop()