允许在 DEoptim for R 中使用 Rcpp 函数

Question

一段时间以来，我一直在使用 DEoptim 来测试水文算法的不同参数。 代码主要是用R写的，但是有一个函数是用Rcpp写的。 如果我在非并行模式下运行 DEoptim，它运行良好，但如果我在并行模式下运行（即 paralleltype=1），代码返回一个错误，说它找不到我的 Rcpp 函数。 所以 Rcpp 函数看起来像这样：

loadcppfunctions <- function() {
eastfunc <<- 'NumericMatrix eastC(NumericMatrix e, NumericMatrix zerocolmatrix, NumericMatrix zerorowmatrix) {
int ecoln = e.ncol();
int ecolnlessone = ecoln - 1;
int erown = e.nrow();
int erownlessone = erown - 1;

NumericMatrix eout(e.nrow(),e.ncol()) ;
for (int j = 0;j < ecoln;j++) {
if (j > 0) {
eout(_,j) = e(_,j-1);
} else {
eout(_,j) = e(_,0);
}
}
eout(_,0) = zerocolmatrix(_,0);
return eout;
}'
eastC <<- cppFunction(eastfunc)
}

然后我只使用：

loadcppfunctions()

稍后在代码中我调用这个函数如下：

movefdrerunoff <- eastC(fdrerunoff, zerocolmatrix, zerorowmatrix)

正如我所说，一切正常 - 但如果我按如下方式运行 DEoptim：

ans <- DEoptimone(Calibrate,lower,upper,DEoptim.control(trace=TRUE,parallelType=1,parVar=c(parVarnames),packages=c("raster","rgdal","maptools","matrixcalc","Rcpp","RcppArmadillo")))

它没有说：

Error in checkForRemoteErrors(val) : 
  7 nodes produced errors; first error: could not find function "eastC"

那么当所有其他基于 R 的函数都很好时，我如何让 DEoptim 看到这个函数。

谢谢，安东尼沃克

Answer 1

我发现通过在主 DEoptim 函数（校准）中添加 Rcpp 函数起作用了。 校准功能如下所示：

Calibrate <- function(x) {
eastfunc <<- 'NumericMatrix eastC(NumericMatrix e, NumericMatrix zerocolmatrix, NumericMatrix zerorowmatrix) {
int ecoln = e.ncol();
int ecolnlessone = ecoln - 1;
int erown = e.nrow();
int erownlessone = erown - 1;

NumericMatrix eout(e.nrow(),e.ncol()) ;
for (int j = 0;j < ecoln;j++) {
if (j > 0) {
eout(_,j) = e(_,j-1);
} else {
eout(_,j) = e(_,0);
}
}
eout(_,0) = zerocolmatrix(_,0);
return eout;
}'
eastC <<- cppFunction(eastfunc)

cmax <<- x[1]
Cr <<- x[2]
Cl <<- x[3]
Crb <<- x[4]
Clb <<- x[5]
returnflowriver <<- x[6]
returnflowland <<- x[7]
kd <<- x[8]
startyear()
-NashSutcliffe
}

然后将 DEoptim 运行为：

ans <- DEoptimone(Calibrate,lower,upper,DEoptim.control(trace=TRUE,parallelType=1,parVar=c(parVarnames),packages=c("raster","rgdal","maptools","matrixcalc","Rcpp","RcppArmadillo","moveCpp")))