c++ - 在 Rcpp 中分配一个大型二维数组 NumericVectors

Question

我正在尝试分配一个大型(大概)2D 数组的 NumericVectors(总内存应该只有 16Mb 左右，除非有开销)，但我遇到了堆栈溢出。这是我能想到的最小的可重现示例:

require(Rcpp)
require(inline)
testfun = cxxfunction(
    signature(x="list"),
    body='using namespace std;
        vector<vector<NumericVector> > rs=as< vector<vector<NumericVector> > >( x );
        return wrap(rs);',
    plugin="Rcpp")
x=lapply(1:1000,function(g) lapply(1:1000, function(h) runif(3)))
testfun(x)

使用 gdb 进行分析只告诉我，我在 libR.so 递归中非常深入。请注意，虽然这里的数组是矩形的，所以我可以在我的实际问题中使用 NumericMatrix，但它们会呈锯齿状。感谢您的任何想法。

Answer 1

编辑:使用(相当新的)ListOf<T> 的解决方案模板类:

require(Rcpp)
require(inline)
testfun = cxxfunction(
  signature(x="list"),
  body='using namespace std;
        ListOf<ListOf<NumericVector> > rs=as< ListOf<ListOf<NumericVector> > >( x );
        return wrap(rs);',
  plugin="Rcpp")
x=lapply(1:1000,function(g) lapply(1:1000, function(h) runif(3)))
testfun(x)

我们只是替换所有对 std::vector 的调用与 Rcpp::ListOf而且，有点神奇的是，它起作用了。最重要的是 - 数据不会复制到此处(假设 x 确实是一个数字 vector 列表)。