我实现了一个算法,在 C++
返回一个巨大的元素数组作为输出。现在,我想在这个函数中实现一个包装器 Rcpp
这样,我就可以通过使用 R
.
我在Makevars文件中指定了以下设置。
PKG_CXXFLAGS = -std=c++11
以便我可以使用C++11版本。
// [[Rcpp::export]]
NumericMatrix compute(int width, int height)
{
vector<data_t> weights(width * height);
compute_weights(weights);
NumericMatrix mat(height, width);
copy(begin(weights), end(weights), mat.begin());
return mat;
}
如果NumericMatrix在函数返回时被移动,则上述封装函数仍然有效,否则将创建一个新对象。
是否 Rcpp
利用移动语义?如果没有,有没有什么变通的办法来避免复制的构造?
如果NumericMatrix在函数返回时被移动,则上述封装函数仍然有效,否则将创建一个新对象。
... 如果没有,有没有什么变通的办法来避免拷贝的构造?
我认为复制构造函数只创建了一个浅层的副本,所以不应该有任何副本。请看 Rcpp: 如何确保NumericMatrix的深度复制? 和
这个例子也证实了这一点
#include <Rcpp.h>
// [[Rcpp::export]]
Rcpp::NumericVector allocate_the_vec(R_xlen_t n_ele){
Rcpp::NumericVector out(n_ele);
return out;
}
/*** R
# got 16 GB ram on my laptop. 3 x 7 is an issue but 2 x 7 is not
how_large <- as.integer(7 * 10^9 / 8)
the_large_vec_1 <- allocate_the_vec(how_large)
object.size(the_large_vec_1)
the_large_vec_2 <- allocate_the_vec(how_large)
*/