r - foreach не может найти объект внутри функции

Question

У меня есть определенная ниже функция с именем algor, которая переводится из MATLAB в R. Чтобы сделать функцию быстрее, я впервые использую конструкцию foreach. У меня есть полный код функции ниже:

library("ramify")
library("foreach")

algor <- function (vc) {

# initialize A, ybar, and Ia
A <- 0

ybar <- 0

Ia <- 0

# x is the first column of vc   
x <- vc[, 1, drop = FALSE]

# y is the second column of vc  
y <- vc[, 2, drop = FALSE]

# n is the length of x
n <- length(x)

foreach(i = 1:(n-1), .combine = 'c', .export = c("A", "ybar", "Ia", "x", "y")) %do% {
        A <- A + 0.5 * (x[i] - x[i+1]) * (y[i] + y[i+1])
        ybar <- ybar + (1 / 6) * (x[i] - x[i+1]) * (y[i] ^ 2 + y[i] * y[i+1] + y[i+1] ^ 2)
        Ia <- Ia + (1 / 12) * (x[i] - x[i+1]) * (y[i] ^ 3 + y[i] ^ 2 * y[i+1] + y[i] * y[i+1] ^ 2 + y[i+1] ^ 3)
}

props <- mat("A, Ia, ybar", eval = TRUE)

return(props)
}


inner <- mat("0, 300; 300, 300; 300, 695; 0, 695; 0, 300")

algor(inner)

Хотя я экспортировал A, ybar, Ia, x и y, я получаю сообщение об ошибке, что объект A не найден, что ниже:

Error in eval(parse(text = paste0("c(", paste0(char_vals, collapse = ","),  : 
  object 'A' not found
Called from: eval(parse(text = paste0("c(", paste0(char_vals, collapse = ","), 
    ")")))

Как мне заставить foreach распознать определенные объекты: A, ybar, Ia, x и y?

Спасибо.

Answer 1

Попробуйте определить переменные .GlobalEnv в цикле foreach в каждом вызове.

library("ramify")
library("foreach")

algor <- function (vc) {

  # initialize A, ybar, and Ia
  A <- 0

  ybar <- 0

  Ia <- 0

  # x is the first column of vc   
  x <- vc[, 1, drop = FALSE]

  # y is the second column of vc  
  y <- vc[, 2, drop = FALSE]

  # n is the length of x
  n <- length(x)

  foreach(i = 1:(n-1), .combine = 'c', .export = c("A", "ybar", "Ia", "x", "y")) %do% {
    .GlobalEnv$A <- A
    .GlobalEnv$ybar <- ybar
    .GlobalEnv$Ia <- Ia
    A <- A + 0.5 * (x[i] - x[i+1]) * (y[i] + y[i+1])
    ybar <- ybar + (1 / 6) * (x[i] - x[i+1]) * (y[i] ^ 2 + y[i] * y[i+1] + y[i+1] ^ 2)
    Ia <- Ia + (1 / 12) * (x[i] - x[i+1]) * (y[i] ^ 3 + y[i] ^ 2 * y[i+1] + y[i] * y[i+1] ^ 2 + y[i+1] ^ 3)
  }

  props <- mat("A, Ia, ybar", eval = TRUE)

  return(props)
}


inner <- mat("0, 300; 300, 300; 300, 695; 0, 695; 0, 300")

algor(inner)

Возвращает:

       [,1]        [,2]     [,3]
[1,] 118500 30870237500 58953750