5

我在 R 中定义了一个矩阵。我需要将此矩阵传递给 C++ 函数并在 C++ 中执行操作。示例：在 R 中，定义一个矩阵，

A <- matrix(c(9,3,1,6),2,2,byrow=T)
PROTECT( A = AS_NUMERIC(A) );
double* p_A = NUMERIC_POINTER(A);

我需要将此矩阵传递给 C++ 函数，其中类型的变量“数据”vector<vector<double>>将使用矩阵 A 初始化。

我似乎无法弄清楚如何做到这一点。我正在以比我应该的更复杂的方式思考，我敢打赌有一种简单的方法可以做到这一点。

4

2 回答 2

5

正如保罗所说，我建议使用Rcpp这种东西。但这也取决于你想要你vector< vector<double> >的意思。假设你想存储列，你可以像这样处理你的矩阵：

require(Rcpp)
require(inline)

fx <- cxxfunction( signature( x_ = "matrix" ), '
    NumericMatrix x(x_) ;
    int nr = x.nrow(), nc = x.ncol() ;
    std::vector< std::vector<double> > vec( nc ) ;
    for( int i=0; i<nc; i++){
        NumericMatrix::Column col = x(_,i) ;
        vec[i].assign( col.begin() , col.end() ) ;
    }
    // now do whatever with it
    // for show here is how Rcpp::wrap can wrap vector<vector<> >
    // back to R as a list of numeric vectors
    return wrap( vec ) ;
', plugin = "Rcpp" )
fx( A )
# [[1]]
# [1] 9 1
# 
# [[2]]
# [1] 3 6
于 2012-11-10T18:48:33.580 回答
5

您可能想使用 Rcpp。该包允许轻松集成 R 和 C++，包括将对象从 R 传递到 C++。该软件包可在 CRAN 上获得。此外，CRAN 上的一些包使用了 Rcpp，因此它们可以作为灵感。Rcpp的网站在这里：

http://dirk.eddelbuettel.com/code/rcpp.html

其中包括一些教程。

于 2012-11-10T18:31:16.190 回答