code: def _check_bounds(x, bounds):
      
      if bounds.shape == (0, 2):
          return True
      else:
          return ((np.atleast_2d(bounds)[:, 0] <= x).all() and
                  (x <= np.atleast_2d(bounds)[:, 1]).all())
docs: Checks whether `x` is within `bounds`. JIT-compiled in `nopython` mode
      using Numba.
  
      Parameters
      ----------
      x : ndarray(float, ndim=1)
          1-D array with shape (n,) of independent variables.
  
      bounds: ndarray(float, ndim=2)
          Sequence of (min, max) pairs for each element in x.
  
      Returns
      ----------
      bool
          `True` if `x` is within `bounds`, `False` otherwise.