Merge pull request #125 from devmotion/dw/log1pexp

Use `log1pexp` in `logistic_logjac`

Author: tpapp

src/utilities.jl

@@ -4,7 +4,7 @@
 
 function logistic_logjac(x::Real)
     mx = -abs(x)
-    mx - 2*log1p(exp(mx))
+    mx - 2*log1pexp(mx)
 end
 
 logit_logjac(y) = -log(y) - log1p(-y)