Rcpp特徵映射與MatrixXf錯誤

Question

爲什麼下面的代碼不能編譯？

library(Rcpp) 
cppFunction(' 
    int rows(const NumericMatrix& X) { 
     using Eigen::MatrixXf; 
     typedef Eigen::Map<MatrixXf> MapMat; 
     MapMat X1(as<MapMat>(X)); 
     return X1.rows(); 
}', depends = "RcppEigen") 

它引發以下錯誤：當我用MatrixXd代替

error: no matching function for call to 'Eigen::Map<Eigen::Matrix<float, -1, -1> >::Map(Rcpp::Vector<14, Rcpp::PreserveStorage>::iterator, int&, int&)' 
     OUT get() {return OUT(vec.begin(), d_nrow, d_ncol);} 

相同的代碼工作正常。

謝謝。

Answer 1

NumericMatrix使用類型double（而不是float）。特徵不支持使用不同類型的矩陣之間的隱式類型轉換。您的代碼似乎嘗試讀取doubleNumericMatrix的內存作爲float特徵矩陣。只需使用MatrixXd類型。