hull.cpp
Go to the documentation of this file.
00001 /* -*- mode: C++; c-basic-offset: 2; indent-tabs-mode: nil -*- */ 00002 /* 00003 * Main authors: 00004 * Guido Tack <tack@gecode.org> 00005 * Christian Schulte <schulte@gecode.org> 00006 * 00007 * Contributing authors: 00008 * Gabor Szokoli <szokoli@gecode.org> 00009 * 00010 * Copyright: 00011 * Guido Tack, 2004 00012 * Christian Schulte, 2004 00013 * Gabor Szokoli, 2004 00014 * 00015 * Last modified: 00016 * $Date: 2010-07-28 16:13:53 +0200 (Wed, 28 Jul 2010) $ by $Author: schulte $ 00017 * $Revision: 11292 $ 00018 * 00019 * This file is part of Gecode, the generic constraint 00020 * development environment: 00021 * http://www.gecode.org 00022 * 00023 * Permission is hereby granted, free of charge, to any person obtaining 00024 * a copy of this software and associated documentation files (the 00025 * "Software"), to deal in the Software without restriction, including 00026 * without limitation the rights to use, copy, modify, merge, publish, 00027 * distribute, sublicense, and/or sell copies of the Software, and to 00028 * permit persons to whom the Software is furnished to do so, subject to 00029 * the following conditions: 00030 * 00031 * The above copyright notice and this permission notice shall be 00032 * included in all copies or substantial portions of the Software. 00033 * 00034 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, 00035 * EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF 00036 * MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND 00037 * NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE 00038 * LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION 00039 * OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION 00040 * WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE. 00041 * 00042 */ 00043 00044 #include <gecode/set/convex.hh> 00045 00046 namespace Gecode { namespace Set { namespace Convex { 00047 00048 Actor* 00049 ConvexHull::copy(Space& home, bool share) { 00050 return new (home) ConvexHull(home,share,*this); 00051 } 00052 00053 ExecStatus 00054 ConvexHull::propagate(Space& home, const ModEventDelta&) { 00055 //x1 is the convex hull of x0 00056 00057 GECODE_ME_CHECK( x1.cardMin(home,x0.cardMin()) ); 00058 GECODE_ME_CHECK( x0.cardMax(home,x1.cardMax()) ); 00059 00060 do { 00061 00062 //intersect x1 with (x0.lubMin(),x0.lubMax()) 00063 //This empties x1 if x0.ub is empty. twice. 00064 GECODE_ME_CHECK( x1.exclude(home,Limits::min, 00065 x0.lubMin()-1) ); 00066 GECODE_ME_CHECK( x1.exclude(home,x0.lubMax()+1, 00067 Limits::max) ); 00068 00069 int minElement = std::min(x1.glbMin(),x0.glbMin()); 00070 int maxElement = std::max(x1.glbMax(),x0.glbMax()); 00071 00072 if (minElement<maxElement) { 00073 GECODE_ME_CHECK( x1.include(home, minElement, maxElement)); 00074 } 00075 00076 unsigned int cardMin = x1.cardMin(); 00077 00078 Region r(home); 00079 LubRanges<SetView> ubRangeIt(x1); 00080 Iter::Ranges::Cache< LubRanges<SetView> > ubRangeItC(r,ubRangeIt); 00081 for (;ubRangeItC();++ubRangeItC){ 00082 if (ubRangeItC.width() < cardMin 00083 || ubRangeItC.min() > minElement //No need to test for empty lb. 00084 || ubRangeItC.max() < maxElement 00085 ) { 00086 GECODE_ME_CHECK( x1.exclude(home, 00087 ubRangeItC.min(), ubRangeItC.max()) ); 00088 } 00089 } 00090 00091 LubRanges<SetView> ubRangeIt2(x1); 00092 GECODE_ME_CHECK(x0.intersectI(home,ubRangeIt2) ); 00093 00094 if (x1.lubMin()!=BndSet::MIN_OF_EMPTY) { 00095 if(x1.lubMin()==x1.glbMin()) { 00096 GECODE_ME_CHECK(x0.include(home,x1.lubMin())); 00097 } 00098 if(x1.lubMax()==x1.glbMax()) { 00099 GECODE_ME_CHECK(x0.include(home,x1.lubMax())); 00100 } 00101 } 00102 } while(x0.assigned()&&!x1.assigned()); 00103 00104 //If x0 is assigned, x1 should be too. 00105 assert(x1.assigned() || !x0.assigned()); 00106 00107 if (x1.assigned()) { 00108 return home.ES_SUBSUMED(*this); 00109 } 00110 00111 return ES_NOFIX; 00112 } 00113 00114 }}} 00115 00116 // STATISTICS: set-prop