I want to use the black-box likelihood function with NUTS sampler. Is there a way of using libraries like Autograd for automatic differentiation to get the gradients that the NUTS sampler requires.
In summary, I want to know if there is a feasible implementation of automatic differentiation that enables the NUTS sampler to work with black-box likelihoods.
Tensorflow Probability implements such an automatic differentiation mechanism when using black-box likelihoods.