RCPP：應用lfactorial到NumericVector

Question

什麼是對lfactorial功能適用於RCPP一個NumericVector的最好/最簡單的方法。在R I將有這樣的事情：

> vec <- c(1,2,3,4) 
> lfactorial(vec) 
[1] 0.0000000 0.6931472 1.7917595 3.1780538 

在RCPP，假設我有：

NumericVector vec = NumericVector::create(1,2,3,4); 

什麼是應用Rcpp::internal::lfactorial到vec每個元素的最好方法是什麼？

感謝您的任何幫助。

Answer 1

RCPP，就像R，被矢量：

R> cppFunction('NumericVector lfac(NumericVector x) { \ 
             return lfactorial(x); }') 
R> lfac(1:4) 
[1] 0.000000 0.693147 1.791759 3.178054 
R> 

也就是說代碼我只是打破了更好的闡述這裏一行。在複製和粘貼之前刪除反斜槓。你可以把它寫無論是作爲一個NumericVector或IntegerVector簽名。