这个问题既不是this one也不this one的重复,这是有关返回外部指针的函数。
这里的问题。在下文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)
}
输出:
$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)
}
我相信有更好的方法。非标评价仍然是一个有点陌生的我...