diff --git a/GPy/core/gp_grid.py b/GPy/core/gp_grid.py
index 89b6cdec..bdb1b614 100644
--- a/GPy/core/gp_grid.py
+++ b/GPy/core/gp_grid.py
@@ -73,7 +73,7 @@ class GpGrid(GP):
                 G[d] = len(A[d])
             N = np.prod(G)
             for d in range(D-1, -1, -1):
-                X = np.reshape(x, (G[d], round(N/G[d])), order='F')
+                X = np.reshape(x, (G[d], np.round(N/G[d])), order='F')
                 Z = np.dot(A[d], X)
                 Z = Z.T
                 x = np.reshape(Z, (-1, 1), order='F')
diff --git a/GPy/inference/latent_function_inference/gaussian_grid_inference.py b/GPy/inference/latent_function_inference/gaussian_grid_inference.py
index 04337e04..abdd78a3 100644
--- a/GPy/inference/latent_function_inference/gaussian_grid_inference.py
+++ b/GPy/inference/latent_function_inference/gaussian_grid_inference.py
@@ -41,7 +41,7 @@ class GaussianGridInference(LatentFunctionInference):
             G[d] = len(A[d])
         N = np.prod(G)
         for d in range(D-1, -1, -1):
-            X = np.reshape(x, (G[d], round(N/G[d])), order='F')
+            X = np.reshape(x, (G[d], np.round(N/G[d])), order='F')
             Z = np.dot(A[d], X)
             Z = Z.T
             x = np.reshape(Z, (-1, 1), order='F')