1 #ifndef STAN_MATH_PRIM_SCAL_FUN_BINOMIAL_COEFFICIENT_LOG_HPP 2 #define STAN_MATH_PRIM_SCAL_FUN_BINOMIAL_COEFFICIENT_LOG_HPP 8 #include <boost/math/tools/promotion.hpp> 62 template <
typename T_N,
typename T_n>
63 inline typename boost::math::tools::promote_args<T_N, T_n>::type
66 const double CUTOFF = 1000;
68 const T_N N_plus_1 = N + 1;
71 typename boost::math::tools::promote_args<T_N, T_n>::type N_minus_n = N - n;
72 const double one_twelfth =
inv(12);
74 + one_twelfth / N - n - one_twelfth / N_minus_n -
lgamma(n + 1);
fvar< T > binomial_coefficient_log(const fvar< T > &x1, const fvar< T > &x2)
fvar< T > lgamma(const fvar< T > &x)
Return the natural logarithm of the gamma function applied to the specified argument.
fvar< T > log(const fvar< T > &x)
fvar< T > multiply_log(const fvar< T > &x1, const fvar< T > &x2)
fvar< T > inv(const fvar< T > &x)