This summer I focused primarily on improvements to the PyMC3 Gaussian process module. These items include support for multiple observations, additional capabilities for larger data sets, a "Latent" Kronecker implementation, and additional supporting documentation and in-depth examples.
Currently, the Latent Kronecker implementation is ready to merge pending a final review. All tests are passing and the example runs smoothly. The pull request is here, with the demo here.