exact inference for N>D of Y

GPy/inference/latent_function_inference/exact_gaussian_inference.py

@@ -29,10 +29,12 @@ class ExactGaussianInference(object):
         """
         N, D = Y.shape
         if (N>D):
+            print "WARNING: D>N we still need caching of L, such that L*L^T = Y, although fine here"
             return Y
         else:
             #if Y in self.cache, return self.Cache[Y], else store Y in cache and return L.
-            raise NotImplementedError, 'TODO' #TODO
+            print "WARNING: N>D of Y, we need caching of L, such that L*L^T = Y, returning Y still!"
+            return Y
 
     def inference(self, kern, X, likelihood, Y, Y_metadata=None):
         """