Generated on Sat May 25 2013 18:00:31 for Gecode by doxygen 1.8.3.1
archimedean-spiral.cpp
Go to the documentation of this file.
1 /* -*- mode: C++; c-basic-offset: 2; indent-tabs-mode: nil -*- */
2 /*
3  * Main authors:
4  * Vincent Barichard <Vincent.Barichard@univ-angers.fr>
5  *
6  * Copyright:
7  * Vincent Barichard, 2012
8  *
9  * Last modified:
10  * $Date: 2013-02-04 22:40:26 +0100 (Mon, 04 Feb 2013) $ by $Author: schulte $
11  * $Revision: 13263 $
12  *
13  * This file is part of Gecode, the generic constraint
14  * development environment:
15  * http://www.gecode.org
16  *
17  * Permission is hereby granted, free of charge, to any person obtaining
18  * a copy of this software and associated documentation files (the
19  * "Software"), to deal in the Software without restriction, including
20  * without limitation the rights to use, copy, modify, merge, publish,
21  * distribute, sublicense, and/or sell copies of the Software, and to
22  * permit persons to whom the Software is furnished to do so, subject to
23  * the following conditions:
24  *
25  * The above copyright notice and this permission notice shall be
26  * included in all copies or substantial portions of the Software.
27  *
28  * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
29  * EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
30  * MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
31  * NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE
32  * LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
33  * OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
34  * WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
35  *
36  */
37 
38 #include <gecode/driver.hh>
39 
40 #include <gecode/minimodel.hh>
41 #include <gecode/float.hh>
42 
43 using namespace Gecode;
44 
71 class ArchimedeanSpiral : public Script {
72 protected:
77 public:
80  : f(*this,4,-20,20), step(0.1) {
81  // Post equation
82  FloatVar theta = f[0];
83  FloatVar r = f[3];
84  FloatVar x = f[1];
85  FloatVar y = f[2];
86 
87  rel(*this, theta >= 0);
88  rel(*this, theta <= 6*FloatVal::pi());
89  rel(*this, r >= 0);
90  rel(*this, r*cos(theta) == x);
91  rel(*this, r*sin(theta) == y);
92  rel(*this, r == theta);
93 
94  branch(*this,f[0],FLOAT_VAL_SPLIT_MIN());
95  }
98  : Script(share,p), step(p.step) {
99  f.update(*this,share,p.f);
100  }
102  virtual Space* copy(bool share) {
103  return new ArchimedeanSpiral(share,*this);
104  }
106  virtual void constrain(const Space& _b) {
107  const ArchimedeanSpiral& b = static_cast<const ArchimedeanSpiral&>(_b);
108  rel(*this, f[0] >= (b.f[0].max()+step));
109  }
111  virtual void print(std::ostream& os) const {
112  os << "XY " << f[1].med() << " " << f[2].med()
113  << std::endl;
114  }
115 
116 };
117 
121 int main(int argc, char* argv[]) {
122  Options opt("ArchimedeanSpiral");
123  opt.parse(argc,argv);
124  opt.solutions(0);
125  Script::run<ArchimedeanSpiral,BAB,Options>(opt);
126  return 0;
127 }
128 
129 // STATISTICS: example-any