gecoder-with-gecode 0.7.1

data/ext/gecode-1.3.1/examples/steiner.cc

@@ -0,0 +1,148 @@
+/*
+ *  Main authors:
+ *     Guido Tack <tack@gecode.org>
+ *     Christian Schulte <schulte@gecode.org>
+ *
+ *  Copyright:
+ *     Guido Tack, 2004
+ *     Christian Schulte, 2004
+ *
+ *  Last modified:
+ *     $Date: 2006-08-04 16:06:52 +0200 (Fri, 04 Aug 2006) $ by $Author: schulte $
+ *     $Revision: 3517 $
+ *
+ *  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/set.hh"
+#include "examples/support.hh"
+
+/**
+ * \brief %Example: Steiner triples
+ *
+ * This is problem 044 from csplib.
+ *
+ * \ingroup Example
+ *
+ */
+class Steiner : public Example {
+public:
+
+  int n;
+  int n1;
+  int n1n1;
+  int len;
+
+  SetVarArray root;
+
+  Steiner(const Options& o)
+    : n(o.size),
+      n1(n+1), n1n1(n1*n1), len((n*(n-1))/6),
+      root(this,len, IntSet::empty, 1, n, 3, 3) {
+
+    for (int i=0; i<len; i++) {
+      for (int j=i+1; j<len; j++) {
+        SetVar x = root[i];
+        SetVar y = root[j];
+
+	SetVar atmostOne(this,IntSet::empty,1,n,0,1);
+	cardinality(this, atmostOne,0,1);
+	rel(this, x, SOT_INTER, y, SRT_EQ, atmostOne);
+
+	IntVar x1(this,1,n);
+	IntVar x2(this,1,n);
+	IntVar x3(this,1,n);
+	IntVar y1(this,1,n);
+	IntVar y2(this,1,n);
+	IntVar y3(this,1,n);
+
+	if (o.naive) {
+	
+	  /* Naive alternative:
+	   * just including the ints in the set
+	   */
+	  rel(this, x, SRT_SUP, x1);
+	  rel(this, x, SRT_SUP, x2);
+	  rel(this, x, SRT_SUP, x3);
+	  rel(this, y, SRT_SUP, y1);
+	  rel(this, y, SRT_SUP, y2);
+	  rel(this, y, SRT_SUP, y3);
+	} else {
+
+	  /* Smart alternative:
+	   * Using matching constraints
+	   */
+	
+	  IntVarArgs xargs(3);
+	  xargs[0] = x1; xargs[1] = x2; xargs[2] = x3;
+	  match(this, x,xargs);
+	
+	  IntVarArgs yargs(3);
+	  yargs[0] = y1; yargs[1] = y2; yargs[2] = y3;
+	  match(this, y,yargs);
+	
+	}
+	
+	/* Breaking symmetries */
+	
+	rel(this, x1,IRT_LE,x2);
+	rel(this, x2,IRT_LE,x3);
+	rel(this, x1,IRT_LE,x3);
+	
+	rel(this, y1,IRT_LE,y2);
+	rel(this, y2,IRT_LE,y3);
+	rel(this, y1,IRT_LE,y3);
+	
+	IntArgs ia(6,n1n1,n1,1,-n1n1,-n1,-1);
+	IntVarArgs iva(6);
+	iva[0]=x1; iva[1]=x2; iva[2]=x3;
+	iva[3]=y1; iva[4]=y2; iva[5]=y3;
+	linear(this, ia, iva, IRT_LE, 0);
+      }
+    }
+
+
+    branch(this, root, SETBVAR_NONE, SETBVAL_MIN);
+  }
+
+  Steiner(bool share, Steiner& s) : Example(share,s),
+				    n(s.n), n1(s.n1), n1n1(s.n1n1), len(s.len)
+  {
+    root.update(this, share, s.root);
+  }
+
+  virtual Space*
+  copy(bool share) {
+    return new Steiner(share,*this);
+  }
+
+  virtual void
+  print(void) {
+    for (int i=0; i<len; i++) {
+      std::cout << "\t[" << i << "]" << root[i] << std::endl;
+    }
+  }
+
+};
+
+int
+main(int argc, char** argv) {
+  Options o("Steiner");
+  o.size = 9;
+  o.naive = true;
+  o.iterations = 20;
+  o.parse(argc,argv);
+  Example::run<Steiner,DFS>(o);
+  return 0;
+}
+
+
+// STATISTICS: example-any
+

data/ext/gecode-1.3.1/examples/stress-domain.cc

