WIP, successor to Soss.jl
For a high-level description of Tilde's design, check out this blog post
julia> m = @model (x, s) begin
σ ~ Exponential()
@inbounds x[1] ~ Normal(σ = σ)
n = length(x)
@inbounds for j = 2:n
x[j] ~ StudentT(1.5, x[j - 1], σ)
end
end;
julia> x = zeros(3);
julia> r = rand(m(x,10))
(σ = 9.096155145583953, x = [14.876934886768867, 6.612967845518229, 2.045770246490428])
julia> x
3-element Vector{Float64}:
14.876934886768867
6.612967845518229
2.045770246490428
julia> ℓ = logdensityof(m(x, 1.0) | (;x), (σ = 1.0,))
-122.91114458882001