code: def contains(bounds, points):
      
      # make sure we have correct input types
      bounds = np.asanyarray(bounds, dtype=np.float64)
      points = np.asanyarray(points, dtype=np.float64)
  
      if len(bounds) != 2:
          raise ValueError('bounds must be (2,dimension)!')
      if not util.is_shape(points, (-1, bounds.shape[1])):
          raise ValueError('bounds shape must match points!')
  
      # run the simple check
      points_inside = np.logical_and(
          (points > bounds[0]).all(axis=1),
          (points < bounds[1]).all(axis=1))
  
      return points_inside
docs: Do an axis aligned bounding box check on a list of points.
  
      Parameters
      -----------
      bounds : (2, dimension) float
         Axis aligned bounding box
      points : (n, dimension) float
         Points in space
  
      Returns
      -----------
      points_inside : (n,) bool
        True if points are inside the AABB