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

Could you please try to upgrade your pymc3 and theano?