将 Rcpp::NumericVector 转换为 Eigen::VectorXd

Question

之前有人问过类似的问题，这里：

Converting between NumericVector/Matrix and VectorXd/MatrixXd in Rcpp(Eigen) to perform Cholesky solve

问题是，我在 fastLm.cpp（最后）中找到的代码对我不起作用。

Rcpp::NumericVector X( (SEXP) R.parseEval("x <- 1:10"));

Eigen::Map<Eigen::VectorXd> XS(Rcpp::as<Eigen::Map<Eigen::VectorXd>>(X)); //!!not working

它给出了以下错误：

error: no matching constructor for initialization of 'Eigen::Map<Eigen::Matrix<double, -1, 1, 0, -1, 1>, 0, Eigen::Stride<0, 0> >'
                    Exporter( SEXP x ) : t(x){}

和（除其他外）：

note: candidate constructor (the implicit copy constructor) not viable: cannot convert argument of incomplete type 'SEXP' (aka 'SEXPREC *') to 'const Eigen::Map<Eigen::Matrix<double, -1, 1, 0, -1, 1>, 0, Eigen::Stride<0, 0> >'
template<typename PlainObjectType, int MapOptions, typename StrideType> class Map

可能与SEXP对象有关？

Answer 1

如果我在 >>:

之间再添加一个 space 对我有用

R> cppFunction("bool conv(NumericVector X) { \
      Eigen::Map<Eigen::VectorXd> \
      XS(Rcpp::as<Eigen::Map<Eigen::VectorXd> >(X)); return true; } ", 
           depends="RcppEigen")
R> conv(1:4)
[1] TRUE
R>

此外，我认为你不能只说

X((SEXP) R.parseEval("x <- 1:10"))

除非您使用 RInside，它有自己的 Eigen 示例——为此只需执行

cd examples/eigen
make
./rinside_sample0
./rinside_sample1

两者在我的安装中仍然运行。

将 Rcpp::NumericVector 转换为 Eigen::VectorXd

convert Rcpp::NumericVector to Eigen::VectorXd

c++

rcpp

eigen