added a symetric round repairer around 0.5, plus a rounding at a given precision

This commit is contained in:
nojhan 2011-11-08 18:40:28 +01:00
commit 0e9d9f8c81
2 changed files with 63 additions and 1 deletions

View file

@ -32,6 +32,7 @@ Authors:
#include "edoRepairerApply.h"
/** A repairer that calls "floor" on each items of a solution
*
* Just a proxy to "edoRepairerApplyUnary<EOT, EOT::AtomType(EOT::AtomType)> rep( std::floor);"
@ -45,6 +46,7 @@ public:
edoRepairerFloor() : edoRepairerApplyUnary<EOT>( std::floor ) {}
};
/** A repairer that calls "ceil" on each items of a solution
*
* @see edoRepairerFloor
@ -59,4 +61,45 @@ public:
};
// FIXME find a way to put this function as a member of edoRepairerRoundDecimals
template< typename ArgType >
ArgType edoRound( ArgType val, ArgType prec = 1.0 )
{
return (val > 0.0) ?
floor(val * prec + 0.5) / prec :
ceil(val * prec - 0.5) / prec ;
}
/** A repairer that round values at a given a precision.
*
* e.g. if prec=10, 8.06 will be rounded to 8.1
*
* @see edoRepairerFloor
* @see edoRepairerCeil
*
* @ingroup Repairers
*/
template < typename EOT >
class edoRepairerRoundDecimals : public edoRepairerApplyBinary<EOT>
{
public:
typedef typename EOT::AtomType ArgType;
edoRepairerRoundDecimals( ArgType prec ) : edoRepairerApplyBinary<EOT>( edoRound<ArgType>, prec ) {}
};
/** A repairer that do a rounding around val+0.5
*
* @see edoRepairerRoundDecimals
*
* @ingroup Repairers
*/
template < typename EOT >
class edoRepairerRound : public edoRepairerRoundDecimals<EOT>
{
public:
edoRepairerRound() : edoRepairerRoundDecimals<EOT>( 1.0 ) {}
};
#endif // !_edoRepairerRound_h

View file

@ -38,10 +38,18 @@ int main(void)
sol.push_back(1.1);
sol.push_back(3.9);
sol.push_back(3.9);
// we expect {1,2,3,4}
sol.push_back(5.4);
sol.push_back(5.6);
sol.push_back(7.011);
sol.push_back(8.09);
sol.push_back(8.21);
std::cout << "expect: INVALID 9 1 2 3 4 5 6 7 8.1 8.2" << std::endl;
edoRepairer<EOT>* rep1 = new edoRepairerFloor<EOT>();
edoRepairer<EOT>* rep2 = new edoRepairerCeil<EOT>();
edoRepairer<EOT>* rep3 = new edoRepairerRound<EOT>();
edoRepairer<EOT>* rep4 = new edoRepairerRoundDecimals<EOT>( 10 );
std::vector<unsigned int> indexes1;
indexes1.push_back(0);
@ -51,8 +59,19 @@ int main(void)
indexes2.push_back(1);
indexes2.push_back(3);
std::vector<unsigned int> indexes3;
indexes3.push_back(4);
indexes3.push_back(5);
std::vector<unsigned int> indexes4;
indexes4.push_back(6);
indexes4.push_back(7);
indexes4.push_back(8);
edoRepairerDispatcher<EOT> repare( indexes1, rep1 );
repare.add( indexes2, rep2 );
repare.add( indexes3, rep3 );
repare.add( indexes4, rep4 );
repare(sol);