Gen: General-purpose probabilistic programming model with programmable inference
dl.acm.org
dl.acm.org
And then the dynamic compiler mode uses AD to work with Julia packages, so you can easily set it up to do Bayesian inference on neural (Flux.jl) partial differential equations with high order adaptive implicit integrators with DifferentialEquations.jl, which is why I am starting to use this library :). At the same time, there is also a static compiled sublanguage of Gen which is extremely fast (but cannot use arbitrary Julia packages of course).
That's because JAGS only supports models that can be expressed as graphical models.