Generated on Tue Feb 11 2025 17:33:26 for Gecode by doxygen 1.12.0
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
39using namespace Gecode;
40
69protected:
73 double step;
74public:
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
92 }
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
117int main(int argc, char* argv[]) {
118 Options opt("DescartesFolium");
119 opt.solutions(0);
120 opt.step(0.1);
121 opt.parse(argc,argv);
123 return 0;
124}
125
126// STATISTICS: example-any
int p
Number of positive literals for node type.
Example: Folium of Descartes
int main(int argc, char *argv[])
Main-function.
DescartesFolium(DescartesFolium &p)
Constructor for cloning p.
FloatVarArray f
The numbers.
DescartesFolium(const Options &opt)
Actual model.
virtual FloatVar cost(void) const
Cost function.
double step
Minimum distance between two solutions.
virtual Space * copy(void)
Copy during cloning.
virtual void print(std::ostream &os) const
Print solution coordinates.
Parametric base-class for scripts.
Definition driver.hh:729
static void run(const Options &opt, Script *s=NULL)
Definition script.hpp:290
Float value type.
Definition float.hh:334
Float variable array.
Definition float.hh:1030
Float variables.
Definition float.hh:870
Options for scripts
Definition driver.hh:366
Computation spaces.
Definition core.hpp:1742
void update(Space &home, VarArray< Var > &a)
Update array to be a clone of array a.
Definition array.hpp:1013
void parse(int argc, char *argv[])
Parse commandline arguments.
Definition test.cpp:120
GECODE_FLOAT_EXPORT void trace(Home home, const FloatVarArgs &x, TraceFilter tf, int te=(TE_INIT|TE_PRUNE|TE_FIX|TE_FAIL|TE_DONE), FloatTracer &t=StdFloatTracer::def)
Create a tracer for float variables.
Definition trace.cpp:39
void branch(Home home, const FloatVarArgs &x, FloatVarBranch vars, FloatValBranch vals, FloatBranchFilter bf=nullptr, FloatVarValPrint vvp=nullptr)
Branch over x with variable selection vars and value selection vals.
Definition branch.cpp:39
void rel(Home home, FloatVar x0, FloatRelType frt, FloatVar x1)
Post propagator for .
Definition rel.cpp:68
Gecode toplevel namespace
void sqr(Home home, FloatVar x0, FloatVar x1)
Post propagator for .
Post propagator for SetVar SetOpType SetVar y
Definition set.hh:767
void pow(Home home, FloatVar x0, int n, FloatVar x1)
Post propagator for for $n\geq 0$.
Post propagator for SetVar x
Definition set.hh:767
FloatValBranch FLOAT_VAL_SPLIT_MIN(void)
Select values not greater than mean of smallest and largest value.
Definition val.hpp:55