vendor/tomotopy/src/Utils/math.h in tomoto-0.1.2 vs vendor/tomotopy/src/Utils/math.h in tomoto-0.1.3

- old
+ new

@@ -152,10 +152,10 @@ inline float lgammaT(float x) { return detail::LUT_lgamma::get(x); } inline float digammaT(float x) { return detail::LUT_digamma::get(x); } template<class _T> - inline _T lgammaApprox(_T z) + inline auto lgammaApprox(_T z) -> decltype((z + 2.5)* log(z + 3) - (z + 3) + 0.91893853 + 1. / 12. / (z + 3) - log(z * (z + 1) * (z + 2))) { // approximation : lgamma(z) ~= (z+2.5)ln(z+3) - z - 3 + 0.5 ln (2pi) + 1/12/(z + 3) - ln (z(z+1)(z+2)) return (z + 2.5) * log(z + 3) - (z + 3) + 0.91893853 + 1. / 12. / (z + 3) - log(z * (z + 1) * (z + 2)); } \ No newline at end of file