testBasicPointFunctors.cpp

#include <cstdio>
#include <cmath>
#include <iostream>
#include <fstream>
#include <vector>
#include "DGtal/base/Common.h"
#include "DGtal/base/CUnaryFunctor.h"
#include "DGtal/kernel/PointVector.h"
#include "DGtal/kernel/SpaceND.h"
#include "DGtal/kernel/BasicPointFunctors.h"

using namespace DGtal;
using namespace std;

template <typename TFunctor, typename TArg, typename TRes >
void checkingConcepts()
{
  BOOST_CONCEPT_ASSERT(( CUnaryFunctor<TFunctor, TArg, TRes > ));
}

bool testProjector()
{
  unsigned int nbok = 0;
  unsigned int nb = 0;

  trace.beginBlock ( "Checking projection on a subspace" );
  {
  //a 3d point
  PointVector<3,int> p(0,1,2); 

  //projectors
  typedef Projector<SpaceND<2,int> > Projector2D; 
  std::vector<Dimension> v1, v2; 
  v1.push_back(0);   v1.push_back(2); 
  v2.push_back(2);   v2.push_back(1);
  Projector2D proj1, proj2, proj3; 
  proj1.init(v1.begin(), v1.end()); 
  proj2.init(v2.begin(), v2.end()); 
 
  //comparison
  PointVector<2,int> res1(0,2); 
  trace.info() << "p " << p << " => " << proj1(p) << " == " << res1 << std::endl;
  nbok += ( proj1(p) == res1 ) ? 1 : 0; 
  nb++;  

  PointVector<2,int> res2(2,1);
  trace.info() << "p " << p << " => " << proj2(p) << " == " << res2 << std::endl;
  nbok += ( proj2(p) == res2 ) ? 1 : 0; 
  nb++;  

  PointVector<2,int> res3(0,1);
  trace.info() << "p " << p << " => " << proj3(p) << " == " << res3 << std::endl;
  nbok += ( proj3(p) == res3 ) ? 1 : 0; 
  nb++;  
  }
  trace.endBlock();

  trace.beginBlock ( "Checking projection on a greater space" );
  {
  //a 2d point
  PointVector<2,int> p(5,2); 

  //projectors
  typedef Projector<SpaceND<3,int> > Projector3D; 
  std::vector<Dimension> v1, v2, v4; 
  v1.push_back(0);   v1.push_back(2); v1.push_back(1);  
  v2.push_back(1);   v2.push_back(0);
  v4.push_back(1);
  Projector3D proj1, proj2, proj3; 
  proj1.init(v1.begin(), v1.end()); 
  proj2.init(v2.begin(), v2.end());
  Projector3D proj4(-1);  
  proj4.init(v4.begin(), v4.end()); 

  //comparison
  PointVector<3,int> res1(5,0,2); 
  trace.info() << "p " << p << " => " << proj1(p) << " == " << res1 << std::endl;
  nbok += ( proj1(p) == res1 ) ? 1 : 0; 
  nb++;  

  PointVector<3,int> res2(2,5,0);
  trace.info() << "p " << p << " => " << proj2(p) << " == " << res2 << std::endl;
  nbok += ( proj2(p) == res2 ) ? 1 : 0; 
  nb++;  

  PointVector<3,int> res3(5,2,0);
  trace.info() << "p " << p << " => " << proj3(p) << " == " << res3 << std::endl;
  nbok += ( proj3(p) == res3 ) ? 1 : 0; 
  nb++;  

  PointVector<3,int> res4(2,-1,-1);
  trace.info() << "p " << p << " => " << proj4(p) << " == " << res4
               << "(-1 as default value)" << std::endl;
  nbok += ( proj4(p) == res4 ) ? 1 : 0; 
  nb++;  
  }
  trace.endBlock();

  return nbok == nb;
}

// Standard services - public :

int main( int argc, char** argv )
{
  trace.beginBlock ( "Testing basic point functors" );
  trace.info() << "Args:";
  for ( int i = 0; i < argc; ++i )
    trace.info() << " " << argv[ i ];
  trace.info() << endl;


  checkingConcepts<Projector<SpaceND<2,int> >, PointVector<6,int>, PointVector<2,int> >(); 
  //for instance, this does not compile because 
  //the point of dim 6 is projected on a point of dim 2 (and not 3)
  //checkingConcepts<Projector<SpaceND<2,int> >, PointVector<6,int>, PointVector<3,int> >(); 

  bool res = testProjector();

  trace.emphase() << ( res ? "Passed." : "Error." ) << endl;
  trace.endBlock();
  return res ? 0 : 1;
}