linear.cpp
Go to the documentation of this file.
/* -*- mode: C++; c-basic-offset: 2; indent-tabs-mode: nil -*- */
/*
* Main authors:
* Christian Schulte <schulte@gecode.org>
* Vincent Barichard <Vincent.Barichard@univ-angers.fr>
*
* Copyright:
* Christian Schulte, 2005
* Vincent Barichard, 2012
*
* This file is part of Gecode, the generic constraint
* development environment:
* http://www.gecode.org
*
* Permission is hereby granted, free of charge, to any person obtaining
* a copy of this software and associated documentation files (the
* "Software"), to deal in the Software without restriction, including
* without limitation the rights to use, copy, modify, merge, publish,
* distribute, sublicense, and/or sell copies of the Software, and to
* permit persons to whom the Software is furnished to do so, subject to
* the following conditions:
*
* The above copyright notice and this permission notice shall be
* included in all copies or substantial portions of the Software.
*
* THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
* EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
* MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
* NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE
* LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
* OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
* WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
*
*/
#include "test/float.hh"
#include <gecode/minimodel.hh>
namespace Test { namespace Float {
namespace Linear {
bool one(const Gecode::FloatValArgs& a) {
for (int i=a.size(); i--; )
if (a[i] != 1)
return false;
return true;
}
class FloatFloat : public Test {
protected:
Gecode::FloatValArgs a;
Gecode::FloatRelType frt;
Gecode::FloatNum c;
public:
FloatFloat(const std::string& s, const Gecode::FloatVal& d,
const Gecode::FloatValArgs& a0, Gecode::FloatRelType frt0,
Gecode::FloatNum c0, Gecode::FloatNum st)
: Test("Linear::Float::"+
str(frt0)+"::"+s+"::"+str(c0)+"::"
+str(a0.size()),
a0.size(),d,st,CPLT_ASSIGNMENT,true),
a(a0), frt(frt0), c(c0) {
using namespace Gecode;
testfix = false;
if ((frt == FRT_NQ) || (frt == FRT_LE) || (frt == FRT_GR) || reified)
testsubsumed = false;
}
virtual MaybeType solution(const Assignment& x) const {
Gecode::FloatVal e = 0.0;
for (int i=x.size(); i--; )
e += a[i]*x[i];
switch (cmp(e, frt, Gecode::FloatVal(c))) {
case MT_FALSE: {
Gecode::FloatVal eError = e;
for (int i=x.size(); i--; )
eError -= a[i]*x[i];
if (cmp(e+eError, frt, Gecode::FloatVal(c)) == MT_FALSE)
return MT_FALSE;
else
return MT_MAYBE;
}
case MT_TRUE:
return MT_TRUE;
case MT_MAYBE:
return MT_MAYBE;
}
GECODE_NEVER;
return MT_FALSE;
}
virtual void post(Gecode::Space& home, Gecode::FloatVarArray& x) {
if (one(a))
Gecode::linear(home, x, frt, c);
else
Gecode::linear(home, a, x, frt, c);
}
virtual void post(Gecode::Space& home, Gecode::FloatVarArray& x,
Gecode::Reify r) {
if (one(a))
Gecode::linear(home, x, frt, c, r);
else
Gecode::linear(home, a, x, frt, c, r);
}
};
class FloatVar : public Test {
protected:
Gecode::FloatValArgs a;
Gecode::FloatRelType frt;
public:
FloatVar(const std::string& s, const Gecode::FloatVal& d,
const Gecode::FloatValArgs& a0, Gecode::FloatRelType frt0, Gecode::FloatNum st)
: Test("Linear::Var::"+
str(frt0)+"::"+s+"::"+str(a0.size()),
a0.size()+1,d,st,CPLT_ASSIGNMENT,true),
a(a0), frt(frt0) {
using namespace Gecode;
testfix = false;
if ((frt == FRT_NQ) || (frt == FRT_LE) || (frt == FRT_GR) || reified)
testsubsumed = false;
}
virtual MaybeType solution(const Assignment& x) const {
Gecode::FloatVal e = 0.0;
for (int i=a.size(); i--; )
e += a[i]*x[i];
switch (cmp(e, frt, x[a.size()])) {
case MT_FALSE: {
Gecode::FloatVal eError = e;
for (int i=a.size(); i--; )
eError -= a[i]*x[i];
if (cmp(e+eError, frt, x[a.size()]) == MT_FALSE)
return MT_FALSE;
else
return MT_MAYBE;
}
case MT_TRUE:
return MT_TRUE;
case MT_MAYBE:
return MT_MAYBE;
}
GECODE_NEVER;
return MT_FALSE;
}
virtual void post(Gecode::Space& home, Gecode::FloatVarArray& x) {
int n = a.size();
Gecode::FloatVarArgs y(n);
for (int i=n; i--; )
y[i] = x[i];
if (one(a))
Gecode::linear(home, y, frt, x[n]);
else
Gecode::linear(home, a, y, frt, x[n]);
}
virtual void post(Gecode::Space& home, Gecode::FloatVarArray& x,
Gecode::Reify r) {
int n = a.size();
Gecode::FloatVarArgs y(n);
for (int i=n; i--; )
y[i] = x[i];
if (one(a))
Gecode::linear(home, y, frt, x[n], r);
else
Gecode::linear(home, a, y, frt, x[n], r);
}
};
class Create {
public:
Create(void) {
using namespace Gecode;
{
FloatNum step = 0.7;
FloatVal f1(-2,2);
FloatVal f2(-3,-1);
FloatVal f3(3,8);
FloatValArgs a1(1);
a1[0] = 0.0;
for (FloatRelTypes frts; frts(); ++frts) {
(void) new FloatFloat("11",f1,a1,frts.frt(),0.0,step);
(void) new FloatVar("11",f1,a1,frts.frt(),step);
(void) new FloatFloat("21",f2,a1,frts.frt(),0.0,step);
(void) new FloatVar("21",f2,a1,frts.frt(),step);
(void) new FloatFloat("31",f3,a1,frts.frt(),1.0,step);
}
const FloatVal av2[4] = {1.0,1.0,1.0,1.0};
const FloatVal av3[4] = {1.0,-1.0,-1.0,1.0};
const FloatVal av4[4] = {2.0,3.0,5.0,7.0};
const FloatVal av5[4] = {-2.0,3.0,-5.0,7.0};
for (int i=1; i<=4; i++) {
FloatValArgs a2(i, av2);
FloatValArgs a3(i, av3);
FloatValArgs a4(i, av4);
FloatValArgs a5(i, av5);
for (FloatRelTypes frts; frts(); ++frts) {
(void) new FloatFloat("12",f1,a2,frts.frt(),0.0,step);
(void) new FloatFloat("13",f1,a3,frts.frt(),0.0,step);
(void) new FloatFloat("14",f1,a4,frts.frt(),0.0,step);
(void) new FloatFloat("15",f1,a5,frts.frt(),0.0,step);
(void) new FloatFloat("22",f2,a2,frts.frt(),0.0,step);
(void) new FloatFloat("23",f2,a3,frts.frt(),0.0,step);
(void) new FloatFloat("24",f2,a4,frts.frt(),0.0,step);
(void) new FloatFloat("25",f2,a5,frts.frt(),0.0,step);
(void) new FloatFloat("32",f3,a2,frts.frt(),1.0,step);
if (i < 4) {
(void) new FloatVar("12",f1,a2,frts.frt(),step);
(void) new FloatVar("13",f1,a3,frts.frt(),step);
(void) new FloatVar("14",f1,a4,frts.frt(),step);
(void) new FloatVar("15",f1,a5,frts.frt(),step);
(void) new FloatVar("22",f2,a2,frts.frt(),step);
(void) new FloatVar("23",f2,a3,frts.frt(),step);
(void) new FloatVar("24",f2,a4,frts.frt(),step);
(void) new FloatVar("25",f2,a5,frts.frt(),step);
}
}
}
}
}
};
Create c;
}
}}
// STATISTICS: test-float
Gecode::FloatValArgs
Passing float arguments.
Definition float.hh:950
Gecode::FloatVal
Float value type.
Definition float.hh:334
Gecode::FloatVal::size
FloatNum size(void) const
Return size of float value (distance between maximum and minimum).
Definition val.hpp:78
Gecode::FloatVarArgs
Passing float variables.
Definition float.hh:979
Gecode::FloatVarArray
Float variable array.
Definition float.hh:1030
Gecode::Reify
Reification specification.
Definition int.hh:876
Gecode::Space
Computation spaces.
Definition core.hpp:1742
Test::Float::Assignment
Base class for assignments
Definition float.hh:80
Test::Float::Assignment::size
int size(void) const
Return number of variables.
Definition float.hpp:48
Test::Float::FloatRelTypes
Iterator for float relation types.
Definition float.hh:339
Test::Float::Linear::Create
Help class to create and register tests.
Definition linear.cpp:191
Test::Float::Linear::Create::Create
Create(void)
Perform creation and registration.
Definition linear.cpp:194
Test::Float::Linear::FloatFloat
Test linear relation over float variables
Definition linear.cpp:59
Test::Float::Linear::FloatFloat::solution
virtual MaybeType solution(const Assignment &x) const
Test whether x is solution
Definition linear.cpp:83
Test::Float::Linear::FloatFloat::FloatFloat
FloatFloat(const std::string &s, const Gecode::FloatVal &d, const Gecode::FloatValArgs &a0, Gecode::FloatRelType frt0, Gecode::FloatNum c0, Gecode::FloatNum st)
Create and register test.
Definition linear.cpp:69
Test::Float::Linear::FloatFloat::post
virtual void post(Gecode::Space &home, Gecode::FloatVarArray &x)
Post constraint on x.
Definition linear.cpp:106
Test::Float::Linear::FloatFloat::post
virtual void post(Gecode::Space &home, Gecode::FloatVarArray &x, Gecode::Reify r)
Post reified constraint on x for r.
Definition linear.cpp:113
Test::Float::Linear::FloatFloat::frt
Gecode::FloatRelType frt
Float relation type to propagate.
Definition linear.cpp:64
Test::Float::Linear::FloatFloat::a
Gecode::FloatValArgs a
Coefficients.
Definition linear.cpp:62
Test::Float::Linear::FloatFloat::c
Gecode::FloatNum c
Result.
Definition linear.cpp:66
Test::Float::Linear::FloatVar
Test linear relation over float variables
Definition linear.cpp:123
Test::Float::Linear::FloatVar::frt
Gecode::FloatRelType frt
Float relation type to propagate.
Definition linear.cpp:128
Test::Float::Linear::FloatVar::post
virtual void post(Gecode::Space &home, Gecode::FloatVarArray &x)
Post constraint on x.
Definition linear.cpp:166
Test::Float::Linear::FloatVar::solution
virtual MaybeType solution(const Assignment &x) const
Test whether x is solution
Definition linear.cpp:143
Test::Float::Linear::FloatVar::post
virtual void post(Gecode::Space &home, Gecode::FloatVarArray &x, Gecode::Reify r)
Post reified constraint on x for r.
Definition linear.cpp:177
Test::Float::Linear::FloatVar::a
Gecode::FloatValArgs a
Coefficients.
Definition linear.cpp:126
Test::Float::Linear::FloatVar::FloatVar
FloatVar(const std::string &s, const Gecode::FloatVal &d, const Gecode::FloatValArgs &a0, Gecode::FloatRelType frt0, Gecode::FloatNum st)
Create and register test.
Definition linear.cpp:131
Test::Float::Test::cmp
static MaybeType cmp(Gecode::FloatVal x, Gecode::FloatRelType r, Gecode::FloatVal y)
Compare x and y with respect to r.
Definition float.hpp:235
Test::Float::Test::testsubsumed
bool testsubsumed
Whether to test for subsumption.
Definition float.hh:263
Test::Float::Test::testfix
bool testfix
Whether to perform fixpoint test.
Definition float.hh:261
Test::Float::Test::str
static std::string str(Gecode::FloatRelType frt)
Map float relation to string.
Definition float.hpp:194
Test::Float::Test::reified
bool reified
Does the constraint also exist as reified constraint.
Definition float.hh:255
Gecode::linear
void linear(Home home, const FloatVarArgs &x, FloatRelType frt, FloatVal c)
Post propagator for .
Definition linear.cpp:41
Gecode::FloatNum
double FloatNum
Floating point number base type.
Definition float.hh:106
Gecode::FloatRelType
FloatRelType
Relation types for floats.
Definition float.hh:1068
Gecode::FRT_GR
@ FRT_GR
Greater ( ).
Definition float.hh:1074
Gecode::FRT_NQ
@ FRT_NQ
Disequality ( ).
Definition float.hh:1070
Gecode::FRT_LE
@ FRT_LE
Less ( ).
Definition float.hh:1072
minimodel.hh
Gecode
Gecode toplevel namespace
Test::Float::Linear
Tests for linear constraints
Definition linear.cpp:43
Test::Float::Linear::c
Create c
Definition linear.cpp:249
Test::Float::Linear::one
bool one(const Gecode::FloatValArgs &a)
Check whether has only one coefficients.
Definition linear.cpp:46
Test::Float
Testing domain floats.
Definition float.cpp:43
Test::Float::CPLT_ASSIGNMENT
@ CPLT_ASSIGNMENT
Definition float.hh:62
Test::Float::MaybeType
MaybeType
Type for comparisons and solutions.
Definition float.hh:51
Test::Float::MT_MAYBE
@ MT_MAYBE
Definition float.hh:54
Test::Float::MT_FALSE
@ MT_FALSE
Definition float.hh:52
Test::Float::MT_TRUE
@ MT_TRUE
Definition float.hh:53
Test
General test support.
Definition afc.cpp:39
Test::r
Region r
Definition region.cpp:65
GECODE_NEVER
#define GECODE_NEVER
Assert that this command is never executed.
Definition macros.hpp:56
float.hh