Определите матрицу в R и передайте ее переполнению стека

У меня есть матрица, определенная в 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>> будет инициализирован с матрицей А.

Я не мог понять, как это сделать. Я думаю, что более сложным, чем я должен быть, держу пари, что есть простой способ сделать это.

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
5

Другие решения

Вы, вероятно, хотите использовать Rcpp. Этот пакет позволяет легко интегрировать R и C ++, включая передачу объектов из R в C ++. Пакет доступен на CRAN. Кроме того, ряд пакетов в CRAN использует Rcpp, поэтому они могут служить источником вдохновения. Веб-сайт Rcpp находится здесь:

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

который включает в себя несколько уроков.

5

По вопросам рекламы [email protected]