Definire una matrice in R e passarla a C++

Question

Ho una matrice definita in R. Ho bisogno di passare questa matrice a una funzione C++ e fare operazioni in C++. Esempio: In R, definire una matrice,

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

devo passare questa matrice a una funzione C++ dove verrà inizializzata variabile 'dati' di tipo vector<vector<double>> con la matrice A.

couldn' t sembra capire come fare questo. Sto pensando in un modo più complicato, quindi dovrei esserlo, scommetto che c'è un modo semplice per farlo.

Answer 1

Come ha detto Paul, mi consiglia di utilizzare Rcpp per questo tipo di cose. Ma dipende anche da cosa vuoi che il tuo vector< vector<double> > significhi. Supponendo che si desidera memorizzare le colonne, si potrebbe elaborare la matrice in questo modo:

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  

Answer 2

Probabilmente si desidera utilizzare Rcpp. Questo pacchetto consente una facile integrazione di R e C++, incluso il passaggio di oggetti da R a C++. Il pacchetto è disponibile su CRAN. Inoltre, un certo numero di pacchetti su CRAN utilizzano Rcpp, in modo che possano fungere da ispirazione. Il sito web di Rcpp è qui:

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

che include un qualche tutorial.