并行化采用外部指针的函数 (XPtr)

这个问题既不是重复的this one https://stackoverflow.com/questions/37167479也不属于this one https://stackoverflow.com/questions/6074310，这是关于返回外部指针的函数。

问题就在这里。此后的 Rcpp 代码定义了两个函数，一个用于创建 XPtr，另一个可以在 XPtr 上工作。

#include <Rcpp.h>
using namespace Rcpp;

//[[Rcpp::export]]
SEXP f(int n) {
  std::vector<int> * v = new std::vector<int>;

  for(int i = 0; i < n; i++)
    v->push_back(i);

  XPtr< std::vector<int> > p(v, true);
  return p;
}

//[[Rcpp::export]]
int g(XPtr< std::vector<int> > p, int i) {
  return (*p)[i];

而且效果很好：

> x <- f(100)
> g(x, 45)
[1] 45

让我们尝试并行调用g。这有效：

require(parallel)
test1 <- function(a) {
  cl <- makeForkCluster(nnodes=2)
  r <- parLapply(cl, 1:5, function(i) g(a,i) )
  stopCluster(cl)
  return(r)
}

预期行为：

> unlist( test1(x) )
[1] 1 2 3 4 5

但这不起作用：

test2 <- function(a) {
  cl <- makeForkCluster(nnodes=2)

  p <- g(a, 0)
  r <- parLapply(cl, 1:5, function(i) g(a,i) )
  stopCluster(cl)
  return(r)
}

意外行为：

> test2(x)
Error in checkForRemoteErrors(val) : 
  2 nodes produced errors; first error: external pointer is not valid

问题似乎是由于在调用集群中的从属设备之前在函数中使用了一次外部指针而引起的。如何解释这种行为？有解决方法吗？提前谢谢了。

在你的函数开始时，a是一个承诺，即表示在特定环境中评估特定表达式。当您访问变量时，将计算表达式，所以现在a是一个指针，并且该指针特定于特定的 R 实例。您可以使用以下命令查看此内容pryr::promise_info:

test2 <- function(a) {
  cl <- makeForkCluster(nnodes = 2)
  print(pryr::promise_info(a))
  p <- g(a, 0)
  print(pryr::promise_info(a))
  stopCluster(cl)
  return(r)
}

Output:

$code
x

$env
<environment: R_GlobalEnv>

$evaled
[1] FALSE

$value
NULL

$code
x

$env
NULL

$evaled
[1] TRUE

$value
<pointer: 0x565295e3a410>

一种解决方法是使用eval(substitute(a)):

test2 <- function(a) {
  cl <- makeForkCluster(nnodes = 2)
  print(pryr::promise_info(a))
  p <- g(eval(substitute(a)), 0)
  print(pryr::promise_info(a))
  r <- parLapply(cl, 1:5, function(i) g(a,i) )
  stopCluster(cl)
  return(r)
}

我确信有更好的方法。非标准的评价对我来说还是有点陌生​​……