If you follow the step in Frequently Asked Questions - #11 by junpenglao and there is no problem of the model test point, the derivative zero error is usually when you have prior being too flat - changing them to more informative prior usually works
There is a logsumexp function in pymc3.