main page
modules
namespaces
classes
files
Gecode home
Generated on Tue Oct 22 2013 00:49:01 for Gecode by
doxygen
1.8.4
gecode
float
transcendental
exp-log.hpp
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: 2012-07-27 10:25:47 +0200 (Fri, 27 Jul 2012) $ by $Author: vbarichard $
11
* $Revision: 12971 $
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
namespace
Gecode {
namespace
Float {
namespace
Transcendental {
39
40
/*
41
* Bounds consistent exponential operator
42
*
43
*/
44
45
template
<
class
A,
class
B>
46
forceinline
47
Exp<A,B>::Exp
(
Home
home, A x0, B x1)
48
:
MixBinaryPropagator
<A,
PC_FLOAT_BND
,B,
PC_FLOAT_BND
>(home,x0,x1) {}
49
50
template
<
class
A,
class
B>
51
ExecStatus
52
Exp<A,B>::post
(
Home
home, A x0, B x1) {
53
if
(
same
(x0,x1)) {
54
return
ES_FAILED
;
55
}
else
{
56
GECODE_ME_CHECK
(x1.gq(home,0.0));
57
}
58
59
(void)
new
(home)
Exp<A,B>
(home,x0,x1);
60
return
ES_OK
;
61
}
62
63
64
template
<
class
A,
class
B>
65
forceinline
66
Exp<A,B>::Exp
(
Space
& home,
bool
share,
Exp<A,B>
&
p
)
67
:
MixBinaryPropagator
<A,
PC_FLOAT_BND
,B,
PC_FLOAT_BND
>(home,share,p) {}
68
69
template
<
class
A,
class
B>
70
Actor
*
71
Exp<A,B>::copy
(
Space
& home,
bool
share) {
72
return
new
(home)
Exp<A,B>
(home,share,*
this
);
73
}
74
75
template
<
class
A,
class
B>
76
ExecStatus
77
Exp<A,B>::propagate
(
Space
& home,
const
ModEventDelta
&) {
78
GECODE_ME_CHECK
(x1.eq(home,
exp
(x0.domain())));
79
if
(x1.max() == 0.0)
80
return
ES_FAILED
;
81
GECODE_ME_CHECK
(x0.eq(home,
log
(x1.domain())));
82
return
x0.assigned() ? home.
ES_SUBSUMED
(*
this
) :
ES_FIX
;
83
}
84
85
86
/*
87
* Bounds consistent logarithm operator with base
88
*
89
*/
90
91
template
<
class
A,
class
B>
92
forceinline
93
Pow<A,B>::Pow
(
Home
home,
FloatNum
base0, A x0, B x1)
94
:
MixBinaryPropagator
<A,
PC_FLOAT_BND
,B,
PC_FLOAT_BND
>(home,x0,x1),
95
base(base0) {}
96
97
template
<
class
A,
class
B>
98
ExecStatus
99
Pow<A,B>::post
(
Home
home,
FloatNum
base, A x0, B x1) {
100
if
(base <= 0)
return
ES_FAILED
;
101
if
(
same
(x0,x1)) {
102
GECODE_ME_CHECK
(x0.eq(home,0.0));
103
}
else
{
104
GECODE_ME_CHECK
(x1.gq(home,0.0));
105
(void)
new
(home)
Pow<A,B>
(home,base,x0,x1);
106
}
107
return
ES_OK
;
108
}
109
110
template
<
class
A,
class
B>
111
forceinline
112
Pow<A,B>::Pow
(
Space
& home,
bool
share,
Pow<A,B>
&
p
)
113
:
MixBinaryPropagator
<A,
PC_FLOAT_BND
,B,
PC_FLOAT_BND
>(home,share,p),
114
base(p.base) {}
115
116
template
<
class
A,
class
B>
117
Actor
*
118
Pow<A,B>::copy
(
Space
& home,
bool
share) {
119
return
new
(home)
Pow<A,B>
(home,share,*
this
);
120
}
121
122
template
<
class
A,
class
B>
123
ExecStatus
124
Pow<A,B>::propagate
(
Space
& home,
const
ModEventDelta
&) {
125
if
(x1.max() == 0.0)
126
return
ES_FAILED
;
127
GECODE_ME_CHECK
(x0.eq(home,
log
(x1.domain())/
log
(base)));
128
GECODE_ME_CHECK
(x1.eq(home,
exp
(x0.domain()*
log
(base))));
129
return
x0.assigned() ? home.
ES_SUBSUMED
(*
this
) :
ES_FIX
;
130
}
131
132
}}}
133
134
// STATISTICS: float-prop
135