Generated on Wed Nov 5 2014 05:18:15 for Gecode by doxygen 1.7.6.1
archimedean-spiral.cpp
Go to the documentation of this file.
00001 /* -*- mode: C++; c-basic-offset: 2; indent-tabs-mode: nil -*- */
00002 /*
00003  *  Main authors:
00004  *     Vincent Barichard <Vincent.Barichard@univ-angers.fr>
00005  *
00006  *  Copyright:
00007  *     Vincent Barichard, 2012
00008  *
00009  *  Last modified:
00010  *     $Date: 2013-07-08 14:22:40 +0200 (Mon, 08 Jul 2013) $ by $Author: schulte $
00011  *     $Revision: 13820 $
00012  *
00013  *  This file is part of Gecode, the generic constraint
00014  *  development environment:
00015  *     http://www.gecode.org
00016  *
00017  *  Permission is hereby granted, free of charge, to any person obtaining
00018  *  a copy of this software and associated documentation files (the
00019  *  "Software"), to deal in the Software without restriction, including
00020  *  without limitation the rights to use, copy, modify, merge, publish,
00021  *  distribute, sublicense, and/or sell copies of the Software, and to
00022  *  permit persons to whom the Software is furnished to do so, subject to
00023  *  the following conditions:
00024  *
00025  *  The above copyright notice and this permission notice shall be
00026  *  included in all copies or substantial portions of the Software.
00027  *
00028  *  THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
00029  *  EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
00030  *  MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
00031  *  NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE
00032  *  LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
00033  *  OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
00034  *  WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
00035  *
00036  */
00037 
00038 #include <gecode/driver.hh>
00039 
00040 #include <gecode/minimodel.hh>
00041 #include <gecode/float.hh>
00042 
00043 using namespace Gecode;
00044 
00071 class ArchimedeanSpiral : public Script {
00072 protected:
00074   FloatVarArray f;
00076   FloatNum step;
00077 public:
00079   ArchimedeanSpiral(const Options&) 
00080     : f(*this,4,-20,20), step(0.1) {
00081     // Post equation
00082     FloatVar theta = f[0];
00083     FloatVar r = f[3];
00084     FloatVar x = f[1];
00085     FloatVar y = f[2];
00086 
00087     rel(*this, theta >= 0);
00088     rel(*this, theta <= 6*FloatVal::pi());
00089     rel(*this, r >= 0);
00090     rel(*this, r*cos(theta) == x);
00091     rel(*this, r*sin(theta) == y);
00092     rel(*this, r == theta);
00093 
00094     branch(*this,f[0],FLOAT_VAL_SPLIT_MIN());
00095   }
00097   ArchimedeanSpiral(bool share, ArchimedeanSpiral& p) 
00098     : Script(share,p), step(p.step) {
00099     f.update(*this,share,p.f);
00100   }
00102   virtual Space* copy(bool share) { 
00103     return new ArchimedeanSpiral(share,*this); 
00104   }
00106   virtual void constrain(const Space& _b) {
00107     const ArchimedeanSpiral& b = static_cast<const ArchimedeanSpiral&>(_b);
00108     rel(*this, f[0] >= (b.f[0].max()+step));
00109   }
00111   virtual void print(std::ostream& os) const {
00112     os << "XY " << f[1].med() << " " << f[2].med()
00113        << std::endl;
00114   }
00115   
00116 };
00117 
00121 int main(int argc, char* argv[]) {
00122   Options opt("ArchimedeanSpiral");
00123   opt.parse(argc,argv);
00124   opt.solutions(0);
00125   Script::run<ArchimedeanSpiral,BAB,Options>(opt);
00126   return 0;
00127 }
00128 
00129 // STATISTICS: example-any