Generated on Fri Jan 10 2020 11:38:25 for Gecode by doxygen 1.8.16
descartes-folium.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  * This file is part of Gecode, the generic constraint
10  * development environment:
11  * http://www.gecode.org
12  *
13  * Permission is hereby granted, free of charge, to any person obtaining
14  * a copy of this software and associated documentation files (the
15  * "Software"), to deal in the Software without restriction, including
16  * without limitation the rights to use, copy, modify, merge, publish,
17  * distribute, sublicense, and/or sell copies of the Software, and to
18  * permit persons to whom the Software is furnished to do so, subject to
19  * the following conditions:
20  *
21  * The above copyright notice and this permission notice shall be
22  * included in all copies or substantial portions of the Software.
23  *
24  * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
25  * EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
26  * MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
27  * NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE
28  * LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
29  * OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
30  * WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
31  *
32  */
33 
34 #include <gecode/driver.hh>
35 
36 #include <gecode/minimodel.hh>
37 #include <gecode/float.hh>
38 
39 using namespace Gecode;
40 
69 protected:
73  double step;
74 public:
77  : FloatMaximizeScript(opt), f(*this,3,-20,20) {
78 
79  if (opt.trace() != 0)
80  trace(*this, f, opt.trace());
81 
82  // Post equation
83  FloatVar p = f[0];
84  FloatVar x = f[1];
85  FloatVar y = f[2];
86  rel(*this, 3*p/(1+pow(p,3)) == x);
87  rel(*this, 3*sqr(p)/(1+pow(p,3)) == y);
88  rel(*this, pow(x,3) + pow(y,3) == 3 * x * y);
89  rel(*this, x == FloatVal(-1,2));
90 
91  branch(*this,p,FLOAT_VAL_SPLIT_MIN());
92  }
96  f.update(*this, p.f);
97  }
99  virtual Space* copy(void) {
100  return new DescartesFolium(*this);
101  }
103  virtual FloatVar cost(void) const {
104  return f[0];
105  }
107  virtual void print(std::ostream& os) const {
108  os << "XY " << f[1].med() << " " << f[2].med()
109  << std::endl;
110  }
111 
112 };
113 
117 int main(int argc, char* argv[]) {
118  Options opt("DescartesFolium");
119  opt.solutions(0);
120  opt.step(0.1);
121  opt.parse(argc,argv);
122  FloatMaximizeScript::run<DescartesFolium,BAB,Options>(opt);
123  return 0;
124 }
125 
126 // STATISTICS: example-any
void pow(Home home, FloatVar x0, int n, FloatVar x1)
Post propagator for for $n\geq 0$.
Definition: arithmetic.cpp:109
Post propagator for SetVar x
Definition: set.hh:767
Post propagator for SetVar SetOpType SetVar y
Definition: set.hh:767
DescartesFolium(DescartesFolium &p)
Constructor for cloning p.
int main(int argc, char *argv[])
Main-function.
void sqr(Home home, FloatVar x0, FloatVar x1)
Post propagator for .
Definition: arithmetic.cpp:95
virtual Space * copy(void)
Copy during cloning.
Computation spaces.
Definition: core.hpp:1742
void branch(Home home, const FloatVarArgs &x, FloatVarBranch vars, FloatValBranch vals, FloatBranchFilter bf, FloatVarValPrint vvp)
Branch over x with variable selection vars and value selection vals.
Definition: branch.cpp:39
Gecode toplevel namespace
Options opt
The options.
Definition: test.cpp:97
DescartesFolium(const Options &opt)
Actual model.
Options for scripts
Definition: driver.hh:366
FloatValBranch FLOAT_VAL_SPLIT_MIN(void)
Select values not greater than mean of smallest and largest value.
Definition: val.hpp:55
Parametric base-class for scripts.
Definition: driver.hh:729
void step(double s)
Set default step value.
Definition: options.hpp:256
FloatVarArray f
The numbers.
void parse(int &argc, char *argv[])
Parse options from arguments argv (number is argc)
Definition: options.cpp:548
virtual void print(std::ostream &os) const
Print solution coordinates.
virtual FloatVar cost(void) const
Cost function.
Float value type.
Definition: float.hh:334
void rel(Home home, FloatVar x0, FloatRelType frt, FloatVal n)
Propagates .
Definition: rel.cpp:43
Float variables.
Definition: float.hh:870
Float variable array.
Definition: float.hh:1030
void trace(Home home, const FloatVarArgs &x, TraceFilter tf, int te, FloatTracer &t)
Create a tracer for float variables.
Definition: trace.cpp:39
Example: Folium of Descartes
void solutions(unsigned int n)
Set default number of solutions to search for.
Definition: options.hpp:283
double step
Minimum distance between two solutions.
Post propagator for f(x \diamond_{\mathit{op}} y) \sim_r z \f$ void rel(Home home
int p
Number of positive literals for node type.
Definition: bool-expr.cpp:232