gecoder-with-gecode 0.7.1 → 0.8.0

data/ext/gecode-1.3.1/gecode/int/linear/bool-view.icc

@@ -1,311 +0,0 @@
-/*
- *  Main authors:
- *     Christian Schulte <schulte@gecode.org>
- *
- *  Copyright:
- *     Christian Schulte, 2004
- *
- *  Last modified:
- *     $Date: 2006-09-07 11:03:52 +0200 (Thu, 07 Sep 2006) $ by $Author: schulte $
- *     $Revision: 3609 $
- *
- *  This file is part of Gecode, the generic constraint
- *  development environment:
- *     http://www.gecode.org
- *
- *  See the file "LICENSE" for information on usage and
- *  redistribution of this file, and for a
- *     DISCLAIMER OF ALL WARRANTIES.
- *
- */
-
-namespace Gecode { namespace Int { namespace Linear {
-
-  /*
-   * Base-class
-   *
-   */
-  template <class XV, class YV>
-  LinBoolView<XV,YV>::LinBoolView(Space* home,
-				  ViewArray<XV>& x0, YV y0, int c0)
-    :  Propagator(home), x(x0), y(y0), c(c0) {
-    x.subscribe(home,this,PC_INT_VAL);
-    y.subscribe(home,this,PC_INT_BND);
-  }
-
-  template <class XV, class YV>
-  size_t
-  LinBoolView<XV,YV>::dispose(Space* home) {
-    assert(!home->failed());
-    x.cancel(home,this,PC_INT_VAL);
-    y.cancel(home,this,PC_INT_BND);
-    (void) Propagator::dispose(home);
-    return sizeof(*this);
-  }
-
-  template <class XV, class YV>
-  forceinline
-  LinBoolView<XV,YV>::LinBoolView(Space* home, bool share, LinBoolView& p)
-    : Propagator(home,share,p), c(p.c) {
-    x.update(home,share,p.x);
-    y.update(home,share,p.y);
-  }
-
-  template <class XV, class YV>
-  PropCost
-  LinBoolView<XV,YV>::cost(void) const {
-    return cost_lo(x.size(),PC_LINEAR_LO);
-  }
-
-
-
-  /*
-   * Equality propagator
-   *
-   */
-  template <class XV, class YV>
-  forceinline
-  EqBoolView<XV,YV>::EqBoolView(Space* home, ViewArray<XV>& x, YV y, int c)
-    : LinBoolView<XV,YV>(home,x,y,c) {}
-
-  template <class XV, class YV>
-  ExecStatus
-  EqBoolView<XV,YV>::post(Space* home, ViewArray<XV>& x, YV y, int c) {
-    if (y.assigned())
-      return EqBoolInt<XV>::post(home,x,y.val()+c);
-    int n = x.size();
-    for (int i = n; i--; )
-      if (x[i].one()) {
-	x[i]=x[--n]; c--;
-      } else if (x[i].zero()) {
-	x[i]=x[--n];
-      }
-    x.size(n);
-    GECODE_ME_CHECK(y.lq(home,n-c));
-    GECODE_ME_CHECK(y.gq(home,-c));
-    if (n == 0)
-      return ES_OK;
-    if (y.min()+c == n) {
-      assert(y.assigned());
-      for (int i = n; i--; )
-	x[i].t_one_none(home);
-      return ES_OK;
-    }
-    if (y.max()+c == 0) {
-      assert(y.assigned());
-      for (int i = n; i--; )
-	x[i].t_zero_none(home);
-      return ES_OK;
-    }
-    (void) new (home) EqBoolView<XV,YV>(home,x,y,c);
-    return ES_OK;
-  }
-
-  template <class XV, class YV>
-  forceinline
-  EqBoolView<XV,YV>::EqBoolView(Space* home, bool share, EqBoolView<XV,YV>& p)
-    : LinBoolView<XV,YV>(home,share,p) {}
-
-  template <class XV, class YV>
-  Actor*
-  EqBoolView<XV,YV>::copy(Space* home, bool share) {
-    return new (home) EqBoolView<XV,YV>(home,share,*this);
-  }
-
-  template <class XV, class YV>
-  ExecStatus
-  EqBoolView<XV,YV>::propagate(Space* home) {
-    int n = x.size();
-    for (int i = n; i--; )
-      if (x[i].one()) {
-	x[i]=x[--n]; c--;
-      } else if (x[i].zero()) {
-	x[i]=x[--n];
-      }
-    x.size(n);
-    GECODE_ME_CHECK(y.lq(home,n-c));
-    GECODE_ME_CHECK(y.gq(home,-c));
-    if (n == 0)
-      return ES_SUBSUMED;
-    if (y.min()+c == n) {
-      assert(y.assigned());
-      for (int i = n; i--; )
-	x[i].t_one_none(home);
-      return ES_SUBSUMED;
-    }
-    if (y.max()+c == 0) {
-      assert(y.assigned());
-      for (int i = n; i--; )
-	x[i].t_zero_none(home);
-      return ES_SUBSUMED;
-    }
-    if (y.assigned()) {
-      GECODE_ES_CHECK(EqBoolInt<XV>::post(home,x,y.val()+c));
-      return ES_SUBSUMED;
-    }
-    return ES_FIX;
-  }
-
-
-  /*
-   * Disequality propagator
-   *
-   */
-  template <class XV, class YV>
-  forceinline
-  NqBoolView<XV,YV>::NqBoolView(Space* home, ViewArray<XV>& x, YV y, int c)
-    : LinBoolView<XV,YV>(home,x,y,c) {}
-
-  template <class XV, class YV>
-  ExecStatus
-  NqBoolView<XV,YV>::post(Space* home, ViewArray<XV>& x, YV y, int c) {
-    if (y.assigned())
-      return NqBoolInt<XV>::post(home,x,y.val()+c);
-    int n = x.size();
-    for (int i = n; i--; )
-      if (x[i].one()) {
-	x[i]=x[--n]; c--;
-      } else if (x[i].zero()) {
-	x[i]=x[--n];
-      }
-    x.size(n);
-    if ((n-c < y.min() ) || (-c > y.max()))
-      return ES_OK;
-    if (n == 0) {
-      GECODE_ME_CHECK(y.nq(home,-c));
-      return ES_OK;
-    }
-    if ((n == 1) && y.assigned()) {
-      if (y.val()+c == 1) {
-	x[0].t_zero_none(home);
-      } else {
-	assert(y.val()+c == 0);
-	x[0].t_one_none(home);
-      }
-      return ES_OK;
-    }
-    (void) new (home) NqBoolView<XV,YV>(home,x,y,c);
-    return ES_OK;
-  }
-
-
-  template <class XV, class YV>
-  forceinline
-  NqBoolView<XV,YV>::NqBoolView(Space* home, bool share, NqBoolView<XV,YV>& p)
-    : LinBoolView<XV,YV>(home,share,p) {}
-
-  template <class XV, class YV>
-  Actor*
-  NqBoolView<XV,YV>::copy(Space* home, bool share) {
-    return new (home) NqBoolView<XV,YV>(home,share,*this);
-  }
-
-
-  template <class XV, class YV>
-  ExecStatus
-  NqBoolView<XV,YV>::propagate(Space* home) {
-    int n = x.size();
-    for (int i = n; i--; )
-      if (x[i].one()) {
-	x[i]=x[--n]; c--;
-      } else if (x[i].zero()) {
-	x[i]=x[--n];
-      }
-    x.size(n);
-    if ((n-c < y.min() ) || (-c > y.max()))
-      return ES_SUBSUMED;
-    if (n == 0) {
-      GECODE_ME_CHECK(y.nq(home,-c));
-      return ES_SUBSUMED;
-    }
-    if ((n == 1) && y.assigned()) {
-      if (y.val()+c == 1) {
-	x[0].t_zero_none(home);
-      } else {
-	assert(y.val()+c == 0);
-	x[0].t_one_none(home);
-      }
-      return ES_SUBSUMED;
-    }
-    return ES_FIX;
-  }
-
-
-  /*
-   * Greater or equal propagator
-   *
-   */
-  template <class XV, class YV>
-  forceinline
-  GqBoolView<XV,YV>::GqBoolView(Space* home, ViewArray<XV>& x, YV y, int c)
-    : LinBoolView<XV,YV>(home,x,y,c) {}
-
-  template <class XV, class YV>
-  ExecStatus
-  GqBoolView<XV,YV>::post(Space* home, ViewArray<XV>& x, YV y, int c) {
-    if (y.assigned())
-      return GqBoolInt<XV>::post(home,x,y.val()+c);
-    // Eliminate assigned views
-    int n = x.size();
-    for (int i = n; i--; )
-      if (x[i].one()) {
-	x[i]=x[--n]; c--;
-      } else if (x[i].zero()) {
-	x[i]=x[--n];
-      }
-    x.size(n);
-    GECODE_ME_CHECK(y.lq(home,n-c));
-    if (-c >= y.max())
-      return ES_OK;
-    if (y.min()+c == n) {
-      for (int i = n; i--; )
-	x[i].t_one_none(home);
-      return ES_OK;
-    }
-    (void) new (home) GqBoolView<XV,YV>(home,x,y,c);
-    return ES_OK;
-  }
-
-
-  template <class XV, class YV>
-  forceinline
-  GqBoolView<XV,YV>::GqBoolView(Space* home, bool share, GqBoolView<XV,YV>& p)
-    : LinBoolView<XV,YV>(home,share,p) {}
-
-  template <class XV, class YV>
-  Actor*
-  GqBoolView<XV,YV>::copy(Space* home, bool share) {
-    return new (home) GqBoolView<XV,YV>(home,share,*this);
-  }
-
-
-  template <class XV, class YV>
-  ExecStatus
-  GqBoolView<XV,YV>::propagate(Space* home) {
-    int n = x.size();
-    for (int i = n; i--; )
-      if (x[i].one()) {
-	x[i]=x[--n]; c--;
-      } else if (x[i].zero()) {
-	x[i]=x[--n];
-      }
-    x.size(n);
-    GECODE_ME_CHECK(y.lq(home,n-c));
-    if (-c >= y.max())
-      return ES_SUBSUMED;
-    if (y.min()+c == n) {
-      for (int i = n; i--; )
-	x[i].t_one_none(home);
-      return ES_SUBSUMED;
-    }
-    if (y.assigned()) {
-      GECODE_ES_CHECK(GqBoolInt<XV>::post(home,x,y.val()+c));
-      return ES_SUBSUMED;
-    }
-    return ES_FIX;
-  }
-
-}}}
-
-// STATISTICS: int-prop
-

