re-subset.icc

Go to the documentation of this file.

/* -*- mode: C++; c-basic-offset: 2; indent-tabs-mode: nil -*- */
/*
 *  Main authors:
 *     Guido Tack <tack@gecode.org>
 *     Christian Schulte <schulte@gecode.org>
 *
 *  Copyright:
 *     Guido Tack, 2004
 *     Christian Schulte, 2004
 *
 *  Last modified:
 *     $Date: 2008-07-11 09:31:51 +0200 (Fri, 11 Jul 2008) $ by $Author: tack $
 *     $Revision: 7288 $
 *
 *  This file is part of Gecode, the generic constraint
 *  development environment:
 *     http://www.gecode.org
 *
 *  Permission is hereby granted, free of charge, to any person obtaining
 *  a copy of this software and associated documentation files (the
 *  "Software"), to deal in the Software without restriction, including
 *  without limitation the rights to use, copy, modify, merge, publish,
 *  distribute, sublicense, and/or sell copies of the Software, and to
 *  permit persons to whom the Software is furnished to do so, subject to
 *  the following conditions:
 *
 *  The above copyright notice and this permission notice shall be
 *  included in all copies or substantial portions of the Software.
 *
 *  THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
 *  EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
 *  MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
 *  NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE
 *  LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
 *  OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
 *  WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
 *
 */

namespace Gecode { namespace Set { namespace Rel {

  template <class View0, class View1>
  forceinline
  ReSubset<View0,View1>::ReSubset(Space* home, View0 y0,
                                  View1 y1, Gecode::Int::BoolView y2)
    : Propagator(home), x0(y0), x1(y1), b(y2) {
    b.subscribe(home,this, Gecode::Int::PC_INT_VAL);
    x0.subscribe(home,this, PC_SET_ANY);
    x1.subscribe(home,this, PC_SET_ANY);
  }

  template <class View0, class View1>
  forceinline
  ReSubset<View0,View1>::ReSubset(Space* home, bool share, ReSubset& p)
    : Propagator(home,share,p) {
    x0.update(home,share,p.x0);
    x1.update(home,share,p.x1);
    b.update(home,share,p.b);
  }

  template <class View0, class View1>
  PropCost
  ReSubset<View0,View1>::cost(ModEventDelta) const {
    return PC_TERNARY_LO;
  }

  template <class View0, class View1>
  Support::Symbol
  ReSubset<View0,View1>::ati(void) {
    return Reflection::mangle<View0,View1>("Gecode::Set::Rel::ReSubset");
  }

  template <class View0, class View1>
  Reflection::ActorSpec
  ReSubset<View0,View1>::spec(const Space* home,
                              Reflection::VarMap& m) const {
    Reflection::ActorSpec s(ati());
    return s << x0.spec(home, m)
             << x1.spec(home, m)
             << b.spec(home, m);
  }

  template <class View0, class View1>
  void
  ReSubset<View0,View1>::post(Space* home, Reflection::VarMap& vars,
                              const Reflection::ActorSpec& spec) {
    spec.checkArity(3);
    View0 x0(home, vars, spec[0]);
    View1 x1(home, vars, spec[1]);
    Gecode::Int::BoolView b(home, vars, spec[2]);
    (void) new (home) ReSubset(home,x0,x1,b);
  }

  template <class View0, class View1>
  size_t
  ReSubset<View0,View1>::dispose(Space* home) {
    assert(!home->failed());
    b.cancel(home,this, Gecode::Int::PC_INT_VAL);
    x0.cancel(home,this, PC_SET_ANY);
    x1.cancel(home,this, PC_SET_ANY);
    (void) Propagator::dispose(home);
    return sizeof(*this);
  }

  template <class View0, class View1>
  ExecStatus
  ReSubset<View0,View1>::post(Space* home, View0 x0, View1 x1,
                              Gecode::Int::BoolView b) {
    (void) new (home) ReSubset<View0,View1>(home,x0,x1,b);
    return ES_OK;
  }

  template <class View0, class View1>
  Actor*
  ReSubset<View0,View1>::copy(Space* home, bool share) {
    return new (home) ReSubset<View0,View1>(home,share,*this);
  }

  template <class View0, class View1>
  ExecStatus
  ReSubset<View0,View1>::propagate(Space* home, ModEventDelta) {
    if (b.one())
      GECODE_REWRITE(this,(SubSet<View0,View1>::post(home,x0,x1)));
    if (b.zero())
      GECODE_REWRITE(this,(NoSubSet<View0,View1>::post(home,x0,x1)));

    // check whether cardinalities still allow subset
    if (x0.cardMin() > x1.cardMax()) {
      GECODE_ME_CHECK(b.zero_none(home));
      return ES_SUBSUMED(this,home);
    }

    // check lub(x0) subset glb(x1)
    {
      LubRanges<View0> x0ub(x0);
      GlbRanges<View1> x1lb(x1);
      Iter::Ranges::Diff<LubRanges<View0>,GlbRanges<View1> > d(x0ub,x1lb);
      if (!d()) {
        GECODE_ME_CHECK(b.one_none(home));
        return ES_SUBSUMED(this,home);
      }
    }

    // check glb(x0) subset lub(x1)
    {
      GlbRanges<View0> x0lb(x0);
      LubRanges<View1> x1ub(x1);
      Iter::Ranges::Diff<GlbRanges<View0>,LubRanges<View1> > d(x0lb,x1ub);
      if (d()) {
        GECODE_ME_CHECK(b.zero_none(home));
        return ES_SUBSUMED(this,home);
      } else if (x0.assigned() && x1.assigned()) {
        GECODE_ME_CHECK(b.one_none(home));
        return ES_SUBSUMED(this,home);
      }
    }
    
    if (x0.cardMin() > 0) {
      LubRanges<View0> x0ub(x0);
      LubRanges<View1> x1ub(x1);
      Iter::Ranges::Inter<LubRanges<View0>,LubRanges<View1> >
        i(x0ub,x1ub);
      if (!i()) {
        GECODE_ME_CHECK(b.zero_none(home));
        return ES_SUBSUMED(this,home);
      }
    }
    
    return ES_FIX;
  }

}}}

// STATISTICS: set-prop