If anyone is so inclined, there is a fair gaussian process method that is written in Matlab that would be great to port to Python. https://github.com/ztanml/fgp