1

次のコードを使用します。

  library(foreach)
  library(doParallel)
  N<-5
  cl<-makeCluster(8)
  registerDoParallel(cl)

  #loop
  s8 <- foreach(i=1:N, .combine='rbind') %:%
    foreach(j=1:N, .combine='c') %dopar% {

        dis <-as.numeric (i+j)  } ## In reality there something more complicated

  stopCluster(cl)

次の結果が得られます。

> s8
         [,1] [,2] [,3] [,4] [,5]
result.1    2    3    4    5    6
result.2    3    4    5    6    7
result.3    4    5    6    7    8
result.4    5    6    7    8    9
result.5    6    7    8    9   10

次のように上対角行列を取得したいと思います。

> s8
         [,1] [,2] [,3] [,4] [,5]
result.1    2    3    4    5    6
result.2    0    4    5    6    7
result.3    0    0    6    7    8
result.4    0    0    0    8    9
result.5    0    0    0    0   10

内側のループのイテレータを

foreach(j=i:N, .combine='c') %dopar% {

必要な結果が得られません。

4

1 に答える 1

1

たぶんこれです（ただし、おそらく外側のループのみを並列化します）：

library(foreach)
library(doParallel)
N<-5
cl<-makeCluster(4)
registerDoParallel(cl)

#loop
#result is a list of vectors
s8 <- foreach(i=1:N) %:%
  foreach(j=i:N, .combine='c') %dopar% {

    as.numeric (i+j)  } ## In reality there something more complicated

stopCluster(cl)

同じ長さを取得するための後処理:

s8 <- t(vapply(s8, function(x) {
  x <- rev(x)
  length(x) <- N
  rev(x)
}, FUN.VALUE = numeric(N)))

s8[lower.tri(s8)] <- 0

#     [,1] [,2] [,3] [,4] [,5]
#[1,]    2    3    4    5    6
#[2,]    0    4    5    6    7
#[3,]    0    0    6    7    8
#[4,]    0    0    0    8    9
#[5,]    0    0    0    0   10
于 2016-04-25T15:16:48.447 に答える