code: def propose_unif(self):
          
  
          while True:
              # Sample a point from the union of ellipsoids.
              # Returns the point `u`, ellipsoid index `idx`, and number of
              # overlapping ellipsoids `q` at position `u`.
              u, idx, q = self.mell.sample(rstate=self.rstate, return_q=True)
  
              # Check if the point is within the unit cube.
              if unitcheck(u, self.nonperiodic):
                  # Accept the point with probability 1/q to account for
                  # overlapping ellipsoids.
                  if q == 1 or self.rstate.rand() < 1.0 / q:
                      break  # if successful, we're done!
  
          return u, self.mell.ells[idx].axes
docs: Propose a new live point by sampling *uniformly* within
          the union of ellipsoids.