Automatic Differentiation
Switching AD Modes
Turing supports four automatic differentiation (AD) packages in the back end during sampling. The default AD backend is ForwardDiff for forward-mode AD. Three reverse-mode AD backends are also supported, namely Tracker, Zygote and ReverseDiff. Zygote
and ReverseDiff
are supported optionally if explicitly loaded by the user with using Zygote
or using ReverseDiff
next to using Turing
.
To switch between the different AD backends, one can call the function Turing.setadbackend(backend_sym)
, where backend_sym
can be :forwarddiff
(ForwardDiff
), :tracker
(Tracker
), :zygote
(Zygote
) or :reversediff
(ReverseDiff.jl
). When using ReverseDiff
, to compile the tape only once and cache it for later use, the user has to call Turing.setrdcache(true)
. However, note that the use of caching in certain types of models can lead to incorrect results and/or errors.
Compiled tapes should only be used if you are absolutely certain that the computation doesn't change between different executions of your model.
Thus, e.g., in the model definition and all im- and explicitly called functions in the model all loops should be of fixed size, and if
-statements should consistently execute the same branches.
For instance, if
-statements with conditions that can be determined at compile time or conditions that depend only on the data will always execute the same branches during sampling (if the data is constant throughout sampling and, e.g., no mini-batching is used).
However, if
-statements that depend on the model parameters can take different branches during sampling; hence, the compiled tape might be incorrect.
Thus you must not use compiled tapes when your model makes decisions based on the model parameters, and you should be careful if you compute functions of parameters that those functions do not have branching which might cause them to execute different code for different values of the parameter.
Compositional Sampling with Differing AD Modes
Turing supports intermixed automatic differentiation methods for different variable spaces. The snippet below shows using ForwardDiff
to sample the mean (m
) parameter and using the Tracker-based TrackerAD
autodiff for the variance (s
) parameter:
using Turing
# Define a simple Normal model with unknown mean and variance.
@model function gdemo(x, y)
s² ~ InverseGamma(2, 3)
m ~ Normal(0, sqrt(s²))
x ~ Normal(m, sqrt(s²))
return y ~ Normal(m, sqrt(s²))
end
# Sample using Gibbs and varying autodiff backends.
c = sample(
gdemo(1.5, 2),
Gibbs(HMC{Turing.ForwardDiffAD{1}}(0.1, 5, :m), HMC{Turing.TrackerAD}(0.1, 5, :s²)),
1000,
)
Chains MCMC chain (1000×3×1 Array{Float64, 3}):
Iterations = 1:1:1000
Number of chains = 1
Samples per chain = 1000
Wall duration = 3.28 seconds
Compute duration = 3.28 seconds
parameters = s², m
internals = lp
Summary Statistics
parameters mean std mcse ess_bulk ess_tail rhat
e ⋯
Symbol Float64 Float64 Float64 Float64 Float64 Float64
⋯
s² 3.1808 6.0544 1.6290 20.6813 13.5415 1.0268
⋯
m 0.7506 1.5547 0.4181 41.1135 12.3437 1.0175
⋯
1 column om
itted
Quantiles
parameters 2.5% 25.0% 50.0% 75.0% 97.5%
Symbol Float64 Float64 Float64 Float64 Float64
s² 0.6041 1.0567 1.5934 2.5979 17.8584
m -4.4143 0.4816 1.0450 1.5431 2.5570
Generally, TrackerAD
is faster when sampling from variables of high dimensionality (greater than 20), and ForwardDiffAD
is more efficient for lower-dimension variables. This functionality allows those who are performance sensitive to fine-tune their automatic differentiation for their specific models.
If the differentiation method is not specified in this way, Turing will default to using whatever the global AD backend is. Currently, this defaults to ForwardDiff
.