c++ - 从 NumericVector 填充 C++ 数组的更好习惯用法

Question

我正在使用 Rcpp 来包装一个用类似 C 的 C++ 编写(不是我)的算法(据我所知，没有 STL，没有提升，什么都没有)。可以看到实现的算法here (我正在包装 kmeans_w_03)。因此，我从 R 中传入一个 numeric vector ，然后需要将其转换为一个 double 数组。

目前我正在逐个元素地循环并从'tother'填充一个，像这样:

SEXP testfn(SEXP weightvec, SEXP cluster_num_k){
    Rcpp::NumericVector weightR(weightvec) ;
    int point_num = weightR.size();
    double weight[point_num] ;
    for(int i = 0; i < point_num; ++i) {
      weight[i] = weightR[i];
    }
}

但是对于单元素数值 vector ，我可以利用 Rcpp 的漂亮的 as 转换功能:

int cluster_num = Rcpp::as<int>(cluster_num_k);

然而，对长度 > 1 的数字 vector 尝试类似的操作会导致崩溃或错误，具体取决于语法的确切变体:

double weight[point_num] = Rcpp::as<double>(weightvec);

我不一定介意循环，但我完全是个新手，怀疑有更好的方法。我已经通读了 Rcpp 介绍、hadley 的维基教程和 RcppExamples，但还没有找到任何解决这个问题的方法，但这并不意味着我没有错过它。我对 Doxygen Rcpp 文档的阅读是，as 可以转换为 STL vector ，但不能转换为数组(但我很难阅读这些文档，所以我怀疑我错了)。如果是这样，我想我可以转换为一个 vector ，然后转换为一个数组....

所以我的问题是:是否有更好的(更少的代码行、更具表现力的代码，甚至可能避免内存分配)方法将 NumericVector 转换为双[]?

Answer 1

Ari，在之前的回答中，John 提出了一个很好的观点。延伸一点，这里有几个问题

• 您正在处理一个具有过时编码标准的库，它需要一个 double[]或 *double
• 您想使用更好的编码标准并利用 Rcpp .

好吧，不用担心，我们有一个简单的解决方案。实例化 Rcpp::NumericVector X(weightvec);像往常一样，然后将其传递给您的函数 foo() (或其他)作为

foo(X.begin())

它为您提供所需的 double* , 如果需要 X.size()提供长度。因为你来自 R，所以你不需要担心范围和生命周期，因为你之后返回到 R。

如果您需要更明确，我还使用了更丑陋的 &(X[0]) .

另外，Rcpp::as<>() caster 也适用于 std::vector<double> , 所以你可以做

std::vector<double> x = Rcpp::as<std::vector<double> >(weightvec);

但除了转换为 C++ 类型之外，它在这里没有任何收获。

最后，如果 Josh 或其他 C 顽固分子在附近，所有这一切都使用了 SEXP 的事实来自 R 的你保证有连续的 C 指针，所以你也可以通过非常古老的学校 REAL(weightvec) .