@@ -264,19 +264,14 @@ class Learner1D(BaseLearner):

         missing_bounds = [b for b in self.bounds if b not in self.data

                           and b not in self.pending_points]

-        if len(missing_bounds) == 2:

-            # First time

+        if missing_bounds:

             loss_improvements = [np.inf] * n

-            points = np.linspace(*self.bounds, n).tolist()

-        elif len(missing_bounds) == 1:

-            loss_improvements = [np.inf] * n

-            points = np.linspace(*self.bounds, n + 1).tolist()

-            if missing_bounds[0] == self.bounds[1]:

-                # Second time, if we previously returned just self.bounds[0]

-                points = points[1:]

-            else:

-                # Rare case in which self.bounds[1] is present before self.bounds[1]

-                points = points[:1]

+            points = np.linspace(*self.bounds, n + 2 - len(missing_bounds)).tolist()

+            if len(missing_bounds) == 1:

+                # If we previously returned just self.bounds[0] we exclude that point.

+                # In the rare case in which self.bounds[1] is present before self.bounds[1]

+                # we exclude that point.

+                points = points[1:] if missing_bounds[0] == self.bounds[1] else points[:-1]

         else:

             def xs(x_left, x_right, n):

                 if n == 1:

@@ -311,7 +306,6 @@ class Learner1D(BaseLearner):

         return points, loss_improvements

-

     def plot(self):

         hv = ensure_holoviews()

         if not self.data: