Using pm.DensityDist and customized likelihood with a black-box function

Questions
9

Did you validate the gradient: http://docs.pymc.io/advanced_theano.html (see near the end)?
Also, you might find recent PR from @mattpitkin helpful: https://github.com/pymc-devs/pymc3/pull/3169