Method propagate

gecode/int/distinct/bnd.hpp:380–457

  template<class View>
  ExecStatus
  Bnd<View>::propagate(Space& home, const ModEventDelta& med) {
    assert(x.size() > 1);

    if (View::me(med) == ME_INT_VAL) {
      ExecStatus es = prop_val<View,false>(home,y);
      GECODE_ES_CHECK(es);
      if (y.size() < 2)
        return home.ES_SUBSUMED(*this);
      if (es == ES_FIX)
        return home.ES_FIX_PARTIAL(*this,View::med(ME_INT_BND));
    }

    if (y.size() == 2)
      GECODE_REWRITE(*this,(Rel::Nq<View,View>::post(home(*this),y[0],y[1])));

    ExecStatus es = prop_bnd<View>(home,x,min_x,max_x);

    GECODE_ES_CHECK(es);
    if (es == ES_NOFIX && View::me(modeventdelta()) ==  ME_INT_VAL)
      return es;

    const int n = x.size();

    if ((n > 2*y.size()) && (n > 6)) {
      // If there are many assigned views, try to eliminate them
      unsigned int d = static_cast<unsigned int>(max_x - min_x) + 1;
      if (d > 2*static_cast<unsigned int>(n)) {
        MinInc<View> min_inc;
        Support::quicksort<View,MinInc<View> >(&x[0], n, min_inc);
      } else {
        Region r;
        int* minbucket = r.alloc<int>(d);
        View* minsorted = r.alloc<View>(n);

        for (unsigned int i=0; i<d; i++)
          minbucket[i]=0;
        for (int i=0; i<n; i++)
          minbucket[x[i].min() - min_x]++;

        int c_min = 0;
        for (unsigned int i=0; i<d; i++) {
          int t_min = minbucket[i];
          minbucket[i] = c_min; c_min += t_min;
        }
        assert(c_min == n);

        for (int i=0; i<n; i++)
          minsorted[minbucket[x[i].min() - min_x]++] = x[i];
        for (int i=0; i<n; i++)
          x[i] = minsorted[i];
      }

      int i   = 0;
      int j   = 0;
      int max = x[0].max()-1;
      while (i < n-1) {
        if (!x[i].assigned() ||