patch v1.2.2 fix domain error by switching logccdf function to old Rm…

…ath version
wzhorton · May 16, 2023 · 38ca7a0 · 38ca7a0 · wzhorton · May 16, 2023
1 parent deecd41
commit 38ca7a0
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/Project.toml b/Project.toml
@@ -1,7 +1,7 @@
 name = "PolyaGammaHybridSamplers"
 uuid = "c636ee4f-4591-4d8c-9fae-2dea21daa433"
 authors = ["W.Z. Horton"]
-version = "1.2.1"
+version = "1.2.2"
 
 [deps]
 Distributions = "31c24e10-a181-5473-b8eb-7969acd0382f"

diff --git a/src/rand_pgsaddlepoint.jl b/src/rand_pgsaddlepoint.jl
@@ -78,7 +78,7 @@ function rand_Jstar_tilted(b::Real, z::Real, rng::AbstractRNG)
          p_igauss(md, inv(rt2rl), b)
 
     wr = exp(0.5 * log(ar) + lcn + (-b * log(b * rr) + b * ir - b * log(md) + loggamma(b)) +
-         logccdf(Gamma(b, inv(b*rr)), md))
+        StatsFuns.RFunctions.gammalogccdf(b, inv(b*rr), md))
 
     wt = wl + wr
     pl = wl * inv(wt)