rcpp：如何将gamma函数应用于标量？

Question

我想函数gamma only works for a vector as the input. Is there a way to apply it to a scalar, say, gamma（3）`？

实际上，如果我将gamma(3)作为我的代码的一部分，我会得到正确的输出，但是有一条警告信息......

#include <Rcpp.h>
using namespace Rcpp;

// [[Rcpp::export]]
List fool(NumericVector vec){
  double res = 0;
  res = sum(gamma(vec)) + gamma(3);
  List result;result["num"] = res;
   return result;
}

这是警告乱码：

  

exp.cpp:7:27：警告：＆＃39; gamma＆＃39;不推荐使用：首先在OS X 10.9 [-Wdeprecated-declarations]中弃用     res = sum（gamma（vec））+ gamma（3）;
                            ^   /usr/include/math.h:720:15：注意：＆＃39; gamma＆＃39;已明确标记为已弃用   extern double gamma（double）__ OSX_AVAILABLE_BUT_DEPRECATED（__ MAC_10_0，__ MAC_10_9，__ IPHONE_NA，__ IPHONE_NA）;                 ^   生成了1个警告。

Answer 1

感谢您发布代码。你成了名字空间粗心的受害者。在Rcpp名称空间中有一个（向量化的）gamma() - 第一个参数，并且在C数学库中有一个标量gamma()（首选：tgamma()）。最好是明确的。

更正了以下代码：

#include <Rcpp.h>

// [[Rcpp::export]]
double fool(Rcpp::NumericVector vec){
    double res =
        Rcpp::sum(Rcpp::gamma(vec)) +   // Rcpp sugar sum() and gamma() on vector
        ::tgamma(3.0);                  // math library tgamma of double
    return res;
}

/*** R
v <- 1:5
fool(v)
sum(gamma(v)) + gamma(3)
*/

输出

R> sourceCpp("/tmp/shuang.cpp")

R> v <- 1:5

R> fool(v)
[1] 36

R> sum(gamma(v)) + gamma(3)
[1] 36
R>