code: def lookup(self, coordinates):
          
          # Convert coordinates to indices
          index_array = np.clip(np.searchsorted(self.bin_edges, coordinates) - 1,
                                0,
                                len(self.bin_edges) - 2)
  
          # Use the index array to slice the histogram
          return self.histogram[index_array]
docs: Lookup values at coordinates.
          coordinates: arraylike of coordinates.
  
          Clips if out of range!! TODO: option to throw exception instead.
          TODO: Needs tests!!