Vanilla implementation of Gaussian Process in PyMC3?

Could you glance at the code in the new rethinking PR I linked and confirm it’s doing the wrong thing?