this repo has no description
dekker.one
1/* -*- mode: C++; c-basic-offset: 2; indent-tabs-mode: nil -*- */
2/*
3 * Main authors:
4 * Guido Tack <tack@gecode.org>
5 * Christian Schulte <schulte@gecode.org>
6 *
7 * Contributing authors:
8 * Gabor Szokoli <szokoli@gecode.org>
9 *
10 * Copyright:
11 * Guido Tack, 2004
12 * Christian Schulte, 2004
13 * Gabor Szokoli, 2004
14 *
15 * This file is part of Gecode, the generic constraint
16 * development environment:
17 * http://www.gecode.org
18 *
19 * Permission is hereby granted, free of charge, to any person obtaining
20 * a copy of this software and associated documentation files (the
21 * "Software"), to deal in the Software without restriction, including
22 * without limitation the rights to use, copy, modify, merge, publish,
23 * distribute, sublicense, and/or sell copies of the Software, and to
24 * permit persons to whom the Software is furnished to do so, subject to
25 * the following conditions:
26 *
27 * The above copyright notice and this permission notice shall be
28 * included in all copies or substantial portions of the Software.
29 *
30 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
31 * EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
32 * MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
33 * NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE
34 * LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
35 * OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
36 * WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
37 *
38 */
39
40namespace Gecode { namespace Set { namespace RelOp {
41
42 template<class View0, class View1, class View2>
43 Actor*
44 SubOfUnion<View0,View1,View2>::copy(Space& home) {
45 return new (home) SubOfUnion(home,*this);
46 }
47
48 template<class View0, class View1, class View2>
49 ExecStatus
50 SubOfUnion<View0,View1,View2>::propagate(Space& home, const ModEventDelta& med) {
51
52 bool allassigned = x0.assigned() && x1.assigned() && x2.assigned();
53
54 ModEvent me0 = View0::me(med);
55 ModEvent me1 = View1::me(med);
56 ModEvent me2 = View2::me(med);
57
58 bool modified=false;
59
60 do {
61 // lub(x2) <= lub(x0) u lub(x1)
62 if (modified || Rel::testSetEventUB(me0,me1))
63 {
64 LubRanges<View0> ub0(x0);
65 LubRanges<View1> ub1(x1);
66 Iter::Ranges::Union<LubRanges<View0>, LubRanges<View1> > u(ub0,ub1);
67 GECODE_ME_CHECK(x2.intersectI(home,u));
68 }
69
70 // x1 <= glb(x2)-lub(x0)
71 // x0 <= glb(x2)-lub(x1)
72 if (modified || Rel::testSetEventAnyB(me0,me1,me2)) {
73 {
74 modified = false;
75 GlbRanges<View2> lb2(x2);
76 LubRanges<View0> ub0(x0);
77 Iter::Ranges::Diff<GlbRanges<View2>, LubRanges<View0> >
78 diff(lb2,ub0);
79 GECODE_ME_CHECK_MODIFIED(modified, x1.includeI(home,diff));
80 }
81 {
82 GlbRanges<View2> lb2(x2);
83 LubRanges<View1> ub1(x1);
84 Iter::Ranges::Diff<GlbRanges<View2>, LubRanges<View1> >
85 diff(lb2,ub1);
86 GECODE_ME_CHECK_MODIFIED(modified, x0.includeI(home,diff));
87 }
88 } else {
89 modified = false;
90 }
91
92 // cardinality propagation
93 if (modified ||
94 Rel::testSetEventCard(me0,me1,me2) ||
95 Rel::testSetEventLB(me0,me1)
96 ) {
97 GlbRanges<View0> lb0c(x0);
98 GlbRanges<View1> lb1c(x1);
99 Iter::Ranges::Inter<GlbRanges<View0>, GlbRanges<View1> >
100 inter(lb0c,lb1c);
101
102 unsigned int m = Iter::Ranges::size(inter);
103
104 if (m < x0.cardMax()+x1.cardMax()) {
105 GECODE_ME_CHECK_MODIFIED(modified,
106 x2.cardMax( home,
107 x0.cardMax()+x1.cardMax() - m ) );
108 }
109 if (m + x2.cardMin() > x1.cardMax()) {
110 GECODE_ME_CHECK_MODIFIED(modified,
111 x0.cardMin( home,
112 m+x2.cardMin()-x1.cardMax() ) );
113 }
114 if (m + x2.cardMin() > x0.cardMax()) {
115 GECODE_ME_CHECK_MODIFIED(modified,
116 x1.cardMin( home,
117 m+x2.cardMin()-x0.cardMax() ) );
118 }
119 }
120
121 } while (modified);
122
123 if (shared(x0,x1,x2)) {
124 if (allassigned) {
125 return home.ES_SUBSUMED(*this);
126 } else {
127 return ES_NOFIX;
128 }
129 } else {
130 if (x0.assigned() + x1.assigned() + x2.assigned() >= 2) {
131 return home.ES_SUBSUMED(*this);
132 } else {
133 return ES_FIX;
134 }
135 }
136
137 }
138
139 template<class View0, class View1, class View2>
140 forceinline
141 SubOfUnion<View0,View1,View2>::SubOfUnion(Home home, View0 y0,
142 View1 y1, View2 y2)
143 : MixTernaryPropagator<View0,PC_SET_ANY,View1,PC_SET_ANY,
144 View2,PC_SET_ANY>(home,y0,y1,y2) {}
145
146 template<class View0, class View1, class View2>
147 forceinline
148 SubOfUnion<View0,View1,View2>::SubOfUnion
149 (Space& home, SubOfUnion<View0,View1,View2>& p)
150 : MixTernaryPropagator<View0,PC_SET_ANY,View1,PC_SET_ANY,
151 View2,PC_SET_ANY>(home,p) {}
152
153 template<class View0, class View1, class View2>
154 ExecStatus SubOfUnion<View0,View1,View2>::post
155 (Home home, View0 x0, View1 x1, View2 x2) {
156 (void) new (home) SubOfUnion<View0,View1,View2>(home,x0, x1, x2);
157 return ES_OK;
158 }
159
160}}}
161
162// STATISTICS: set-prop