RcppArmadillo sample.h编译错误

1 回答

• 1

  此功能适用于 arma::vec 或 NumericVector 中的预子集数据 . 始终具有且始终为 . 不要将此与从子集操作（例如 .col() ， .cols() 或 .submat() ）获得的中间向量一起使用 .

  ---

  您遇到的问题是've decided to subset the data within the call to sample. (You'已省略了诊断此部分的代码，因此我在此推测 . ）由于 sample() 需要同时使用Rcpp和Armadillo数据类型，因此从未调用过Armadillo特定大小的成员函数 . 相反，库选择调用 STL 容器的 .size() 成员函数，armadillo supported，因为它在两个对象之间共享 . 但是，犰狳限制了成员函数在"active"数据结构中的实现位置，而不是临时数据 . 因此， .size() 成员函数未实现 subview_col . 所以，我们最终得到的错误是：

  错误：'const struct arma :: subview_col'没有名为'size'的成员

  ---

  要解决此限制并节省内存，请使用将重用内存的advanced vec ctor，从而避免创建中间件 arma::subview_col .

  // [[Rcpp::depends(RcppArmadillo)]]
  #include <RcppArmadillo.h>
  
  #include <RcppArmadilloExtensions/sample.h>
  
  // [[Rcpp::export]]
  void adv_rnd(int nrow, int ncol, bool replace = true){
      // Create a matrix of given dimensions
      arma::mat X(nrow, ncol);
      X.randn();
  
      // Show state before randomization
      Rcpp::Rcout << "Before Randomization:" << std::endl << X << std::endl;
  
      // Randomize each column
      for(int i = 0; i < ncol; ++i){
          arma::vec Y(X.colptr(i), nrow, false, true);
          X.col(i) = Rcpp::RcppArmadillo::sample(Y, nrow, replace);
      }
  
      // Show state after randomization
      Rcpp::Rcout << "After Randomization:" << std::endl << X << std::endl;
  
  }
  

  样本输出：

  > adv_rnd(3,3)
  Before Randomization:
    -0.7197   1.2590  -0.5898
     0.0253   0.1493  -0.0685
    -0.6074   1.3843   0.0400
  
  After Randomization:
    -0.7197   1.2590   0.0400
    -0.6074   1.2590  -0.5898
    -0.6074   0.1493  -0.0685
  

  回复于 2024-04-19T07:53:21+08:00