Skip to content

Unify transition also in external samplers #2030

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
merged 12 commits into from
Jul 10, 2023
1 change: 0 additions & 1 deletion Project.toml
Original file line number Diff line number Diff line change
Expand Up @@ -39,7 +39,6 @@ Tracker = "9f7883ad-71c0-57eb-9f7f-b5c9e6d3789c"
[compat]
AbstractMCMC = "4"
AdvancedHMC = "0.3.0, 0.4"
AdvancedMH = "0.6.8, 0.7"
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Just noticed now, seems this was removed by accident?

AdvancedPS = "0.4"
AdvancedVI = "0.2"
BangBang = "0.3"
Expand Down
17 changes: 1 addition & 16 deletions src/contrib/inference/abstractmcmc.jl
Original file line number Diff line number Diff line change
Expand Up @@ -3,28 +3,13 @@ struct TuringState{S,F}
logdensity::F
end

struct TuringTransition{T,NT<:NamedTuple,F<:AbstractFloat}
θ::T
lp::F
stat::NT
end

function TuringTransition(vi::AbstractVarInfo, t)
theta = tonamedtuple(vi)
lp = getlogp(vi)
return TuringTransition(theta, lp, getstats(t))
end

metadata(t::TuringTransition) = merge((lp = t.lp,), t.stat)
DynamicPPL.getlogp(t::TuringTransition) = t.lp

state_to_turing(f::DynamicPPL.LogDensityFunction, state) = TuringState(state, f)
function transition_to_turing(f::DynamicPPL.LogDensityFunction, transition)
θ = getparams(transition)
varinfo = DynamicPPL.unflatten(f.varinfo, θ)
# TODO: `deepcopy` is overkill; make more efficient.
varinfo = DynamicPPL.invlink!!(deepcopy(varinfo), f.model)
return TuringTransition(varinfo, transition)
return Transition(varinfo, transition)
end

# NOTE: Only thing that depends on the underlying sampler.
Expand Down
11 changes: 6 additions & 5 deletions src/inference/Inference.jl
Original file line number Diff line number Diff line change
Expand Up @@ -123,6 +123,9 @@ end
######################
# Default Transition #
######################
# Default
# Extended in contrib/inference/abstractmcmc.jl
getstats(t) = nothing

struct Transition{T, F<:AbstractFloat, S<:Union{NamedTuple, Nothing}}
θ :: T
Expand All @@ -132,10 +135,10 @@ end

Transition(θ, lp) = Transition(θ, lp, nothing)

function Transition(vi::AbstractVarInfo; nt::NamedTuple=NamedTuple())
function Transition(vi::AbstractVarInfo, t=nothing; nt::NamedTuple=NamedTuple())
θ = merge(tonamedtuple(vi), nt)
lp = getlogp(vi)
return Transition(θ, lp, nothing)
return Transition(θ, lp, getstats(t))
end

function metadata(t::Transition)
Expand Down Expand Up @@ -664,9 +667,7 @@ function transitions_from_chain(
model(rng, vi, sampler)

# Convert `VarInfo` into `NamedTuple` and save.
theta = DynamicPPL.tonamedtuple(vi)
lp = Turing.getlogp(vi)
Transition(theta, lp)
Transition(vi)
end

return transitions
Expand Down