332 lines
7.7 KiB
C++
332 lines
7.7 KiB
C++
//-----------------------------------------------------------------------------
|
|
// t-moeoQuickUnboundedArchive.cpp
|
|
//-----------------------------------------------------------------------------
|
|
|
|
#include <eo>
|
|
#include <moeo>
|
|
|
|
//-----------------------------------------------------------------------------
|
|
|
|
class ObjectiveVectorTraits : public moeoObjectiveVectorTraits
|
|
{
|
|
public:
|
|
static bool minimizing (int i)
|
|
{
|
|
return true;
|
|
}
|
|
static bool maximizing (int i)
|
|
{
|
|
return false;
|
|
}
|
|
static unsigned int nObjectives ()
|
|
{
|
|
return 2;
|
|
}
|
|
};
|
|
|
|
class ObjectiveVectorTraits3d : public moeoObjectiveVectorTraits
|
|
{
|
|
public:
|
|
static bool minimizing (int i)
|
|
{
|
|
return true ;
|
|
}
|
|
static bool maximizing (int i)
|
|
{
|
|
return false ;
|
|
}
|
|
static unsigned int nObjectives ()
|
|
{
|
|
return 3;
|
|
}
|
|
};
|
|
|
|
|
|
typedef moeoRealObjectiveVector < ObjectiveVectorTraits > ObjectiveVector;
|
|
typedef moeoRealObjectiveVector < ObjectiveVectorTraits3d > ObjectiveVector3d;
|
|
|
|
typedef MOEO < ObjectiveVector, double, double > Solution;
|
|
typedef MOEO < ObjectiveVector3d, double, double > Solution3d;
|
|
|
|
//-----------------------------------------------------------------------------
|
|
|
|
int main()
|
|
{
|
|
std::cout << "[moeoArchive]\t=>\t";
|
|
// First test is just to verify behavior of moeoQuickUnboundedArchive
|
|
// objective vectors
|
|
ObjectiveVector obj0, obj1, obj2, obj3, obj4, obj5,obj10;
|
|
obj0[0] = 12226;
|
|
obj0[1] = 427894;
|
|
obj1[0] = 12170;
|
|
obj1[1] = 431736;
|
|
obj2[0] = 11965;
|
|
obj2[1] = 435193;
|
|
obj3[0] = 11893;
|
|
obj3[1] = 441839;
|
|
obj4[0] = 11870;
|
|
obj4[1] = 450770;
|
|
obj5[0] = 11769;
|
|
obj5[1] = 460005;
|
|
obj10[0] = 11769;
|
|
obj10[1] = 46005;
|
|
// population
|
|
eoPop < Solution > pop;
|
|
pop.resize(6);
|
|
pop[0].objectiveVector(obj0);
|
|
pop[1].objectiveVector(obj1);
|
|
pop[2].objectiveVector(obj2);
|
|
pop[3].objectiveVector(obj3);
|
|
pop[4].objectiveVector(obj4);
|
|
pop[5].objectiveVector(obj5);
|
|
|
|
// archive
|
|
moeoQuadTree<Solution> index;
|
|
moeoIndexedArchive <Solution> arch(index);
|
|
// moeoQuickUnboundedArchive< Solution > arch;
|
|
arch(pop);
|
|
// arch.printIndex();
|
|
|
|
// size
|
|
if (arch.size() > 6)
|
|
{
|
|
std::cout << "ERROR1 (too much solutions)" << std::endl;
|
|
return EXIT_FAILURE;
|
|
}
|
|
// obj0 must be in
|
|
if (! arch.contains(obj0))
|
|
{
|
|
std::cout << "ERROR2 (obj0 not in)" << std::endl;
|
|
return EXIT_FAILURE;
|
|
}
|
|
// obj1 must be in
|
|
if (! arch.contains(obj1))
|
|
{
|
|
std::cout << "ERROR3 (obj1 not in)" << std::endl;
|
|
return EXIT_FAILURE;
|
|
}
|
|
// obj2 must be in
|
|
if (! arch.contains(obj2))
|
|
{
|
|
std::cout << "ERROR4 (obj2 not in)" << std::endl;
|
|
return EXIT_FAILURE;
|
|
}
|
|
// obj3 must be in
|
|
if (! arch.contains(obj3))
|
|
{
|
|
std::cout << "ERROR5 (obj3 not in)" << std::endl;
|
|
return EXIT_FAILURE;
|
|
}
|
|
// obj4 must be in
|
|
if (!arch.contains(obj4))
|
|
{
|
|
std::cout << "ERROR6 (obj4 not in)! " << obj4<< std::endl;
|
|
// arch.printIndex();
|
|
return EXIT_FAILURE;
|
|
}
|
|
// obj5 must be in
|
|
if (! arch.contains(obj5))
|
|
{
|
|
std::cout << "ERROR7 (obj5 not in)" << std::endl;
|
|
return EXIT_FAILURE;
|
|
}
|
|
|
|
// Second test is to verify behavior with negative and comparator
|
|
ObjectiveVector obj6, obj7;
|
|
obj6[0] = -12170;
|
|
obj6[1] = 427894;
|
|
obj7[0] = -12226;
|
|
obj7[1] = 431736;
|
|
eoPop < Solution > pop2;
|
|
pop2.resize(2);
|
|
pop2[0].objectiveVector(obj6);
|
|
pop2[1].objectiveVector(obj7);
|
|
moeoQuadTree< Solution > index2;
|
|
moeoIndexedArchive< Solution > arch2(index2);
|
|
arch2(pop2);
|
|
|
|
// size
|
|
if (arch2.size() != 2)
|
|
{
|
|
std::cout << "ERROR8 (too much solutions)" << std::endl;
|
|
return EXIT_FAILURE;
|
|
}
|
|
|
|
//Third test is with two equals values
|
|
ObjectiveVector obj8, obj9;
|
|
obj8[0] = 10;
|
|
obj8[1] = 10;
|
|
obj9[0] = 10;
|
|
obj9[1] = 10;
|
|
eoPop < Solution > pop3;
|
|
pop3.resize(2);
|
|
pop3[0].objectiveVector(obj8);
|
|
pop3[1].objectiveVector(obj9);
|
|
|
|
moeoQuadTree< Solution > index3;
|
|
moeoIndexedArchive< Solution > arch3(index3);
|
|
arch3(pop3);
|
|
|
|
if (arch3.size() != 1)
|
|
{
|
|
std::cout << "ERROR9 (too much solutions)" << std::endl;
|
|
return EXIT_FAILURE;
|
|
}
|
|
eoPop < Solution > pop4;
|
|
pop4.resize(6);
|
|
pop4[0].objectiveVector(obj0);
|
|
pop4[1].objectiveVector(obj1);
|
|
pop4[2].objectiveVector(obj2);
|
|
pop4[3].objectiveVector(obj3);
|
|
pop4[4].objectiveVector(obj4);
|
|
pop4[5].objectiveVector(obj10);
|
|
moeoQuadTree< Solution > index4;
|
|
moeoIndexedArchive< Solution > arch4(index4);
|
|
arch4(pop4);
|
|
if (arch4.size() != 1)
|
|
{
|
|
std::cout << "ERROR10 (too much solutions)" << std::endl;
|
|
return EXIT_FAILURE;
|
|
}
|
|
|
|
std::vector<bool>bobj;
|
|
for (unsigned int i=0;i<3;i++){
|
|
bobj.push_back(false);
|
|
}
|
|
|
|
ObjectiveVectorTraits3d::setup(3,bobj);
|
|
moeoQuadTree < Solution3d> index5;
|
|
moeoIndexedArchive < Solution3d> arch5(index5);
|
|
eoPop<Solution3d> pop5;
|
|
pop5.resize(6);
|
|
ObjectiveVector3d obj3d0, obj3d1, obj3d2, obj3d3, obj3d4, obj3d5, obj3d6,obj3d10;
|
|
obj3d0[0] = 12226;
|
|
obj3d0[1] = 427894;
|
|
obj3d0[2] = 10;
|
|
|
|
obj3d1[0] = 12170;
|
|
obj3d1[1] = 431736;
|
|
obj3d1[2] = 10;
|
|
|
|
obj3d2[0] = 11965;
|
|
obj3d2[1] = 435193;
|
|
obj3d2[2] = 10;
|
|
|
|
obj3d3[0] = 11893;
|
|
obj3d3[1] = 441839;
|
|
obj3d3[2] = 10;
|
|
|
|
obj3d4[0] = 11870;
|
|
obj3d4[1] = 450770;
|
|
obj3d4[2] = 10;
|
|
|
|
obj3d5[0] = 11769;
|
|
obj3d5[1] = 46005;
|
|
obj3d5[2] = 2;
|
|
|
|
obj3d10[0] = 11769;
|
|
obj3d10[1] = 46005;
|
|
obj3d10[2] = 10;
|
|
|
|
obj3d6[0] = 11769;
|
|
obj3d6[1] = 460005;
|
|
obj3d6[2] = 1;
|
|
|
|
|
|
pop5[0].objectiveVector(obj3d0);
|
|
pop5[1].objectiveVector(obj3d1);
|
|
pop5[2].objectiveVector(obj3d2);
|
|
pop5[3].objectiveVector(obj3d3);
|
|
pop5[4].objectiveVector(obj3d4);
|
|
pop5[5].objectiveVector(obj3d10);
|
|
|
|
|
|
|
|
arch5(pop5);
|
|
pop5.resize(7);
|
|
pop5[6].objectiveVector(obj3d5);
|
|
arch5(pop5[6]);
|
|
// index5.printTree();
|
|
/* for (unsigned int i=0;i<arch5.size();i++){
|
|
std::cout<<i<<" "<<arch5[i].objectiveVector()<<std::endl;
|
|
}*/
|
|
assert(arch5.size()==1);
|
|
pop5.resize(8);
|
|
pop5[7].objectiveVector(obj3d6);
|
|
arch5(pop5[7]);
|
|
assert(arch5.size()==2);
|
|
|
|
moeoQuadTree < Solution3d> index6;
|
|
moeoIndexedArchive < Solution3d> arch6(index6);
|
|
eoPop<Solution3d> pop6;
|
|
|
|
ObjectiveVector3d jojo, jojo1, jojo2, jojo3, jojo4, jojo5, jojo6, jojo7, jojo8, jojo9 ;
|
|
jojo[0]=10;
|
|
jojo[1]=10;
|
|
jojo[2]=10;
|
|
|
|
jojo1[0]=5;
|
|
jojo1[1]=5;
|
|
jojo1[2]=23;
|
|
|
|
jojo2[0]=3;
|
|
jojo2[1]=25;
|
|
jojo2[2]=16;
|
|
|
|
jojo3[0]=14;
|
|
jojo3[1]=18;
|
|
jojo3[2]=6;
|
|
|
|
jojo4[0]=100;
|
|
jojo4[1]=100;
|
|
jojo4[2]=100;
|
|
|
|
jojo5[0]=4;
|
|
jojo5[1]=8;
|
|
jojo5[2]=12;
|
|
|
|
|
|
jojo6[0]=11;
|
|
jojo6[1]=14;
|
|
jojo6[2]=9;
|
|
|
|
|
|
jojo7[0]=40;
|
|
jojo7[1]=12;
|
|
jojo7[2]=3;
|
|
|
|
jojo8[0]=12;
|
|
jojo8[1]=8;
|
|
jojo8[2]=4;
|
|
|
|
jojo9[0]=4;
|
|
jojo9[1]=4;
|
|
jojo9[2]=22;
|
|
|
|
pop6.resize(8);
|
|
pop6[0].objectiveVector(jojo);
|
|
pop6[1].objectiveVector(jojo1);
|
|
pop6[2].objectiveVector(jojo2);
|
|
pop6[3].objectiveVector(jojo3);
|
|
pop6[4].objectiveVector(jojo4);
|
|
pop6[5].objectiveVector(jojo5);
|
|
pop6[6].objectiveVector(jojo6);
|
|
pop6[7].objectiveVector(jojo7);
|
|
Solution3d add;
|
|
Solution3d add2;
|
|
add.objectiveVector(jojo8);
|
|
add2.objectiveVector(jojo9);
|
|
/* pop[1].objectiveVector(jojo1);
|
|
pop[1].objectiveVector(jojo1);
|
|
pop[1].objectiveVector(jojo1);*/
|
|
arch6(pop6);
|
|
|
|
assert(arch6.size()==7);
|
|
bool res=arch6(add);
|
|
assert(res && arch6.size()==7);
|
|
res=arch6(add2);
|
|
assert(res && arch6.size()==7);
|
|
std::cout << "OK" << std::endl;
|
|
return EXIT_SUCCESS;
|
|
}
|
|
|