code: def in_domain(self, points):
          
          array_view = self.to_regular_array(points)
          non_negative = np.all(np.greater_equal(array_view, 0))
          correct_sum = np.all(np.sum(array_view, axis=-1) == self.n_meas)
          return non_negative and correct_sum
docs: Returns ``True`` if all of the given points are in the domain,
          ``False`` otherwise.
  
          :param np.ndarray points: An `np.ndarray` of type `self.dtype`.
  
          :rtype: `bool`