code: def acquisition_function(self,x):
          
          f_acqu = self._compute_acq(x)
          cost_x, _ = self.cost_withGradients(x)
          return -(f_acqu*self.space.indicator_constraints(x))/cost_x
docs: Takes an acquisition and weights it so the domain and cost are taken into account.