5

我主要使用 R，但最终想使用 Rcpp 与一些接收和返回二维数值数组的 C++ 函数交互。因此，要开始使用 C++ 和 Rcpp，我想我只需创建一个小函数，将我的可变长度数值向量的 R 列表转换为 C++ 等效项，然后再返回。

require(inline)
require(Rcpp)

test1 = cxxfunction(signature(x='List'), body = 
'
  using namespace std;
  List xlist(x);
  int xlen = xlist.size();
  vector< vector<int> > xx;
  for(int i=0; i<xlen; i++) {
    vector<int> test = as<vector<int> > (xlist[i]);
    xx.push_back(test);
  }
  return(wrap(xx));
'
, plugin='Rcpp')

这就像我期望的那样工作：

> test1(list(1:2, 4:6))
[[1]]
[1] 1 2

[[2]]
[1] 4 5 6

诚然，我只是完成了非常详尽的文档的一部分，但是有没有比使用 for 循环更好（即更像 Rcpp）的方法来进行 R -> C++ 转换？我认为可能不是，因为文档提到（至少使用内置方法）as“提供的灵活性较低，并且目前处理将 R 对象转换为原始类型”，但我想检查一下，因为我是个新手在这个区域。

4

1 回答 1

7

我会给你一个可重复的例子的奖励积分，当然还有使用 Rcpp :) 然后我会因为没有在 rcpp-devel 列表上询问而把它们拿走......

至于转换 STL 类型：您不必这样做，但是当您决定这样做时，as<>()成语是正确的。我能想到的唯一“更好的方法”是像在 R 本身中那样进行名称查找：

require(inline)
require(Rcpp)

set.seed(42)
xl <- list(U=runif(4), N=rnorm(4), T2df=rt(4,2))

fun <- cxxfunction(signature(x="list"), plugin="Rcpp", body = '
  Rcpp::List xl(x);         
  std::vector<double> u = Rcpp::as<std::vector<double> >(xl["U"]);
  std::vector<double> n = Rcpp::as<std::vector<double> >(xl["N"]);
  std::vector<double> t2 = Rcpp::as<std::vector<double> >(xl["T2df"]);
  // do something clever here
  return(R_NilValue);
')

希望有帮助。否则，列表总是打开的......

PS至于二维数组，这更棘手，因为没有原生 C++ 二维数组。如果你真的想做线性代数，看看RcppArmadillo和RcppEigen。

于 2011-11-16T23:53:19.370 回答