使用Rcpp :: Nullable与RcppArmadillo类型错误

Question

我在R包中使用RcppArmadillo，我想在参数列表中使用Rcpp :: Nullable。

NumericVector d_snb(NumericVector& x, 
       Nullable<arma::mat> size_param = R_NilValue, 
       const int& infinite = 100000, const bool& log_v = false) 

这给了如下错误：

Error: package or namespace load failed for ‘packagex’ in dyn.load(file, DLLpath = DLLpath, ...): unable to load shared object '/usr/local/lib/R/3.4/site-library/packagex/libs/packagex.so': dlopen(/usr/local/lib/R/3.4/site-library/packagex/libs/packagex.so, 6): 
Symbol not found: __Z5d_snbRN4Rcpp6VectorILi14ENS_15PreserveStorageEEENS_8NullableIS2_EES5_S5_RKiRKb 
Referenced from: /usr/local/lib/R/3.4/site-library/packagex/libs/packagex.so 
Expected in: flat namespace in /usr/local/lib/R/3.4/site-library/packagex/libs/packagex.so 
Error: loading failed 
Execution halted 

惟一可行的办法，现在似乎对参数作为NumericVector然后获取内容，然后将它扔在犰狳的类型。

NumericVector d_snb(NumericVector& x, 
       Nullable<NumericVector> size_param = R_NilValue ...) 
{ 
... 

    if(size_param.isNotNull()) { 
    arma::mat test(NumericVector(size_param)); 
    param_check++; 
    } 
... 
} 

这给出了一个警告：

d_snb.cpp:36:21: warning: parentheses were disambiguated as a function declaration [-Wvexing-parse] 
    arma::mat test(NumericVector(size_param)); 
       ^~~~~~~~~~~~~~~~~~~~~~~~~~~ 
d_snb.cpp:36:22: note: add a pair of parentheses to declare a variable 
    arma::mat test(NumericVector(size_param)); 
       ^
       (      ) 
1 warning generated. 

什么是去了解它的最好方法？

Answer 1

是的，Rcpp::Nullable<>只适用于基于SEXP的Rcpp类型。

所以你必须做你发现的两步过程。但我认为你可以做一个简单：arma::mat test = Rcpp::as<arma::mat>(size_param);

这里是一个编译一个完整的例子（但“无”）：

#include <RcppArmadillo.h> 
// [[Rcpp::depends(RcppArmadillo)]] 

// [[Rcpp::export]] 
void d_snb(Rcpp::NumericVector& x, 
      Rcpp::Nullable<Rcpp::NumericVector> size_param = R_NilValue) { 

    if (size_param.isNotNull()) { 
    arma::vec test = Rcpp::as<arma::vec>(size_param); 
    test.print("vector"); 
    } 
    Rcpp::Rcout << x << std::endl; 

} 

和快速演示：

R> sourceCpp("/tmp/so44102346.cpp") 
R> d_snb(c(1,2,3)) 
1 2 3 
R> d_snb(c(1,2,3), c(4,5,6)) 
vector 
    4.0000 
    5.0000 
    6.0000 
1 2 3 
R>