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 38 namespace Gecode { namespace Float { namespace Arithmetic { 39 40 template<class A, class B> 41 forceinline Abs(Home home,A x0,B x1)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 post(Home home,A x0,B x1)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 Abs(Space & home,Abs<A,B> & p)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* copy(Space & home)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 propagate(Space & home,const ModEventDelta &)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 87