Vanilla implementation of Gaussian Process in PyMC3?

Thanks, very helpful!