main page
modules
namespaces
classes
files
Gecode home
Generated on Tue Oct 22 2013 00:48:58 for Gecode by
doxygen
1.8.4
examples
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
* Last modified:
10
* $Date: 2013-07-08 14:22:40 +0200 (Mon, 08 Jul 2013) $ by $Author: schulte $
11
* $Revision: 13820 $
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
72
class
DescartesFolium
:
public
Script
{
73
protected
:
75
FloatVarArray
f
;
77
double
step
;
78
public
:
80
DescartesFolium
(
const
Options
&)
81
: f(*this,3,-20,20),
step
(0.1) {
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
}
94
DescartesFolium
(
bool
share,
DescartesFolium
&
p
)
95
:
Script
(share,p),
step
(p.
step
) {
96
f.update(*
this
,share,p.
f
);
97
}
99
virtual
Space
*
copy
(
bool
share) {
100
return
new
DescartesFolium
(share,*
this
);
101
}
103
virtual
void
constrain
(
const
Space
&
_b
) {
104
const
DescartesFolium
&
b
=
static_cast<
const
DescartesFolium
&
>
(
_b
);
105
rel
(*
this
, f[0] >= (b.
f
[0].max()+
step
));
106
}
108
virtual
void
print
(std::ostream& os)
const
{
109
os <<
"XY "
<< f[1].med() <<
" "
<< f[2].med()
110
<< std::endl;
111
}
112
113
};
114
118
int
main
(
int
argc,
char
* argv[]) {
119
Options
opt
(
"DescartesFolium"
);
120
opt.
parse
(argc,argv);
121
opt.
solutions
(0);
122
Script::run<DescartesFolium,BAB,Options>(
opt
);
123
return
0;
124
}
125
126
// STATISTICS: example-any