我可以将Eigen::VectorXi转换为Eigen::VectorXd吗？

Question

我在R上使用RcppEigen，我想从R层取一个双向量：

RcppExport SEXP testfunc (SEXP a) {
    const Eigen::Map<Eigen::VectorXd> vecD(as<Eigen::Map<Eigen::VectorXd> >(a));

但是，number是整型的，所以我需要像这样获取整型：

RcppExport SEXP testfunc (SEXP a) {
    const Eigen::Map<Eigen::VectorXi> vecD(as<Eigen::Map<Eigen::VectorXi> >(a));

所以，我需要把它转换成双倍数。我可以将Eigen::VectorXi转换为Eigen::VectorXd吗？

Answer 1

当你指的是numeric，也就是double，或者反之亦然时，你确定你没有混淆从R调用integer吗？没有理由不使用两个函数，或者在函数内部进行调度。

例如，考虑这里的代码(使用Rcpp 0.10.0特性)

#include <RcppEigen.h>

// [[Rcpp::depends(RcppEigen)]]

// [[Rcpp::export]]
double vecdSum(SEXP x) {
    const Eigen::Map<Eigen::VectorXd> 
       vec(Rcpp::as<Eigen::Map<Eigen::VectorXd> >(x));
    return vec.sum();
}

// [[Rcpp::export]]
int veciSum(SEXP x) {
    const Eigen::Map<Eigen::VectorXi> 
       vec(Rcpp::as<Eigen::Map<Eigen::VectorXi> >(x));
    return vec.sum();
}

我们可以很容易地通过

R> sourceCpp('/tmp/vecsums.cpp')    # Rcpp 0.10.0 adds this
R> veciSum(c(1L, 2L, 3L))
[1] 6
R> vecdSum(c(1, 2, 3))
[1] 6
R>