Auto-diff when using NUTS sampler with blackbox likelihoods

All of pymc3 is using automatic differentiation. :slight_smile:
We are using theano, so you can use arbitrary theano code in your model. If you wanted to use autograd, you would have to link theano and autograd together. Quite possible, but some work and knowledge of theano and autograd would be required for that.