try stablelize kuma sample

xukai92 · Apr 28, 2019 · bc4823b · bc4823b
1 parent dc13268
commit bc4823b
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/src/distributions/beta.jl b/src/distributions/beta.jl
@@ -15,7 +15,8 @@ end
 function _u2logkumaraswamysample(T, u, kuma::BatchKumaraswamy)
     _one = one(T)
     _eps = eps(T)
-    return log.(_one .- exp.(log.(u .+ _eps) ./ kuma.b)) ./ kuma.a
+    # return log.(_one .- exp.(log.(u .+ _eps) ./ kuma.b)) ./ kuma.a
+    return log1p.(-exp.(log.(u .+ _eps) ./ (kuma.b .+ _eps))) ./ (kuma.a .+ _eps)
 end
 
 """