Using pymc3 NUTS sampler for an "external model" having its derivatives with respect to its parameters

Thank you. I am working on that and will update here.