gecoder-with-gecode 0.9.0 → 0.9.1

data/ext/gecode-2.2.0/gecode/int/arithmetic.cc

@@ -0,0 +1,237 @@
+/* -*- mode: C++; c-basic-offset: 2; indent-tabs-mode: nil -*- */
+/*
+ *  Main authors:
+ *     Christian Schulte <schulte@gecode.org>
+ *
+ *  Copyright:
+ *     Christian Schulte, 2002
+ *
+ *  Last modified:
+ *     $Date: 2008-07-11 10:52:51 +0200 (Fri, 11 Jul 2008) $ by $Author: tack $
+ *     $Revision: 7352 $
+ *
+ *  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.
+ *
+ */
+
+#include "gecode/int/arithmetic.hh"
+
+namespace Gecode {
+
+  using namespace Int;
+
+  void
+  abs(Space* home, IntVar x0, IntVar x1, IntConLevel icl, PropKind) {
+    if (home->failed()) return;
+    if (icl == ICL_DOM) {
+      GECODE_ES_FAIL(home,Arithmetic::AbsDom<IntView>::post(home,x0,x1));
+    } else {
+      GECODE_ES_FAIL(home,Arithmetic::AbsBnd<IntView>::post(home,x0,x1));
+    }
+  }
+
+
+  void
+  max(Space* home, IntVar x0, IntVar x1, IntVar x2, 
+      IntConLevel icl, PropKind) {
+    if (home->failed()) return;
+    if (icl == ICL_DOM) {
+      GECODE_ES_FAIL(home,Arithmetic::MaxDom<IntView>::post(home,x0,x1,x2));
+    } else {
+      GECODE_ES_FAIL(home,Arithmetic::MaxBnd<IntView>::post(home,x0,x1,x2));
+    }
+  }
+
+  void
+  max(Space* home, const IntVarArgs& x, IntVar y,
+      IntConLevel icl, PropKind) {
+    if (x.size() == 0)
+      throw TooFewArguments("Int::max");
+    if (home->failed()) return;
+    ViewArray<IntView> xv(home,x);
+    if (icl == ICL_DOM) {
+      GECODE_ES_FAIL(home,Arithmetic::NaryMaxDom<IntView>::post(home,xv,y));
+    } else {
+      GECODE_ES_FAIL(home,Arithmetic::NaryMaxBnd<IntView>::post(home,xv,y));
+    }
+  }
+
+
+  void
+  min(Space* home, IntVar x0, IntVar x1, IntVar x2, 
+      IntConLevel icl, PropKind) {
+    if (home->failed()) return;
+    MinusView m0(x0); MinusView m1(x1); MinusView m2(x2);
+    if (icl == ICL_DOM) {
+      GECODE_ES_FAIL(home,Arithmetic::MaxDom<MinusView>::post(home,m0,m1,m2));
+    } else {
+      GECODE_ES_FAIL(home,Arithmetic::MaxBnd<MinusView>::post(home,m0,m1,m2));
+    }
+  }
+
+  void
+  min(Space* home, const IntVarArgs& x, IntVar y, 
+      IntConLevel icl, PropKind) {
+    if (x.size() == 0)
+      throw TooFewArguments("Int::min");
+    if (home->failed()) return;
+    ViewArray<MinusView> m(home,x.size());
+    for (int i=x.size(); i--; )
+      m[i].init(x[i]);
+    MinusView my(y);
+    if (icl == ICL_DOM) {
+      GECODE_ES_FAIL(home,Arithmetic::NaryMaxDom<MinusView>::post(home,m,my));
+    } else {
+      GECODE_ES_FAIL(home,Arithmetic::NaryMaxBnd<MinusView>::post(home,m,my));
+    }
+  }
+
+
+  void
+  mult(Space* home, IntVar x0, IntVar x1, IntVar x2, 
+       IntConLevel icl, PropKind) {
+    if (home->failed()) return;
+    if (icl == ICL_DOM) {
+      GECODE_ES_FAIL(home,Arithmetic::MultDom<IntView>::post(home,x0,x1,x2));
+    } else {
+      GECODE_ES_FAIL(home,Arithmetic::MultBnd<IntView>::post(home,x0,x1,x2));
+    }
+  }
+
+
+  void
+  sqr(Space* home, IntVar x0, IntVar x1, IntConLevel icl, PropKind) {
+    if (home->failed()) return;
+    if (icl == ICL_DOM) {
+      GECODE_ES_FAIL(home,Arithmetic::SqrDom<IntView>::post(home,x0,x1));
+    } else {
+      GECODE_ES_FAIL(home,Arithmetic::SqrBnd<IntView>::post(home,x0,x1));
+    }
+  }
+
+  void
+  sqrt(Space* home, IntVar x0, IntVar x1, IntConLevel icl, PropKind) {
+    if (home->failed()) return;
+    if (icl == ICL_DOM) {
+      GECODE_ES_FAIL(home,Arithmetic::SqrtDom<IntView>::post(home,x0,x1));
+    } else {
+      GECODE_ES_FAIL(home,Arithmetic::SqrtBnd<IntView>::post(home,x0,x1));
+    }
+  }
+
+  void
+  divmod(Space* home, IntVar x0, IntVar x1, IntVar x2, IntVar x3,
+         IntConLevel, PropKind) {
+    if (home->failed()) return;
+
+    IntVar prod(home, Int::Limits::min, Int::Limits::max);
+    GECODE_ES_FAIL(home, 
+                   Arithmetic::MultBnd<IntView>::post(home,x1,x2,prod));
+    Linear::Term<IntView> t[3];
+    t[0].a = 1; t[0].x = prod;
+    t[1].a = 1; t[1].x = x3;
+    int min, max;
+    Linear::estimate(t,2,0,min,max);
+    IntView x0v(x0);
+    GECODE_ME_FAIL(home, x0v.gq(home,min));
+    GECODE_ME_FAIL(home, x0v.lq(home,max));
+    t[2].a=-1; t[2].x=x0;
+    Linear::post(home,t,3,IRT_EQ,0);
+    if (home->failed()) return;
+    IntView x1v(x1);
+    GECODE_ES_FAIL(home,
+      Arithmetic::DivMod<IntView>::post(home,x1,x3));
+  }
+
+  void
+  div(Space* home, IntVar x0, IntVar x1, IntVar x2,
+      IntConLevel, PropKind) {
+    if (home->failed()) return;
+    GECODE_ES_FAIL(home,
+      (Arithmetic::DivBnd<IntView>::post(home,x0,x1,x2)));
+  }
+
+  void
+  mod(Space* home, IntVar x0, IntVar x1, IntVar x2,
+      IntConLevel icl, PropKind pk) {
+    if (home->failed()) return;
+    IntVar _div(home, Int::Limits::min, Int::Limits::max);
+    divmod(home, x0, x1, _div, x2, icl, pk);
+  }
+
+  namespace {
+    GECODE_REGISTER1(Arithmetic::AbsBnd<IntView>);
+    GECODE_REGISTER1(Arithmetic::AbsDom<IntView>);
+
+    GECODE_REGISTER1(Arithmetic::MaxBnd<IntView>);
+    GECODE_REGISTER1(Arithmetic::MaxDom<IntView>);
+    GECODE_REGISTER1(Arithmetic::MaxBnd<MinusView>);
+    GECODE_REGISTER1(Arithmetic::MaxDom<MinusView>);
+    GECODE_REGISTER1(Arithmetic::NaryMaxBnd<IntView>);
+    GECODE_REGISTER1(Arithmetic::NaryMaxDom<IntView>);
+    GECODE_REGISTER1(Arithmetic::NaryMaxBnd<MinusView>);
+    GECODE_REGISTER1(Arithmetic::NaryMaxDom<MinusView>);
+
+    GECODE_REGISTER2(Arithmetic::MultZeroOne<IntView,PC_INT_BND>);
+    GECODE_REGISTER2(Arithmetic::MultZeroOne<IntView,PC_INT_DOM>);
+    GECODE_REGISTER4(Arithmetic::MultPlusBnd<double,IntView,IntView,IntView>);
+    GECODE_REGISTER4(Arithmetic::MultPlusBnd<double,MinusView,IntView,MinusView>);
+    GECODE_REGISTER4(Arithmetic::MultPlusBnd<double,MinusView,MinusView,IntView>);
+    GECODE_REGISTER4(Arithmetic::MultPlusBnd<int,IntView,IntView,IntView>);
+    GECODE_REGISTER4(Arithmetic::MultPlusBnd<int,MinusView,IntView,MinusView>);
+    GECODE_REGISTER4(Arithmetic::MultPlusBnd<int,MinusView,MinusView,IntView>);
+    GECODE_REGISTER1(Arithmetic::MultBnd<IntView>);
+    GECODE_REGISTER4(Arithmetic::MultPlusDom<double,IntView,IntView,IntView>);
+    GECODE_REGISTER4(Arithmetic::MultPlusDom<double,MinusView,IntView,MinusView>);
+    GECODE_REGISTER4(Arithmetic::MultPlusDom<double,MinusView,MinusView,IntView>);
+    GECODE_REGISTER4(Arithmetic::MultPlusDom<int,IntView,IntView,IntView>);
+    GECODE_REGISTER4(Arithmetic::MultPlusDom<int,MinusView,IntView,MinusView>);
+    GECODE_REGISTER4(Arithmetic::MultPlusDom<int,MinusView,MinusView,IntView>);
+    GECODE_REGISTER1(Arithmetic::MultDom<IntView>);
+
+    GECODE_REGISTER1(Arithmetic::SqrBnd<IntView>);
+    GECODE_REGISTER2(Arithmetic::SqrPlusBnd<IntView,IntView>);
+    GECODE_REGISTER2(Arithmetic::SqrPlusBnd<MinusView,IntView>);
+    GECODE_REGISTER1(Arithmetic::SqrDom<IntView>);
+    GECODE_REGISTER2(Arithmetic::SqrPlusDom<IntView,IntView>);
+    GECODE_REGISTER2(Arithmetic::SqrPlusDom<MinusView,IntView>);
+
+    GECODE_REGISTER1(Arithmetic::SqrtBnd<IntView>);
+    GECODE_REGISTER1(Arithmetic::SqrtDom<IntView>);
+    GECODE_REGISTER1(Arithmetic::DivMod<IntView>);
+    
+    GECODE_REGISTER4(Arithmetic::DivPlusBnd<int,IntView,IntView,IntView>);
+    GECODE_REGISTER4(Arithmetic::DivPlusBnd<double,IntView,IntView,IntView>);
+    GECODE_REGISTER4(Arithmetic::DivPlusBnd<double,MinusView,
+                                            MinusView,IntView>);
+    GECODE_REGISTER5(Arithmetic::DivPlusBnd<double,IntView,MinusView,
+                                            MinusView,false>);
+    GECODE_REGISTER5(Arithmetic::DivPlusBnd<double,MinusView,IntView,
+                                            MinusView,false>);
+    GECODE_REGISTER1(Arithmetic::DivBnd<IntView>);
+  }
+}
+
+// STATISTICS: int-post
+

