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