scikit-learn
diff --git a/‎sklearn/gaussian_process/kernels.py
Copy file name to clipboardExpand all lines: sklearn/gaussian_process/kernels.py
+3-3Lines changed: 3 additions & 3 deletions b/‎sklearn/gaussian_process/kernels.py
Copy file name to clipboardExpand all lines: sklearn/gaussian_process/kernels.py
+3-3Lines changed: 3 additions & 3 deletions
@@ -1686,9 +1686,9 @@ def __call__(self, X, Y=None, eval_gradient=False):
                 D = squareform(dists**2)[:, :, np.newaxis]
 
             if self.nu == 0.5:
-                K_gradient = K[..., np.newaxis] * D \
-                    / np.sqrt(D.sum(2))[:, :, np.newaxis]
-                K_gradient[~np.isfinite(K_gradient)] = 0
+                denominator = np.sqrt(D.sum(axis=2))[:, :, np.newaxis]
+                K_gradient = K[..., np.newaxis] * \
+                    np.divide(D, denominator, where=denominator != 0)
             elif self.nu == 1.5:
                 K_gradient = \
                     3 * D * np.exp(-np.sqrt(3 * D.sum(-1)))[..., np.newaxis]