From dd27b285d6376b56f77d5969dfb13961e78d5444 Mon Sep 17 00:00:00 2001
From: James Hensman <james.hensman@gmail.com>
Date: Mon, 15 Apr 2013 17:53:26 +0100
Subject: [PATCH] rbf computation of psi2 now works if there's only one datum

---
 GPy/kern/rbf.py | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/GPy/kern/rbf.py b/GPy/kern/rbf.py
index a654cd0f..9ff7a93e 100644
--- a/GPy/kern/rbf.py
+++ b/GPy/kern/rbf.py
@@ -239,13 +239,13 @@ class rbf(kernpart):
         psi2 = np.empty((N,M,M))
 
         psi2_Zdist_sq = self._psi2_Zdist_sq
-        half_log_psi2_denom = 0.5*np.log(self._psi2_denom).squeeze()
+        _psi2_denom = self._psi2_denom.squeeze().reshape(N,self.D)
+        half_log_psi2_denom = 0.5*np.log(self._psi2_denom).squeeze().reshape(N,self.D)
         variance_sq = float(np.square(self.variance))
         if self.ARD:
             lengthscale2 = self.lengthscale2
         else:
             lengthscale2 = np.ones(Q)*self.lengthscale2
-        _psi2_denom = self._psi2_denom.squeeze()
         code = """
         double tmp;