Hierarchical Gaussian Process

No, no specific support for multiple outputs or a hierarchical GP has been added (though it would be great if it was there!). Does your code here run fine?