data/ext/gecode-1.3.1/gecode/int/linear/dom.icc

@@ -1,503 +0,0 @@
-/*
- *  Main authors:
- *     Christian Schulte <schulte@gecode.org>
- *
- *  Copyright:
- *     Christian Schulte, 2006
- *
- *  Last modified:
- *     $Date: 2006-08-04 16:03:26 +0200 (Fri, 04 Aug 2006) $ by $Author: schulte $
- *     $Revision: 3512 $
- *
- *  This file is part of Gecode, the generic constraint
- *  development environment:
- *     http://www.gecode.org
- *
- *  See the file "LICENSE" for information on usage and
- *  redistribution of this file, and for a
- *     DISCLAIMER OF ALL WARRANTIES.
- *
- */
-
-#include "gecode/int/linear.hh"
-
-#include "gecode/iter.hh"
-
-namespace Gecode { namespace Int { namespace Linear {
-
-  /**
-   * \brief Set for support information
-   *
-   * Records supported positions of values such that with iteration
-   * the supported values can be reconstructed.
-   *
-   */
-  class SupportSet {
-  private:
-    /// Number of bits per unsigned integer
-    static const unsigned int bpui = sizeof(unsigned int) * 8;
-    /// Array of bits
-    unsigned int* bits;
-  public:
-    /// Initialize support set with cardinality \arg n
-    void init(unsigned int n);
-    /// Record that there is support at position \arg i
-    void support(unsigned int i);
-    /// Check whether position \arg i has support
-    bool supported(unsigned int i) const;
-
-  private:
-    /// Support-based iterator
-    class ResultIter : public ViewValues<IntView> {
-    protected:
-      /// The support set used
-      const SupportSet* s;
-      /// The current position of the value
-      unsigned int p;
-    public:
-      /// Initialize iterator
-      ResultIter(const SupportSet* s0, const IntView& x);
-      /// Increment to next supported value
-      void operator++(void);
-    };
-
-  public:
-    /// Perform tell according to recorded support information on \arg x
-    ModEvent tell(Space* home, IntView& x) const;
-    /// Release memory
-    void dispose(void);
-  };
-
-  /**
-   * \brief Base-class for support-based iterator
-   *
-   */
-  template <class Val>
-  class SupportIter {
-  protected:
-    /// Integer coefficient for view
-    int           a;
-    /// Integer view
-    IntView       x;
-    /// Set of support for values in x
-    SupportSet    s;
-    /// Current value
-    int           c;
-    /// Position of current value
-    unsigned int  p;
-    /// Lower bound information for value
-    Val           l;
-    /// Upper bound information for value
-    Val           u;
-  public:
-    /// Initialize view
-    void init(int a, const IntView& x, Val l, Val u);
-    /// Record value at current position as supported
-    void support(void);
-    /// Tell back new variable domain according to support found
-    ModEvent tell(Space* home);
-    /// Finalize memory for support set
-    void dispose(void);
-  };
-
-
-  /**
-   * \brief Support-based iterator for positive view
-   *
-   */
-  template <class Val>
-  class PosSupportIter : public SupportIter<Val> {
-  private:
-    /// Iterate ranges of integer view in increasing order
-    IntVarImpFwd i;
-    // Using-declarations for dependant names
-    using SupportIter<Val>::a;
-    using SupportIter<Val>::x;
-    using SupportIter<Val>::s;
-    using SupportIter<Val>::c;
-    using SupportIter<Val>::p;
-    using SupportIter<Val>::l;
-    using SupportIter<Val>::u;
-  public:
-    /// Reset iterator to beginning and adjust \arg d accordingly
-    bool reset(Val& d);
-    /// Adjust \arg d and return true if next value still works
-    bool adjust(Val& d);
-  };
-
-
-  /**
-   * \brief Support-based iterator for negative view
-   *
-   */
-  template <class Val>
-  class NegSupportIter : public SupportIter<Val> {
-  private:
-    /// Iterate ranges of integer view in decreasing order
-    IntVarImpBwd i;
-    // Using-declarations for dependant names
-    using SupportIter<Val>::a;
-    using SupportIter<Val>::x;
-    using SupportIter<Val>::s;
-    using SupportIter<Val>::c;
-    using SupportIter<Val>::p;
-    using SupportIter<Val>::l;
-    using SupportIter<Val>::u;
-  public:
-    /// Reset iterator to beginning and adjust \arg d accordingly
-    bool reset(Val& d);
-    /// Adjust \arg d and return true if next value still works
-    bool adjust(Val& d);
-  };
-
-
-  /*
-   * Support set
-   *
-   */
-  forceinline void
-  SupportSet::init(unsigned int n) {
-    bits = Memory::bmalloc<unsigned int>((n / bpui) + 1);
-    for (unsigned int i = (n / bpui) + 1; i--; )
-      bits[i] = 0;
-  }
-  forceinline void
-  SupportSet::support(unsigned int i) {
-    unsigned int p = i / bpui;
-    bits[p] |= 1 << (i-p*bpui);
-  }
-  forceinline bool
-  SupportSet::supported(unsigned int i) const {
-    unsigned int p = i / bpui;
-    return (bits[p] & (1 << (i-p*bpui))) != 0;
-  }
-
-  forceinline
-  SupportSet::ResultIter::ResultIter(const SupportSet* s0, const IntView& x)
-    : ViewValues<IntView>(x), s(s0), p(0) {
-    while (ViewValues<IntView>::operator()() && !s->supported(p)) {
-      ViewValues<IntView>::operator++(); ++p;
-    }
-  }
-  forceinline void
-  SupportSet::ResultIter::operator++(void) {
-    do {
-      ViewValues<IntView>::operator++(); ++p;
-    } while (ViewValues<IntView>::operator()() && !s->supported(p));
-  }
-
-
-  inline ModEvent
-  SupportSet::tell(Space* home, IntView& x) const {
-    unsigned int n = x.size() / bpui;
-    // Check whether all bits are zero: failure
-    for (unsigned int i=n+1; i--; )
-      if (bits[i] != 0)
-	goto all;
-    return ME_INT_FAILED;
-  all:
-    // Check whether all bits are one: nothing changed
-    for (unsigned int i=n; i--; )
-      if (bits[i] != ~0U)
-	goto tell;
-    // Now check the bits in the last word
-    for (unsigned int i=n*bpui; i<x.size(); i++)
-      if (!supported(i))
-	goto tell;
-    return ME_INT_NONE;
-  tell:
-    {
-      ResultIter i(this,x);
-      Iter::Values::ToRanges<ResultIter> r(i);
-      return x.narrow(home,r);
-    }
-  }
-
-  forceinline void
-  SupportSet::dispose(void) {
-    Memory::free(bits);
-  }
-
-
-  /*
-   * Base-class for support-based iterator
-   *
-   */
-  template <class Val>
-  forceinline void
-  SupportIter<Val>::init(int a0, const IntView& x0, Val l0, Val u0) {
-    a=a0; x=x0; l=l0; u=u0;
-    s.init(x.size());
-  }
-  template <class Val>
-  forceinline void
-  SupportIter<Val>::support(void) {
-    s.support(p);
-  }
-  template <class Val>
-  forceinline ModEvent
-  SupportIter<Val>::tell(Space* home) {
-    return s.tell(home,x);
-  }
-  template <class Val>
-  forceinline void
-  SupportIter<Val>::dispose(void) {
-    s.dispose();
-  }
-
-
-  /*
-   * Support-based iterator for positive view
-   *
-   */
-  template <class Val>
-  forceinline bool
-  PosSupportIter<Val>::reset(Val& d) {
-    // Way too small, no value can make it big enough
-    if (d + static_cast<Val>(a)*x.max() < u)
-      return false;
-    // Restart iterator and position of values
-    i.init(x.variable()); p = 0;
-    // Skip all ranges which are too small
-    while (d + static_cast<Val>(a)*i.max() < u) {
-      p += i.width(); ++i;
-    }
-    // There is at least one range left (check of max)
-    assert(i());
-    // Initialize current range and adjust value
-    c = i.min();
-    // Skip all values which are too small
-    while (d + static_cast<Val>(a)*c < u) {
-      p++; c++;
-    }
-    // Adjust to new value
-    d += static_cast<Val>(a) * c;
-    return true;
-  }
-  template <class Val>
-  forceinline bool
-  PosSupportIter<Val>::adjust(Val& d) {
-    // Current value
-    Val v = static_cast<Val>(a) * c;
-    // Subtract current value from d
-    d -= v;
-    // Move to next position (number of value)
-    p++;
-    // Still in the same range
-    if (c < i.max()) {
-      // Decrement current values
-      c += 1; v += a;
-    } else {
-      // Go to next range
-      ++i;
-      if (!i())
-	return false;
-      c = i.min(); v = static_cast<Val>(a) * c;
-    }
-    // Is d with the current value too large?
-    if (d + v > l)
-      return false;
-    // Update d
-    d += v;
-    return true;
-  }
-
-
-  /*
-   * Support-based iterator for negative view
-   *
-   */
-  template <class Val>
-  forceinline bool
-  NegSupportIter<Val>::reset(Val& d) {
-    // Way too small, no value can make it big enough
-    if (d + static_cast<Val>(a)*x.min() < u)
-      return false;
-    // Restart iterator and position of values
-    i.init(x.variable()); p = x.size()-1;
-    // Skip all ranges which are too small
-    while (d + static_cast<Val>(a)*i.min() < u) {
-      p -= i.width(); ++i;
-    }
-    // There is at least one range left (check of max)
-    assert(i());
-    // Initialize current range
-    c = i.max();
-    // Skip all values which are too small
-    while (d + static_cast<Val>(a)*c < u) {
-      p--; c--;
-    }
-    // Adjust to new value
-    d += static_cast<Val>(a) * c;
-    return true;
-  }
-  template <class Val>
-  forceinline bool
-  NegSupportIter<Val>::adjust(Val& d) {
-    // Current value
-    Val v = static_cast<Val>(a) * c;
-    // Subtract current value from d
-    d -= v;
-    // Move to next position (number of value)
-    p--;
-    // Still in the same range
-    if (c > i.min()) {
-      // Decrement current values
-      c -= 1; v -= a;
-    } else {
-      // Go to next range
-      ++i;
-      if (!i())
-	return false;
-      c = i.max(); v = static_cast<Val>(a) * c;
-    }
-    // Is d with the current value too large?
-    if (d + v > l)
-      return false;
-    // Update d
-    d += v;
-    return true;
-  }
-
-
-
-  /*
-   * The domain-consisten equality propagator
-   *
-   */
-  template <class Val, class View>
-  forceinline
-  DomEq<Val,View>::DomEq(Space* home,
-			 ViewArray<View >& x, ViewArray<View >& y,
-			 Val c)
-    : Lin<Val,View,View,PC_INT_DOM>(home,x,y,c) {}
-
-  template <class Val, class View>
-  ExecStatus
-  DomEq<Val,View>::post(Space* home,
-			ViewArray<View>& x, ViewArray<View>& y,
-			Val c) {
-    (void) new (home) DomEq<Val,View>(home,x,y,c);
-    return ES_OK;
-  }
-
-  template <class Val, class View>
-  forceinline
-  DomEq<Val,View>::DomEq(Space* home, bool share, DomEq<Val,View>& p)
-    : Lin<Val,View,View,PC_INT_DOM>(home,share,p) {}
-
-  template <class Val, class View>
-  Actor*
-  DomEq<Val,View>::copy(Space* home, bool share) {
-    return new (home) DomEq<Val,View>(home,share,*this);
-  }
-
-  template <class Val, class View>
-  PropCost
-  DomEq<Val,View>::cost(void) const {
-    return (View::pme(this) != ME_INT_DOM)
-      ? cost_hi(x.size()+y.size(),PC_LINEAR_LO)
-      : cost_hi(x.size()+y.size(),PC_CRAZY_HI);
-  }
-
-  template <class Val, class View>
-  ExecStatus
-  DomEq<Val,View>::propagate(Space* home) {
-    if (View::pme(this) != ME_INT_DOM) {
-      ExecStatus es = prop_bnd<Val,View,View>(this,home,x,y,c);
-      if ((es == ES_SUBSUMED) || (es == ES_FAILED))
-	return es;
-      return ES_FIX_PARTIAL(View::pme(ME_INT_DOM));
-    }
-
-    // Value of equation for partial assignment
-    Val d = -c;
-
-    int n = x.size();
-    int m = y.size();
-
-    // Create support-base iterators
-    GECODE_AUTOARRAY(PosSupportIter<Val>, xp, n);
-    GECODE_AUTOARRAY(NegSupportIter<Val>, yp, m);
-
-    // Initialize views for assignments
-    {
-      Val l = 0;
-      Val u = 0;
-      for (int j=m; j--; ) {
-	yp[j].init(-y[j].scale(),y[j].base(),l,u);
-	l += y[j].max(); u += y[j].min();
-      }
-      for (int i=n; i--; ) {
-	xp[i].init(x[i].scale(),x[i].base(),l,u);
-	l -= x[i].min(); u -= x[i].max();
-      }
-    }
-
-    // Collect support information by iterating assignments
-    {
-      // Force reset of all iterators in first round
-      int i = 0;
-      int j = 0;
-
-    next_i:
-      // Reset all iterators for positive views and update d
-      while (i<n) {
-	if (!xp[i].reset(d)) goto prev_i;
-	i++;
-      }
-    next_j:
-      // Reset all iterators for negative views and update d
-      while (j<m) {
-	if (!yp[j].reset(d)) goto prev_j;
-	j++;
-      }
-      // Check whether current assignment is solution
-      if (d == 0) {
-	// Record support
-	for (int is=n; is--; ) xp[is].support();
-	for (int js=m; js--; ) yp[js].support();
-      }
-    prev_j:
-      // Try iterating to next assignment: negative views
-      while (j>0) {
-	if (yp[j-1].adjust(d)) goto next_j;
-	j--;
-      }
-    prev_i:
-      // Try iterating to next assignment: positive views
-      while (i>0) {
-	if (xp[i-1].adjust(d)) goto next_i;
-	i--;
-      }
-    }
-
-    // Tell back new variable domains
-    ExecStatus es = ES_SUBSUMED;
-    for (int i=n; i--; ) {
-      ModEvent me = xp[i].tell(home);
-      if (me_failed(me)) {
-	es = ES_FAILED; goto dispose;
-      }
-      if (!x[i].assigned())
-	es = ES_FIX;
-    }
-    for (int j=m; j--; ) {
-      ModEvent me = yp[j].tell(home);
-      if (me_failed(me)) {
-	es = ES_FAILED; goto dispose;
-      }
-      if (!y[j].assigned())
-	es = ES_FIX;
-    }
-
-    // Release memory
-  dispose:
-    for (int i=n; i--; ) xp[i].dispose();
-    for (int j=m; j--; ) yp[j].dispose();
-    return es;
-  }
-
-}}}
-
-// STATISTICS: int-prop