@@ -0,0 +1,86 @@
+/*
+ *  Main authors:
+ *     Christian Schulte <schulte@gecode.org>
+ *
+ *  Copyright:
+ *     Christian Schulte, 2001
+ *
+ *  Last modified:
+ *     $Date: 2006-08-03 13:51:17 +0200 (Thu, 03 Aug 2006) $ by $Author: schulte $
+ *     $Revision: 3506 $
+ *
+ *  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 "examples/support.hh"
+
+/**
+ * \brief %Example: Domain stress test
+ *
+ * \ingroup Example
+ *
+ */
+class StressDomain : public Example {
+protected:
+  /// Variables
+  IntVarArray x;
+public:
+  /// The actual problem
+  StressDomain(const Options& opt)
+    : x(this,5,0,5*opt.size) {
+
+    // Cut holes: expand
+    for (int i = 5; i--; ) {
+      for (unsigned int j = 0; j <= 5*opt.size; j++)
+	rel(this, x[i], IRT_NQ, 5*j);
+      for (unsigned int j = 0; j <= 5*opt.size; j++)
+	rel(this, x[i], IRT_NQ, 5*j+2);
+      for (unsigned int j = 0; j <= 5*opt.size; j++)
+	rel(this, x[i], IRT_NQ, 5*j+4);
+    }
+    // Contract
+    for (unsigned int j = 0; j <= 5*opt.size/2; j++)
+      for (unsigned int i = 5; i--; ) {
+	rel(this, x[i], IRT_GQ, 5*j);
+	rel(this, x[i], IRT_LQ, 5*(j + (5*opt.size/2)));
+      }
+  }
+
+  /// Constructor for cloning \a s
+  StressDomain(bool share, StressDomain& s) : Example(share,s) {
+    x.update(this, share, s.x);
+  }
+
+  /// Perform copying during cloning
+  virtual Space*
+  copy(bool share) {
+    return new StressDomain(share,*this);
+  }
+
+  /// Print solution
+  virtual void
+  print(void) {}
+};
+
+/** \brief Main-function
+ *  \relates StressDomain
+ */
+int
+main(int argc, char** argv) {
+  Options opt("StressDomain");
+  opt.iterations = 200;
+  opt.size       = 1000;
+  opt.parse(argc,argv);
+  Example::run<StressDomain,DFS>(opt);
+  return 0;
+}
+
+// STATISTICS: example-any
+

data/ext/gecode-1.3.1/examples/stress-element.cc

@@ -0,0 +1,95 @@
+/*
+ *  Main authors:
+ *     Christian Schulte <schulte@gecode.org>
+ *
+ *  Copyright:
+ *     Christian Schulte, 2006
+ *
+ *  Last modified:
+ *     $Date: 2006-08-04 16:06:52 +0200 (Fri, 04 Aug 2006) $ by $Author: schulte $
+ *     $Revision: 3517 $
+ *
+ *  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 "examples/support.hh"
+
+static const int n = 15;
+static const int m = 90;
+static const int p[n] = {16,35,90,42,88,6,40,42,64,48,46,5,90,29,70};
+
+/**
+ * \brief %Example: Stress test for element constraint (involving integers)
+ *
+ * This stress test is an adaption of the benchmark for element posted
+ * by Neng-Fa Zhou for B-Prolog, November, 2005.
+ *
+ * \ingroup Example
+ *
+ */
+class StressElement : public Example {
+protected:
+  /// Variables
+  IntVarArray x;
+public:
+  /// The actual problem
+  StressElement(const Options& opt)
+    : x(this,n,0,n-1) {
+
+    IntVarArgs s(n);
+    for (int i=0; i<n; i++) {
+      IntVar si(this,0,m); s[i]=si;
+    }
+
+    for (int i=0; i<n-1; i++)
+      rel(this, s[i], IRT_LQ, s[i+1]);
+
+    IntArgs e(n,p);
+
+    for (int i=0; i<n; i++)
+      element(this, e, x[i], s[i]);
+
+    distinct(this, x);
+    branch(this, x, BVAR_SIZE_MIN, BVAL_MIN);
+  }
+
+  /// Constructor for cloning \a s
+  StressElement(bool share, StressElement& s) : Example(share,s) {
+    x.update(this, share, s.x);
+  }
+
+  /// Perform copying during cloning
+  virtual Space*
+  copy(bool share) {
+    return new StressElement(share,*this);
+  }
+
+  /// Print solution
+  virtual void
+  print(void) {
+    std::cout << "\tx[" << n << "] = {";
+    for (int i = 0; i < n; i++)
+      std::cout << x[i] << ((i<n-1)?",":"};\n");
+  }
+};
+
+/** \brief Main-function
+ *  \relates StressElement
+ */
+int
+main(int argc, char** argv) {
+  Options opt("StressElement");
+  opt.parse(argc,argv);
+  Example::run<StressElement,DFS>(opt);
+  return 0;
+}
+
+// STATISTICS: example-any
+

