Gibbs doesn't work with NUTS #2520
Description
Minimal working example
using Turing, Turing.DynamicPPL
@model function coinflips(y)
p ~ Beta(1, 1)
y .~ Bernoulli(p)
end
m = coinflips(rand(Bool, 100))
# This gets stuck!
sample(m, Gibbs(:p => NUTS()), 1000)
# This works fine
sample(m, NUTS(), 1000)
# This also works fine
sample(m, Gibbs(:p => HMC(0.01, 20)), 1000)Description
Gibbs(:p => NUTS()) just totally gets stuck and doesn't move at all. Am I doing something wrong here?
Julia version info
versioninfo()
Julia Version 1.11.3
Commit d63adeda50 (2025-01-21 19:42 UTC)
Build Info:
Official https://julialang.org/ release
Platform Info:
OS: Windows (x86_64-w64-mingw32)
CPU: 8 × 11th Gen Intel(R) Core(TM) i7-1185G7 @ 3.00GHz
WORD_SIZE: 64
LLVM: libLLVM-16.0.6 (ORCJIT, tigerlake)
Threads: 7 default, 0 interactive, 3 GC (on 8 virtual cores)
Environment:
JULIA_EDITOR = code
JULIA_NUM_THREADS = 7
Manifest
]st --manifest
(Paste here)