如何将Rcpp NumericVector映射到静态/固定大小的Eigen向量？

要实现将Rcpp::NumericVector转换为Eigen固定大小容器（比如Eigen::Matrix<double,4,1>），你需要注意两点：一是手动构造Eigen::Map（因为RcppEigen对静态大小Map的自动转换支持有限），二是必须确保输入向量的长度和固定容器的大小一致，避免内存越界。

正确实现代码 #

方法1：构造固定大小的Eigen::Map

直接通过NumericVector的底层指针构造Map，同时添加长度检查：

Rcpp::cppFunction("bool conv(NumericVector X) { 
    // 先检查输入长度是否匹配固定大小
    if (X.size() != 4) {
        Rcpp::stop(\"输入向量长度必须为4！\");
    }
    // 用X的底层数据指针构造固定大小的Map
    Eigen::Map<Eigen::Matrix<double,4,1>> XS(X.begin());
    // 这里可以添加对XS的操作，比如打印验证
    // Rcpp::Rcout << XS << std::endl;
    return true; 
}", depends="RcppEigen")
conv(1:4)

方法2：转换为固定大小的Eigen容器（而非Map）

如果需要将数据拷贝到真正的固定大小Eigen容器（不是映射），可以在构造Map后进行赋值：

Rcpp::cppFunction("bool conv(NumericVector X) { 
    if (X.size() != 4) {
        Rcpp::stop(\"输入向量长度必须为4！\");
    }
    Eigen::Matrix<double,4,1> XS = Eigen::Map<Eigen::Matrix<double,4,1>>(X.begin());
    // 此时XS是独立的容器，可修改操作
    // XS(0) = 100;
    return true; 
}", depends="RcppEigen")
conv(1:4)

原代码问题说明 #

你之前尝试用Rcpp::as<Eigen::Map<Eigen::Matrix<double,4,1>>>(X)失败，是因为RcppEigen的自动转换模板主要针对动态大小的Eigen类型（比如VectorXd），对于静态大小的Map没有预定义的转换规则，因此需要手动通过指针构造。

内容的提问来源于stack exchange，提问作者Taylor