data/ext/gecode-1.3.1/examples/stress-exec.cc

@@ -0,0 +1,74 @@
+/*
+ *  Main authors:
+ *     Christian Schulte <schulte@gecode.org>
+ *
+ *  Copyright:
+ *     Christian Schulte, 2001
+ *
+ *  Last modified:
+ *     $Date: 2006-08-03 13:51:17 +0200 (Thu, 03 Aug 2006) $ by $Author: schulte $
+ *     $Revision: 3506 $
+ *
+ *  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 "examples/support.hh"
+
+/**
+ * \brief %Example: Execution stress test
+ *
+ * \ingroup Example
+ *
+ */
+class StressExec : public Example {
+protected:
+  /// Variables
+  IntVarArray x;
+public:
+  /// The actual problem
+  StressExec(const Options& opt)
+    : x(this,2,0,opt.size) {
+
+    rel(this, x[0], IRT_LE, x[1]);
+    rel(this, x[1], IRT_LE, x[0]);
+
+  }
+
+  /// Constructor for cloning \a s
+  StressExec(bool share, StressExec& s) : Example(share,s) {
+    x.update(this, share, s.x);
+  }
+
+  /// Perform copying during cloning
+  virtual Space*
+  copy(bool share) {
+    return new StressExec(share,*this);
+  }
+
+  /// Print solution
+  virtual void
+  print(void) {}
+};
+
+/** \brief Main-function
+ *  \relates StressExec
+ */
+int
+main(int argc, char** argv) {
+  Options opt("StressExec");
+  opt.iterations = 20;
+  opt.size       = 1000000;
+  opt.parse(argc,argv);
+  Example::run<StressExec,DFS>(opt);
+  return 0;
+}
+
+// STATISTICS: example-any
+

data/ext/gecode-1.3.1/examples/stress-min.cc

@@ -0,0 +1,84 @@
+/*
+ *  Main authors:
+ *     Christian Schulte <schulte@gecode.org>
+ *
+ *  Copyright:
+ *     Christian Schulte, 2004
+ *
+ *  Last modified:
+ *     $Date: 2006-07-25 12:42:29 +0200 (Tue, 25 Jul 2006) $ by $Author: schulte $
+ *     $Revision: 3456 $
+ *
+ *  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 "examples/support.hh"
+
+/**
+ * \brief %Example: Stress test for the minimum constraint
+ *
+ * \ingroup Example
+ *
+ */
+class StressMin : public Example {
+protected:
+  /// Size of problem
+  const int n;
+  /// Variables
+  IntVarArray x;
+public:
+  /// The actual problem
+  StressMin(const Options& opt)
+    : n(opt.size), x(this,n,0,2*n-1) {
+    for (int i=1; i<n; i++) {
+      IntVarArgs y(i);
+      for (int j=0; j<i; j++)
+	y[j]=x[j];
+      IntVar m(this,0,2*n);
+      min(this, y, m);
+      rel(this, m, IRT_GR, x[i]);
+    }
+    branch(this, x, BVAR_NONE, BVAL_SPLIT_MAX);
+  }
+
+  /// Constructor for cloning \a s
+  StressMin(bool share, StressMin& s) : Example(share,s), n(s.n) {
+    x.update(this, share, s.x);
+  }
+
+  /// Perform copying during cloning
+  virtual Space*
+  copy(bool share) {
+    return new StressMin(share,*this);
+  }
+
+  /// Print solution
+  virtual void
+  print(void) {
+    std::cout << "\tx[" << n << "] = {";
+    for (int i = 0; i < n; i++)
+      std::cout << x[i] << ((i<n-1)?",":"};\n");
+  }
+};
+
+
+/** \brief Main-function
+ *  \relates StressMin
+ */
+int
+main(int argc, char** argv) {
+  Options opt("StressMin");
+  opt.parse(argc,argv);
+  opt.size = 200;
+  Example::run<StressMin,DFS>(opt);
+  return 0;
+}
+
+// STATISTICS: example-any