Generated on Fri Oct 19 11:24:48 2018 for Gecode by doxygen 1.6.3

descartes-folium.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  *  This file is part of Gecode, the generic constraint
00010  *  development environment:
00011  *     http://www.gecode.org
00012  *
00013  *  Permission is hereby granted, free of charge, to any person obtaining
00014  *  a copy of this software and associated documentation files (the
00015  *  "Software"), to deal in the Software without restriction, including
00016  *  without limitation the rights to use, copy, modify, merge, publish,
00017  *  distribute, sublicense, and/or sell copies of the Software, and to
00018  *  permit persons to whom the Software is furnished to do so, subject to
00019  *  the following conditions:
00020  *
00021  *  The above copyright notice and this permission notice shall be
00022  *  included in all copies or substantial portions of the Software.
00023  *
00024  *  THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
00025  *  EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
00026  *  MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
00027  *  NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE
00028  *  LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
00029  *  OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
00030  *  WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
00031  *
00032  */
00033 
00034 #include <gecode/driver.hh>
00035 
00036 #include <gecode/minimodel.hh>
00037 #include <gecode/float.hh>
00038 
00039 using namespace Gecode;
00040 
00068 class DescartesFolium : public FloatMaximizeScript {
00069 protected:
00071   FloatVarArray f;
00073   double step;
00074 public:
00076   DescartesFolium(const Options& opt)
00077     : FloatMaximizeScript(opt), f(*this,3,-20,20) {
00078 
00079     if (opt.trace() != 0)
00080       trace(*this, f, opt.trace());
00081 
00082     // Post equation
00083     FloatVar p = f[0];
00084     FloatVar x = f[1];
00085     FloatVar y = f[2];
00086     rel(*this, 3*p/(1+pow(p,3)) == x);
00087     rel(*this, 3*sqr(p)/(1+pow(p,3)) == y);
00088     rel(*this, pow(x,3) + pow(y,3)  == 3 * x * y);
00089     rel(*this, x == FloatVal(-1,2));
00090 
00091     branch(*this,p,FLOAT_VAL_SPLIT_MIN());
00092   }
00094   DescartesFolium(DescartesFolium& p)
00095     : FloatMaximizeScript(p) {
00096     f.update(*this, p.f);
00097   }
00099   virtual Space* copy(void) {
00100     return new DescartesFolium(*this);
00101   }
00103   virtual FloatVar cost(void) const {
00104     return f[0];
00105   }
00107   virtual void print(std::ostream& os) const {
00108     os << "XY " << f[1].med() << " " << f[2].med()
00109        << std::endl;
00110   }
00111 
00112 };
00113 
00117 int main(int argc, char* argv[]) {
00118   Options opt("DescartesFolium");
00119   opt.solutions(0);
00120   opt.step(0.1);
00121   opt.parse(argc,argv);
00122   FloatMaximizeScript::run<DescartesFolium,BAB,Options>(opt);
00123   return 0;
00124 }
00125 
00126 // STATISTICS: example-any