data/ext/gecode-2.2.0/gecode/int/arithmetic.hh

@@ -0,0 +1,827 @@
+/* -*- mode: C++; c-basic-offset: 2; indent-tabs-mode: nil -*- */
+/*
+ *  Main authors:
+ *     Christian Schulte <schulte@gecode.org>
+ *     Guido Tack <tack@gecode.org>
+ *
+ *  Copyright:
+ *     Christian Schulte, 2002
+ *     Guido Tack, 2004
+ *
+ *  Last modified:
+ *     $Date: 2008-07-11 10:52:51 +0200 (Fri, 11 Jul 2008) $ by $Author: tack $
+ *     $Revision: 7352 $
+ *
+ *  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.
+ *
+ */
+
+#ifndef __GECODE_INT_ARITHMETIC_HH__
+#define __GECODE_INT_ARITHMETIC_HH__
+
+#include "gecode/int.hh"
+
+#include "gecode/int/rel.hh"
+#include "gecode/int/linear.hh"
+
+/**
+ * \namespace Gecode::Int::Arithmetic
+ * \brief Numerical (arithmetic) propagators
+ */
+
+namespace Gecode { namespace Int { namespace Arithmetic {
+
+  /**
+   * \brief Bounds consistent absolute value propagator
+   *
+   * Requires \code #include "gecode/int/arithmetic.hh" \endcode
+   * \ingroup FuncIntProp
+   */
+  template <class View>
+  class AbsBnd : public BinaryPropagator<View,PC_INT_BND> {
+  protected:
+    using BinaryPropagator<View,PC_INT_BND>::x0;
+    using BinaryPropagator<View,PC_INT_BND>::x1;
+
+    /// Constructor for cloning \a p
+    AbsBnd(Space* home, bool share, AbsBnd& p);
+    /// Constructor for posting
+    AbsBnd(Space* home, View x0, View x1);
+  public:
+
+    /// Copy propagator during cloning
+    virtual Actor* copy(Space* home, bool share);
+    /**
+     * \brief Cost function
+     *
+     * If a view has been assigned, the cost is PC_UNARY_LO.
+     * Otherwise it is PC_BINARY_LO.
+     */
+    virtual PropCost cost(ModEventDelta med) const;
+    /// Perform propagation
+    virtual ExecStatus  propagate(Space* home, ModEventDelta med);
+    /// Post bounds consistent propagator \f$ |x_0|=x_1\f$
+    static  ExecStatus  post(Space* home, View x0, View x1);
+    /// Post propagator for specification
+    static void post(Space* home, Reflection::VarMap& vars,
+                     const Reflection::ActorSpec& spec);
+    /// Specification for this propagator
+    virtual Reflection::ActorSpec spec(const Space* home,
+                                        Reflection::VarMap& m) const;
+    /// Name of this propagator
+    static Support::Symbol ati(void);
+  };
+
+  /**
+   * \brief Domain consistent absolute value propagator
+   *
+   * Requires \code #include "gecode/int/arithmetic.hh" \endcode
+   * \ingroup FuncIntProp
+   */
+  template <class View>
+  class AbsDom : public BinaryPropagator<View,PC_INT_DOM> {
+  protected:
+    using BinaryPropagator<View,PC_INT_DOM>::x0;
+    using BinaryPropagator<View,PC_INT_DOM>::x1;
+
+    /// Constructor for cloning \a p
+    AbsDom(Space* home, bool share, AbsDom& p);
+    /// Constructor for posting
+    AbsDom(Space* home, View x0, View x1);
+  public:
+    /// Copy propagator during cloning
+    virtual Actor* copy(Space* home, bool share);
+    /**
+     * \brief Cost function
+     *
+     * If a view has been assigned, the cost is PC_UNARY_LO.
+     * If in stage for bounds propagation, the cost is
+     * PC_BINARY_LO. Otherwise it is PC_BINARY_HI.
+     */
+    virtual PropCost cost(ModEventDelta med) const;
+    /// Perform propagation
+    virtual ExecStatus  propagate(Space* home, ModEventDelta med);
+    /// Post domain consistent propagator \f$ |x_0|=x_1\f$
+    static  ExecStatus  post(Space* home, View x0, View x1);
+    /// Post propagator for specification
+    static void post(Space* home, Reflection::VarMap& vars,
+                     const Reflection::ActorSpec& spec);
+    /// Specification for this propagator
+    virtual Reflection::ActorSpec spec(const Space* home,
+                                        Reflection::VarMap& m) const;
+    /// Name of this propagator
+    static Support::Symbol ati(void);
+  };
+
+
+
+  /**
+   * \brief Bounds consistent ternary maximum propagator
+   *
+   * Requires \code #include "gecode/int/arithmetic.hh" \endcode
+   * \ingroup FuncIntProp
+   */
+  template <class View>
+  class MaxBnd : public TernaryPropagator<View,PC_INT_BND> {
+  protected:
+    using TernaryPropagator<View,PC_INT_BND>::x0;
+    using TernaryPropagator<View,PC_INT_BND>::x1;
+    using TernaryPropagator<View,PC_INT_BND>::x2;
+
+    /// Constructor for cloning \a p
+    MaxBnd(Space* home, bool share, MaxBnd& p);
+    /// Constructor for posting
+    MaxBnd(Space* home, View x0, View x1, View x2);
+  public:
+    /// Constructor for rewriting \a p during cloning
+    MaxBnd(Space* home, bool share, Propagator& p, View x0, View x1, View x2);
+    /// Copy propagator during cloning
+    virtual Actor* copy(Space* home, bool share);
+    /// Perform propagation
+    virtual ExecStatus propagate(Space* home, ModEventDelta med);
+    /// Post propagator \f$ \max\{x_0,x_1\}=x_2\f$
+    static  ExecStatus post(Space* home, View x0, View x1, View x2);
+    /// Post propagator for specification
+    static void post(Space* home, Reflection::VarMap& vars,
+                     const Reflection::ActorSpec& spec);
+    /// Specification for this propagator
+    virtual Reflection::ActorSpec spec(const Space* home,
+                                       Reflection::VarMap& m) const;
+    /// Name of this propagator
+    static Support::Symbol ati(void);
+  };
+
+  /**
+   * \brief Bounds consistent n-ary maximum propagator
+   *
+   * Requires \code #include "gecode/int/arithmetic.hh" \endcode
+   * \ingroup FuncIntProp
+   */
+  template <class View>
+  class NaryMaxBnd : public NaryOnePropagator<View,PC_INT_BND> {
+  protected:
+    using NaryOnePropagator<View,PC_INT_BND>::x;
+    using NaryOnePropagator<View,PC_INT_BND>::y;
+
+    /// Constructor for cloning \a p
+    NaryMaxBnd(Space* home, bool share, NaryMaxBnd& p);
+    /// Constructor for posting
+    NaryMaxBnd(Space* home, ViewArray<View>& x, View y);
+  public:
+    /// Copy propagator during cloning
+    virtual Actor* copy(Space* home, bool share);
+    /// Perform propagation
+    virtual ExecStatus propagate(Space* home, ModEventDelta med);
+    /// Post propagator \f$ \max x=y\f$
+    static  ExecStatus post(Space* home, ViewArray<View>& x, View y);
+    /// Post propagator for specification
+    static void post(Space* home, Reflection::VarMap& vars,
+                     const Reflection::ActorSpec& spec);
+    /// Specification for this propagator
+    virtual Reflection::ActorSpec spec(const Space* home,
+                                       Reflection::VarMap& m) const;
+    /// Name of this propagator
+    static Support::Symbol ati(void);
+  };
+
+  /**
+   * \brief Domain consistent ternary maximum propagator
+   *
+   * Requires \code #include "gecode/int/arithmetic.hh" \endcode
+   * \ingroup FuncIntProp
+   */
+  template <class View>
+  class MaxDom : public 
+  MixTernaryPropagator<View,PC_INT_DOM,View,PC_INT_DOM,View,PC_INT_BND> {
+  protected:
+    using MixTernaryPropagator<View,PC_INT_DOM,View,PC_INT_DOM,View,PC_INT_BND>::x0;
+    using MixTernaryPropagator<View,PC_INT_DOM,View,PC_INT_DOM,View,PC_INT_BND>::x1;
+    using MixTernaryPropagator<View,PC_INT_DOM,View,PC_INT_DOM,View,PC_INT_BND>::x2;
+
+    /// Constructor for cloning \a p
+    MaxDom(Space* home, bool share, MaxDom& p);
+    /// Constructor for posting
+    MaxDom(Space* home, View x0, View x1, View x2);
+  public:
+    /// Constructor for rewriting \a p during cloning
+    MaxDom(Space* home, bool share, Propagator& p, View x0, View x1, View x2);
+    /// Copy propagator during cloning
+    virtual Actor* copy(Space* home, bool share);
+    /**
+     * \brief Cost function
+     *
+     * If in stage for bounds propagation, the cost is
+     * PC_TERNARY_LO. Otherwise it is PC_TERNARY_HI.
+     */
+    virtual PropCost cost(ModEventDelta med) const;
+    /// Perform propagation
+    virtual ExecStatus propagate(Space* home, ModEventDelta med);
+    /// Post propagator \f$ \max\{x_0,x_1\}=x_2\f$
+    static  ExecStatus post(Space* home, View x0, View x1, View x2);
+    /// Post propagator for specification
+    static void post(Space* home, Reflection::VarMap& vars,
+                     const Reflection::ActorSpec& spec);
+    /// Specification for this propagator
+    virtual Reflection::ActorSpec spec(const Space* home,
+                                       Reflection::VarMap& m) const;
+    /// Name of this propagator
+    static Support::Symbol ati(void);
+  };
+
+  /**
+   * \brief Domain consistent n-ary maximum propagator
+   *
+   * Requires \code #include "gecode/int/arithmetic.hh" \endcode
+   * \ingroup FuncIntProp
+   */
+  template <class View>
+  class NaryMaxDom 
+    : public MixNaryOnePropagator<View,PC_INT_DOM,View,PC_INT_BND> {
+  protected:
+    using MixNaryOnePropagator<View,PC_INT_DOM,View,PC_INT_BND>::x;
+    using MixNaryOnePropagator<View,PC_INT_DOM,View,PC_INT_BND>::y;
+
+    /// Constructor for cloning \a p
+    NaryMaxDom(Space* home, bool share, NaryMaxDom& p);
+    /// Constructor for posting
+    NaryMaxDom(Space* home, ViewArray<View>& x, View y);
+  public:
+    /// Copy propagator during cloning
+    virtual Actor* copy(Space* home, bool share);
+    /**
+     * \brief Cost function
+     *
+     * If in stage for bounds propagation, the cost is dynamic
+     * PC_LINEAR_LO. Otherwise it is dynamic PC_LINEAR_HI.
+     */
+    virtual PropCost cost(ModEventDelta med) const;
+    /// Perform propagation
+    virtual ExecStatus propagate(Space* home, ModEventDelta med);
+    /// Post propagator \f$ \max x=y\f$
+    static  ExecStatus post(Space* home, ViewArray<View>& x, View y);
+    /// Post propagator for specification
+    static void post(Space* home, Reflection::VarMap& vars,
+                     const Reflection::ActorSpec& spec);
+    /// Specification for this propagator
+    virtual Reflection::ActorSpec spec(const Space* home,
+                                       Reflection::VarMap& m) const;
+    /// Name of this propagator
+    static Support::Symbol ati(void);
+  };
+
+
+
+
+  /**
+   * \brief Bounds consistent positive square propagator
+   *
+   * This propagator provides multiplication for positive views only.
+   */
+  template <class VA, class VB>
+  class SqrPlusBnd : public MixBinaryPropagator<VA,PC_INT_BND,VB,PC_INT_BND> {
+  protected:
+    using MixBinaryPropagator<VA,PC_INT_BND,VB,PC_INT_BND>::x0;
+    using MixBinaryPropagator<VA,PC_INT_BND,VB,PC_INT_BND>::x1;
+    /// Constructor for posting
+    SqrPlusBnd(Space* home, VA x0, VB x1);
+    /// Constructor for cloning \a p
+    SqrPlusBnd(Space* home, bool share, SqrPlusBnd<VA,VB>& p);
+  public:
+    /// Post propagator \f$x_0\cdot x_0=x_1\f$
+    static ExecStatus post(Space* home, VA x0, VB x1);
+    /// Post propagator for specification
+    static void post(Space* home, Reflection::VarMap& vars,
+                     const Reflection::ActorSpec& spec);
+    /// Copy propagator during cloning
+    virtual Actor* copy(Space* home, bool share);
+    /// Perform propagation
+    virtual ExecStatus propagate(Space* home, ModEventDelta med);
+    /// Specification for this propagator
+    virtual Reflection::ActorSpec spec(const Space* home,
+                                       Reflection::VarMap& m) const;
+    /// Name of this propagator
+    static Support::Symbol ati(void);
+  };
+
+  /**
+   * \brief Bounds consistent square propagator
+   *
+   * Requires \code #include "gecode/int/arithmetic.hh" \endcode
+   * \ingroup FuncIntProp
+   */
+  template <class View>
+  class SqrBnd : public BinaryPropagator<View,PC_INT_BND> {
+  protected:
+    using BinaryPropagator<View,PC_INT_BND>::x0;
+    using BinaryPropagator<View,PC_INT_BND>::x1;
+
+    /// Constructor for cloning \a p
+    SqrBnd(Space* home, bool share, SqrBnd<View>& p);
+    /// Constructor for posting
+    SqrBnd(Space* home, View x0, View x1);
+  public:
+    /// Copy propagator during cloning
+    virtual Actor* copy(Space* home, bool share);
+    /// Perform propagation
+    virtual ExecStatus propagate(Space* home, ModEventDelta med);
+    /// Post propagator for specification
+    static void post(Space* home, Reflection::VarMap& vars,
+                     const Reflection::ActorSpec& spec);
+    /// Post propagator \f$x_0\cdot x_0=x_1\f$
+    static ExecStatus post(Space* home, View x0, View x1);
+    /// Specification for this propagator
+    virtual Reflection::ActorSpec spec(const Space* home,
+                                       Reflection::VarMap& m) const;
+    /// Name of this propagator
+    static Support::Symbol ati(void);
+  };
+
+
+
+  /**
+   * \brief Domain consistent positive square propagator
+   *
+   * This propagator provides multiplication for positive views only.
+   */
+  template <class VA, class VB>
+  class SqrPlusDom : public MixBinaryPropagator<VA,PC_INT_DOM,VB,PC_INT_DOM> {
+  protected:
+    using MixBinaryPropagator<VA,PC_INT_DOM,VB,PC_INT_DOM>::x0;
+    using MixBinaryPropagator<VA,PC_INT_DOM,VB,PC_INT_DOM>::x1;
+    /// Constructor for posting
+    SqrPlusDom(Space* home, VA x0, VB x1);
+    /// Constructor for cloning \a p
+    SqrPlusDom(Space* home, bool share, SqrPlusDom<VA,VB>& p);
+  public:
+    /// Post propagator \f$x_0\cdot x_0=x_1\f$
+    static ExecStatus post(Space* home, VA x0, VB x1);
+    /// Post propagator for specification
+    static void post(Space* home, Reflection::VarMap& vars,
+                     const Reflection::ActorSpec& spec);
+    /// Copy propagator during cloning
+    virtual Actor* copy(Space* home, bool share);
+    /**
+     * \brief Cost function
+     *
+     * If a view has been assigned, the cost is PC_UNARY_LO.
+     * If in stage for bounds propagation, the cost is
+     * PC_BINARY_LO. Otherwise it is PC_BINARY_HI.
+     */
+    virtual PropCost cost(ModEventDelta med) const;
+    /// Perform propagation
+    virtual ExecStatus propagate(Space* home, ModEventDelta med);
+    /// Specification for this propagator
+    virtual Reflection::ActorSpec spec(const Space* home,
+                                       Reflection::VarMap& m) const;
+    /// Name of this propagator
+    static Support::Symbol ati(void);
+  };
+
+  /**
+   * \brief Domain consistent square propagator
+   *
+   * Requires \code #include "gecode/int/arithmetic.hh" \endcode
+   * \ingroup FuncIntProp
+   */
+  template <class View>
+  class SqrDom : public BinaryPropagator<View,PC_INT_DOM> {
+  protected:
+    using BinaryPropagator<View,PC_INT_DOM>::x0;
+    using BinaryPropagator<View,PC_INT_DOM>::x1;
+
+    /// Constructor for cloning \a p
+    SqrDom(Space* home, bool share, SqrDom<View>& p);
+    /// Constructor for posting
+    SqrDom(Space* home, View x0, View x1);
+  public:
+    /// Copy propagator during cloning
+    virtual Actor* copy(Space* home, bool share);
+    /// Perform propagation
+    virtual ExecStatus propagate(Space* home, ModEventDelta med);
+    /**
+     * \brief Cost function
+     *
+     * If a view has been assigned, the cost is PC_UNARY_LO.
+     * If in stage for bounds propagation, the cost is
+     * PC_BINARY_LO. Otherwise it is PC_BINARY_HI.
+     */
+    virtual PropCost cost(ModEventDelta med) const;
+    /// Post propagator for specification
+    static void post(Space* home, Reflection::VarMap& vars,
+                     const Reflection::ActorSpec& spec);
+    /// Post propagator \f$x_0\cdot x_0=x_1\f$
+    static ExecStatus post(Space* home, View x0, View x1);
+    /// Specification for this propagator
+    virtual Reflection::ActorSpec spec(const Space* home,
+                                       Reflection::VarMap& m) const;
+    /// Name of this propagator
+    static Support::Symbol ati(void);
+  };
+
+
+
+  /**
+   * \brief Bounds consistent square root propagator
+   *
+   * Requires \code #include "gecode/int/arithmetic.hh" \endcode
+   * \ingroup FuncIntProp
+   */
+  template <class View>
+  class SqrtBnd : public BinaryPropagator<View,PC_INT_BND> {
+  protected:
+    using BinaryPropagator<View,PC_INT_BND>::x0;
+    using BinaryPropagator<View,PC_INT_BND>::x1;
+
+    /// Constructor for cloning \a p
+    SqrtBnd(Space* home, bool share, SqrtBnd<View>& p);
+    /// Constructor for posting
+    SqrtBnd(Space* home, View x0, View x1);
+  public:
+    /// Copy propagator during cloning
+    virtual Actor* copy(Space* home, bool share);
+    /// Perform propagation
+    virtual ExecStatus propagate(Space* home, ModEventDelta med);
+    /// Post propagator for specification
+    static void post(Space* home, Reflection::VarMap& vars,
+                     const Reflection::ActorSpec& spec);
+    /// Post propagator \f$\lfloor\sqrt{x_0}\rfloor=x_1\f$
+    static ExecStatus post(Space* home, View x0, View x1);
+    /// Specification for this propagator
+    virtual Reflection::ActorSpec spec(const Space* home,
+                                       Reflection::VarMap& m) const;
+    /// Name of this propagator
+    static Support::Symbol ati(void);
+  };
+
+  /**
+   * \brief Domain consistent square root propagator
+   *
+   * Requires \code #include "gecode/int/arithmetic.hh" \endcode
+   * \ingroup FuncIntProp
+   */
+  template <class View>
+  class SqrtDom : public BinaryPropagator<View,PC_INT_DOM> {
+  protected:
+    using BinaryPropagator<View,PC_INT_DOM>::x0;
+    using BinaryPropagator<View,PC_INT_DOM>::x1;
+
+    /// Constructor for cloning \a p
+    SqrtDom(Space* home, bool share, SqrtDom<View>& p);
+    /// Constructor for posting
+    SqrtDom(Space* home, View x0, View x1);
+  public:
+    /// Copy propagator during cloning
+    virtual Actor* copy(Space* home, bool share);
+    /// Perform propagation
+    virtual ExecStatus propagate(Space* home, ModEventDelta med);
+    /**
+     * \brief Cost function
+     *
+     * If a view has been assigned, the cost is PC_UNARY_LO.
+     * If in stage for bounds propagation, the cost is
+     * PC_BINARY_LO. Otherwise it is PC_BINARY_HI.
+     */
+    virtual PropCost cost(ModEventDelta med) const;
+    /// Post propagator for specification
+    static void post(Space* home, Reflection::VarMap& vars,
+                     const Reflection::ActorSpec& spec);
+    /// Post propagator \f$\lfloor\sqrt{x_0}\rfloor=x_1\f$
+    static ExecStatus post(Space* home, View x0, View x1);
+    /// Specification for this propagator
+    virtual Reflection::ActorSpec spec(const Space* home,
+                                       Reflection::VarMap& m) const;
+    /// Name of this propagator
+    static Support::Symbol ati(void);
+  };
+
+
+
+  /**
+   * \brief Bounds or domain consistent propagator for \f$x_0\times x_1=x_0\f$
+   *
+   * Requires \code #include "gecode/int/arithmetic.hh" \endcode
+   * \ingroup FuncIntProp
+   */
+  template <class View, PropCond pc>
+  class MultZeroOne : public BinaryPropagator<View,pc> {
+  protected:
+    using BinaryPropagator<View,pc>::x0;
+    using BinaryPropagator<View,pc>::x1;
+
+    /// Constructor for cloning \a p
+    MultZeroOne(Space* home, bool share, MultZeroOne<View,pc>& p);
+    /// Constructor for posting
+    MultZeroOne(Space* home, View x0, View x1);
+    /// Test whether \a x is equal to \a n
+    static RelTest equal(View x, int n);
+  public:
+    /// Copy propagator during cloning
+    virtual Actor* copy(Space* home, bool share);
+    /// Perform propagation
+    virtual ExecStatus propagate(Space* home, ModEventDelta med);
+    /// Post propagator for specification
+    static void post(Space* home, Reflection::VarMap& vars,
+                     const Reflection::ActorSpec& spec);
+    /// Post propagator \f$x_0\cdot x_1=x_0\f$
+    static ExecStatus post(Space* home, View x0, View x1);
+    /// Specification for this propagator
+    virtual Reflection::ActorSpec spec(const Space* home,
+                                       Reflection::VarMap& m) const;
+    /// Name of this propagator
+    static Support::Symbol ati(void);
+  };
+
+
+
+  /**
+   * \brief Bounds consistent positive multiplication propagator
+   *
+   * This propagator provides multiplication for positive views only.
+   */
+  template <class Val, class VA, class VB, class VC>
+  class MultPlusBnd : 
+    public MixTernaryPropagator<VA,PC_INT_BND,VB,PC_INT_BND,VC,PC_INT_BND> {
+  protected:
+    using MixTernaryPropagator<VA,PC_INT_BND,VB,PC_INT_BND,VC,PC_INT_BND>::x0;
+    using MixTernaryPropagator<VA,PC_INT_BND,VB,PC_INT_BND,VC,PC_INT_BND>::x1;
+    using MixTernaryPropagator<VA,PC_INT_BND,VB,PC_INT_BND,VC,PC_INT_BND>::x2;
+  public:
+    /// Constructor for posting
+    MultPlusBnd(Space* home, VA x0, VB x1, VC x2);
+    /// Constructor for cloning \a p
+    MultPlusBnd(Space* home, bool share, MultPlusBnd<Val,VA,VB,VC>& p);
+    /// Post propagator \f$x_0\cdot x_1=x_2\f$
+    static ExecStatus post(Space* home, VA x0, VB x1, VC x2);
+    /// Post propagator for specification
+    static void post(Space* home, Reflection::VarMap& vars,
+                     const Reflection::ActorSpec& spec);
+    /// Copy propagator during cloning
+    virtual Actor* copy(Space* home, bool share);
+    /// Perform propagation
+    virtual ExecStatus propagate(Space* home, ModEventDelta med);
+    /// Specification for this propagator
+    virtual Reflection::ActorSpec spec(const Space* home,
+                                       Reflection::VarMap& m) const;
+    /// Name of this propagator
+    static Support::Symbol ati(void);
+  };
+
+  /**
+   * \brief Bounds consistent multiplication propagator
+   *
+   * Requires \code #include "gecode/int/arithmetic.hh" \endcode
+   *
+   * \ingroup FuncIntProp
+   */
+  template <class View>
+  class MultBnd : public TernaryPropagator<View,PC_INT_BND> {
+  protected:
+    using TernaryPropagator<View,PC_INT_BND>::x0;
+    using TernaryPropagator<View,PC_INT_BND>::x1;
+    using TernaryPropagator<View,PC_INT_BND>::x2;
+
+    /// Constructor for cloning \a p
+    MultBnd(Space* home, bool share, MultBnd<View>& p);
+  public:
+    /// Constructor for posting
+    MultBnd(Space* home, View x0, View x1, View x2);
+    /// Post propagator \f$x_0\cdot x_1=x_2\f$
+    static  ExecStatus post(Space* home, View x0, View x1, View x2);
+    /// Post propagator for specification
+    static void post(Space* home, Reflection::VarMap& vars,
+                     const Reflection::ActorSpec& spec);
+    /// Copy propagator during cloning
+    virtual Actor* copy(Space* home, bool share);
+    /// Perform propagation
+    virtual ExecStatus propagate(Space* home, ModEventDelta med);
+    /// Specification for this propagator
+    virtual Reflection::ActorSpec spec(const Space* home,
+                                        Reflection::VarMap& m) const;
+    /// Name of this propagator
+    static Support::Symbol ati(void);
+  };
+
+
+
+  /**
+   * \brief Domain consistent positive multiplication propagator
+   *
+   * This propagator provides multiplication for positive views only.
+   */
+  template <class Val, class VA, class VB, class VC>
+  class MultPlusDom : 
+    public MixTernaryPropagator<VA,PC_INT_DOM,VB,PC_INT_DOM,VC,PC_INT_DOM> {
+  protected:
+    using MixTernaryPropagator<VA,PC_INT_DOM,VB,PC_INT_DOM,VC,PC_INT_DOM>::x0;
+    using MixTernaryPropagator<VA,PC_INT_DOM,VB,PC_INT_DOM,VC,PC_INT_DOM>::x1;
+    using MixTernaryPropagator<VA,PC_INT_DOM,VB,PC_INT_DOM,VC,PC_INT_DOM>::x2;
+  public:
+    /// Constructor for posting
+    MultPlusDom(Space* home, VA x0, VB x1, VC x2);
+    /// Constructor for cloning \a p
+    MultPlusDom(Space* home, bool share, MultPlusDom<Val,VA,VB,VC>& p);
+    /// Post propagator \f$x_0\cdot x_1=x_2\f$
+    static ExecStatus post(Space* home, VA x0, VB x1, VC x2);
+    /// Post propagator for specification
+    static void post(Space* home, Reflection::VarMap& vars,
+                     const Reflection::ActorSpec& spec);
+    /// Copy propagator during cloning
+    virtual Actor* copy(Space* home, bool share);
+    /**
+     * \brief Cost function
+     *
+     * If in stage for bounds propagation, the cost is
+     * PC_TERNARY_LO. Otherwise it is PC_TERNARY_HI.
+     */
+    virtual PropCost cost(ModEventDelta med) const;
+    /// Perform propagation
+    virtual ExecStatus propagate(Space* home, ModEventDelta med);
+    /// Specification for this propagator
+    virtual Reflection::ActorSpec spec(const Space* home,
+                                       Reflection::VarMap& m) const;
+    /// Name of this propagator
+    static Support::Symbol ati(void);
+  };
+
+  /**
+   * \brief Domain consistent multiplication propagator
+   *
+   * Requires \code #include "gecode/int/arithmetic.hh" \endcode
+   *
+   * \ingroup FuncIntProp
+   */
+  template <class View>
+  class MultDom : public TernaryPropagator<View,PC_INT_DOM> {
+  protected:
+    using TernaryPropagator<View,PC_INT_DOM>::x0;
+    using TernaryPropagator<View,PC_INT_DOM>::x1;
+    using TernaryPropagator<View,PC_INT_DOM>::x2;
+
+    /// Constructor for cloning \a p
+    MultDom(Space* home, bool share, MultDom<View>& p);
+  public:
+    /// Constructor for posting
+    MultDom(Space* home, View x0, View x1, View x2);
+    /// Post propagator \f$x_0\cdot x_1=x_2\f$
+    static  ExecStatus post(Space* home, View x0, View x1, View x2);
+    /// Post propagator for specification
+    static void post(Space* home, Reflection::VarMap& vars,
+                     const Reflection::ActorSpec& spec);
+    /// Copy propagator during cloning
+    virtual Actor* copy(Space* home, bool share);
+    /**
+     * \brief Cost function
+     *
+     * If in stage for bounds propagation, the cost is
+     * PC_TERNARY_LO. Otherwise it is PC_TERNARY_HI.
+     */
+    virtual PropCost cost(ModEventDelta med) const;
+    /// Perform propagation
+    virtual ExecStatus propagate(Space* home, ModEventDelta med);
+    /// Specification for this propagator
+    virtual Reflection::ActorSpec spec(const Space* home,
+                                        Reflection::VarMap& m) const;
+    /// Name of this propagator
+    static Support::Symbol ati(void);
+  };
+
+  /**
+   * \brief Bounds consistent positive division propagator
+   *
+   * This propagator provides division for positive views only.
+   */
+  template <class Val, class VA, class VB, class VC, bool towardsMinInf=true>
+  class DivPlusBnd : 
+    public MixTernaryPropagator<VA,PC_INT_BND,VB,PC_INT_BND,VC,PC_INT_BND> {
+  protected:
+    using MixTernaryPropagator<VA,PC_INT_BND,VB,PC_INT_BND,VC,PC_INT_BND>::x0;
+    using MixTernaryPropagator<VA,PC_INT_BND,VB,PC_INT_BND,VC,PC_INT_BND>::x1;
+    using MixTernaryPropagator<VA,PC_INT_BND,VB,PC_INT_BND,VC,PC_INT_BND>::x2;
+  public:
+    /// Constructor for posting
+    DivPlusBnd(Space* home, VA x0, VB x1, VC x2);
+    /// Constructor for cloning \a p
+    DivPlusBnd(Space* home, bool share,
+               DivPlusBnd<Val,VA,VB,VC,towardsMinInf>& p);
+    /// Post propagator \f$x_0\mathrm{div} x_1=x_2\f$ (rounding to \f$-\infty\f$)
+    static ExecStatus post(Space* home, VA x0, VB x1, VC x2);
+    /// Post propagator for specification
+    static void post(Space* home, Reflection::VarMap& vars,
+                     const Reflection::ActorSpec& spec);
+    /// Copy propagator during cloning
+    virtual Actor* copy(Space* home, bool share);
+    /// Perform propagation
+    virtual ExecStatus propagate(Space* home, ModEventDelta med);
+    /// Specification for this propagator
+    virtual Reflection::ActorSpec spec(const Space* home,
+                                       Reflection::VarMap& m) const;
+    /// Name of this propagator
+    static Support::Symbol ati(void);
+  };
+
+  /**
+   * \brief Bounds consistent division propagator
+   *
+   * Requires \code #include "gecode/int/arithmetic.hh" \endcode
+   *
+   * \ingroup FuncIntProp
+   */
+  template <class View>
+  class DivBnd : public TernaryPropagator<View,PC_INT_BND> {
+  protected:
+    using TernaryPropagator<View,PC_INT_BND>::x0;
+    using TernaryPropagator<View,PC_INT_BND>::x1;
+    using TernaryPropagator<View,PC_INT_BND>::x2;
+
+    /// Constructor for cloning \a p
+    DivBnd(Space* home, bool share, DivBnd<View>& p);
+  public:
+    /// Constructor for posting
+    DivBnd(Space* home, View x0, View x1, View x2);
+    /// Post propagator \f$x_0\mathrm{div} x_1=x_2\f$ (rounding to \f$-\infty\f$)
+    static  ExecStatus post(Space* home, View x0, View x1, View x2);
+    /// Post propagator for specification
+    static void post(Space* home, Reflection::VarMap& vars,
+                     const Reflection::ActorSpec& spec);
+    /// Copy propagator during cloning
+    virtual Actor* copy(Space* home, bool share);
+    /// Perform propagation
+    virtual ExecStatus propagate(Space* home, ModEventDelta med);
+    /// Specification for this propagator
+    virtual Reflection::ActorSpec spec(const Space* home,
+                                        Reflection::VarMap& m) const;
+    /// Name of this propagator
+    static Support::Symbol ati(void);
+  };
+  
+  /**
+   * \brief Integer division/modulo propagator
+   *
+   * This propagator implements the relation between divisor and
+   * modulo of an integer division.
+   *
+   * Requires \code #include "gecode/int/arithmetic.hh" \endcode
+   *
+   * \ingroup FuncIntProp
+   */
+  template <class View>
+  class DivMod : public BinaryPropagator<View,PC_INT_BND> {
+  protected:
+    using BinaryPropagator<View,PC_INT_BND>::x0;
+    using BinaryPropagator<View,PC_INT_BND>::x1;
+    
+    /// Constructor for cloning \a p
+    DivMod(Space* home, bool share, DivMod<View>& p);
+  public:
+    /// Constructor for posting
+    DivMod(Space* home, View x0, View x1);
+    /// Post propagator \f$x_0\neq 0 \land (x_1\neq 0\Rightarrow x_0\times x_1>0) \land \mathrm{abs}(x_1)<\mathrm{abs}(x_0)\f$
+    static  ExecStatus post(Space* home, View x0, View x1);
+    /// Post propagator for specification
+    static void post(Space* home, Reflection::VarMap& vars,
+                     const Reflection::ActorSpec& spec);
+    /// Copy propagator during cloning
+    virtual Actor* copy(Space* home, bool share);
+    /// Perform propagation
+    virtual ExecStatus propagate(Space* home, ModEventDelta med);
+    /// Specification for this propagator
+    virtual Reflection::ActorSpec spec(const Space* home,
+                                       Reflection::VarMap& m) const;
+    /// Name of this propagator
+    static Support::Symbol ati(void);
+  };
+
+}}}
+
+#include "gecode/int/arithmetic/abs.icc"
+#include "gecode/int/arithmetic/max.icc"
+#include "gecode/int/arithmetic/sqr.icc"
+#include "gecode/int/arithmetic/sqrt.icc"
+#include "gecode/int/arithmetic/mult.icc"
+#include "gecode/int/arithmetic/divmod.icc"
+
+#endif
+
+// STATISTICS: int-prop
+