Use to∞ for TransformVariables

JuliaOcean · Nov 12, 2020 · a7551ac · a7551ac · briochemc · Nov 12, 2020
1 parent 07af753
commit a7551ac
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 4 deletions.
diff --git a/Project.toml b/Project.toml
@@ -1,7 +1,7 @@
 name = "AIBECS"
 uuid = "ace601d6-714c-11e9-04e5-89b7fad23838"
 authors = ["Benoit Pasquier <briochemc@gmail.com>"]
-version = "0.7.13"
+version = "0.7.14"
 
 [deps]
 BSON = "fbb218c0-5317-5bc6-957e-2ee96dd4b1f0"

diff --git a/src/Parameters.jl b/src/Parameters.jl
@@ -422,10 +422,9 @@ subfun(::Type{T}, s::Symbol) where {T<:AbstractParameters} = subfun(prior(T,s))
 ∇²subfun(::Type{T}, s::Symbol) where {T<:AbstractParameters} = ∇²subfun(prior(T,s))
 invsubfun(::Type{T}, s::Symbol) where {T<:AbstractParameters} = invsubfun(prior(T,s))
 # using TransformVariables
-lb_domain(d::ContinuousUnivariateDistribution) = (lb = support(d).lb; lb == -Inf ? TransformVariables.Infinity{false}() : lb)
-ub_domain(d::ContinuousUnivariateDistribution) = (ub = support(d).ub; ub == Inf ? TransformVariables.Infinity{true}() : ub)
 import TransformVariables: transform
-transform(d::ContinuousUnivariateDistribution) = as(Real, lb_domain(d), ub_domain(d))
+to∞(x) = isinf(x) ? TransformVariables.Infinity{x>0}() : x # convert Inf to ∞
+transform(d::ContinuousUnivariateDistribution) = as(Real, to∞(minimum(d)), to∞(maximum(d)))
 subfun(d::ContinuousUnivariateDistribution) = x -> transform(d)(x)
 ∇subfun(d::ContinuousUnivariateDistribution) = x -> ForwardDiff.derivative(subfun(d), x)
 ∇²subfun(d::ContinuousUnivariateDistribution) = x -> ForwardDiff.derivative(∇subfun(d), x)