this repo has no description
dekker.one
1/* -*- mode: C++; c-basic-offset: 2; indent-tabs-mode: nil -*- */
2/*
3 * Main authors:
4 * Christian Schulte <schulte@gecode.org>
5 * Guido Tack <tack@gecode.org>
6 * Vincent Barichard <Vincent.Barichard@univ-angers.fr>
7 *
8 * Copyright:
9 * Christian Schulte, 2004
10 * Guido Tack, 2006
11 * Vincent Barichard, 2012
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
38namespace Gecode { namespace Float { namespace Arithmetic {
39
40 template<class A, class B>
41 forceinline
42 Abs<A,B>::Abs(Home home, A x0, B x1)
43 : MixBinaryPropagator<A,PC_FLOAT_BND,B,PC_FLOAT_BND>(home,x0,x1) {}
44
45 template<class A, class B>
46 ExecStatus
47 Abs<A,B>::post(Home home, A x0, B x1) {
48 GECODE_ME_CHECK(x1.eq(home,abs(x0.val())));
49 if (x0.min() >= 0)
50 GECODE_ME_CHECK(x0.eq(home,FloatVal(x1.min(), x1.max())));
51 else if (x0.max() <= 0)
52 GECODE_ME_CHECK(x0.eq(home,FloatVal(-x1.max(), -x1.min())));
53 else
54 GECODE_ME_CHECK(x0.eq(home,FloatVal(-x1.max(), x1.max())));
55 (void) new (home) Abs<A,B>(home,x0,x1);
56 return ES_OK;
57 }
58
59 template<class A, class B>
60 forceinline
61 Abs<A,B>::Abs(Space& home, Abs<A,B>& p)
62 : MixBinaryPropagator<A,PC_FLOAT_BND,B,PC_FLOAT_BND>(home,p) {}
63
64 template<class A, class B>
65 Actor*
66 Abs<A,B>::copy(Space& home) {
67 return new (home) Abs<A,B>(home,*this);
68 }
69
70 template<class A, class B>
71 ExecStatus
72 Abs<A,B>::propagate(Space& home, const ModEventDelta&) {
73 GECODE_ME_CHECK(x1.eq(home,abs(x0.val())));
74 if (x0.min() >= 0)
75 GECODE_ME_CHECK(x0.eq(home,FloatVal(x1.min(), x1.max())));
76 else if (x0.max() <= 0)
77 GECODE_ME_CHECK(x0.eq(home,FloatVal(-x1.max(), -x1.min())));
78 else
79 GECODE_ME_CHECK(x0.eq(home,FloatVal(-x1.max(), x1.max())));
80 return (x0.assigned() && x1.assigned()) ? home.ES_SUBSUMED(*this) : ES_FIX;
81 }
82
83}}}
84
85// STATISTICS: float-prop
86