将 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 获取，然后获取内容，然后将其转换为 Armadillo 类型。

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);

这是一个完整的编译示例（但 'nothing'）：

#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>

将 Rcpp::Nullable 与 RcppArmadillo 类型一起使用时出错

Error using Rcpp::Nullable with RcppArmadillo types

rcpp