r - Rcpp 中 C++ 中的 .nrow() 方法返回零

Question

借助 Rcpp，我在 C++ 中定义了一个矩阵 M。使用 M.nrow()，我们应该能够检索行数。但是，当我尝试将行数作为 IntegerVector 返回时，答案不正确:

set.seed(1100)
M = matrix(sample(18), nrow = 6)
Rcpp::cppFunction('IntegerVector tccp5(IntegerMatrix M) { int x = M.nrow(); return x;}')
tccp5(M)
# [1] 0 0 0 0 0 0

正确答案应该是行数，例如

# [1] 6

你能解释一下发生了什么吗？

Answer 1

虽然@gfgm 的回答是正确的，但我想扩展为什么声明为 int 而不是 IntegerVector 会导致正确的构造...

特别是，每个返回 R 的 return 都使用 wrap() 无缝转换为 SEXP ，或指向该结构的 S 表达式。通过提供 IntegerVector 并返回一个 int，wrap() 必须实例化一个长度为 的新 IntegerVector x 因为 int 下没有 SEXP。另一方面，当定义的返回类型为 int 时，Rcpp 的 wrap() 特性能够正确地强制转换 int 转换为 IntegerVector。

---

为了强调发生的底层“无缝”转换，让我们将参数 verbose = TRUE 添加到 cppFunction() 以查看 C++ 代码是:编译、链接并导入到 R 中。 (注意:我已经截断了编译的输出。)

如果我们考虑:

Rcpp::cppFunction('IntegerVector tccp5(IntegerMatrix M) { int x = M.nrow(); return x;}',
                  verbose = TRUE)

我们得到:

Generated code for function definition: 
--------------------------------------------------------

#include <Rcpp.h>

using namespace Rcpp;

// [[Rcpp::export]]
IntegerVector tccp5(IntegerMatrix M) { int x = M.nrow(); return x;}

Generated extern "C" functions 
--------------------------------------------------------


#include <Rcpp.h>
// tccp5
IntegerVector tccp5(IntegerMatrix M);
RcppExport SEXP sourceCpp_7_tccp5(SEXP MSEXP) {
BEGIN_RCPP
    Rcpp::RObject rcpp_result_gen;
    Rcpp::RNGScope rcpp_rngScope_gen;
    Rcpp::traits::input_parameter< IntegerMatrix >::type M(MSEXP);
    rcpp_result_gen = Rcpp::wrap(tccp5(M));
    return rcpp_result_gen;
END_RCPP
}

相比于:

Rcpp::cppFunction('int tccp5(IntegerMatrix M) { int x = M.nrow(); return x;}', 
                  verbose = TRUE)

给出:

Generated code for function definition: 
--------------------------------------------------------

#include <Rcpp.h>

using namespace Rcpp;

// [[Rcpp::export]]
int tccp5(IntegerMatrix M) { int x = M.nrow(); return x;}

Generated extern "C" functions 
--------------------------------------------------------


#include <Rcpp.h>
// tccp5
int tccp5(IntegerMatrix M);
RcppExport SEXP sourceCpp_9_tccp5(SEXP MSEXP) {
BEGIN_RCPP
    Rcpp::RObject rcpp_result_gen;
    Rcpp::RNGScope rcpp_rngScope_gen;
    Rcpp::traits::input_parameter< IntegerMatrix >::type M(MSEXP);
    rcpp_result_gen = Rcpp::wrap(tccp5(M));
    return rcpp_result_gen;
END_RCPP
}