Vanilla implementation of Gaussian Process in PyMC3?

yup, it looks like its not doing it right. I left a comment at the spot.