testAdjacency.cpp

#include <iostream>
#include <vector>
#include "DGtal/helpers/StdDefs.h"

#include "DGtal/base/Common.h"
#include "DGtal/kernel/SpaceND.h"
#include "DGtal/topology/MetricAdjacency.h"
#include "DGtal/topology/CUndirectedSimpleLocalGraph.h"



using namespace std;
using namespace DGtal;


// Functions for testing class Adjacency.

bool testMetricAdjacency()
{
  unsigned int nbok = 0;
  unsigned int nb = 0;
  
  typedef SpaceND<3> Space3D;
  typedef Space3D::Point Point;
  typedef Z3i::Adj6 Adj6;
  typedef Z3i::Adj18 Adj18;
  typedef Z3i::Adj26 Adj26;
  Point p( 3, -5, 10 );


  trace.beginBlock ( "Testing neighbors of" );
  Adj6::selfDisplay( trace.info() );
  trace.info() << " p = " << p << std::endl;
  vector<Point> neigh6;
  back_insert_iterator< vector<Point> > bii6( neigh6 );
  Adj6::writeNeighbors( bii6, p );
  nbok += neigh6.size() == 6 ? 1 : 0; 
  nb++;
  trace.info() << "(" << nbok << "/" << nb << ") "
         << "Card(6-neigh): " << neigh6.size() 
         << "== 6 ?" << std::endl;
  trace.beginBlock ( "Enumerating neighbors." );
  unsigned int nb_correct = 0;
  for ( unsigned int i = 0; i < neigh6.size(); ++i )
    {
      if ( Adj6::isProperlyAdjacentTo( p, neigh6[ i ] ) )
  {
    trace.info() << neigh6[ i ] << "* " << std::endl;
    ++nb_correct;
  }
      else
  trace.info() << neigh6[ i ] << "- " << std::endl;
    }
  trace.endBlock();
  nbok += nb_correct == 6 ? 1 : 0; 
  nb++;
  trace.info() << "(" << nbok << "/" << nb << ") "
         << "Within, #proper adjacent : " << nb_correct
         << "== 6 ?" << std::endl;
  trace.endBlock();



  trace.beginBlock ( "Testing neighborhood of" );
  Adj18::selfDisplay( trace.info() );
  trace.info() << " p = " << p << std::endl;
  vector<Point> neigh18;
  back_insert_iterator< vector<Point> > bii18( neigh18 );
  Adj18::writeNeighbors( bii18, p );
  nbok += neigh18.size() == 18 ? 1 : 0; 
  nb++;
  trace.info() << "(" << nbok << "/" << nb << ") "
         << "Card(18-neigh): " << neigh18.size() 
         << "== 18 ?" << std::endl;
  trace.beginBlock ( "Enumerating neighbors." );
  nb_correct = 0;
  for ( unsigned int i = 0; i < neigh18.size(); ++i )
    {
      if ( Adj18::isProperlyAdjacentTo( p, neigh18[ i ] ) )
  {
    trace.info() << neigh18[ i ] << "* " << std::endl;
    ++nb_correct;
  }
      else
  trace.info() << neigh18[ i ] << "- " << std::endl;
    }
  trace.endBlock();
  nbok += nb_correct == 18 ? 1 : 0; 
  nb++;
  trace.info() << "(" << nbok << "/" << nb << ") "
         << "Within, #proper adjacent : " << nb_correct
         << "== 18 ?" << std::endl;
  trace.endBlock();


  trace.beginBlock ( "Testing neighborhood of" );
  Adj26::selfDisplay( trace.info() );
  trace.info() << " p = " << p << std::endl;
  vector<Point> neigh26;
  back_insert_iterator< vector<Point> > bii26( neigh26 );
  Adj26::writeNeighbors( bii26, p );
  nbok += neigh26.size() == 26 ? 1 : 0; 
  nb++;
  trace.info() << "(" << nbok << "/" << nb << ") "
         << "Card(26-neigh): " << neigh26.size() 
         << "== 26 ?" << std::endl;
  trace.beginBlock ( "Enumerating neighbors." );
  nb_correct = 0;
  for ( unsigned int i = 0; i < neigh26.size(); ++i )
    {
      if ( Adj26::isProperlyAdjacentTo( p, neigh26[ i ] ) )
  {
    trace.info() << neigh26[ i ] << "* " << std::endl;
    ++nb_correct;
  }
      else
  trace.info() << neigh26[ i ] << "- " << std::endl;
    }
  trace.endBlock();
  nbok += nb_correct == 26 ? 1 : 0; 
  nb++;
  trace.info() << "(" << nbok << "/" << nb << ") "
         << "Within, #proper adjacent : " << nb_correct
         << "== 26 ?" << std::endl;
  trace.endBlock();

  return nbok == nb;
}


bool testLocalGraphModel()
{
  trace.beginBlock ( "Testing graph model" );
  unsigned int nbok=0,nb=0;

  typedef DGtal::MetricAdjacency<SpaceND<6,int>, 2>  Adj;
  BOOST_CONCEPT_ASSERT(( CUndirectedSimpleLocalGraph<Adj> ));
  
  
  nbok += Adj::bestCapacity() == 72 ? 1 : 0; 
  nb++;
  trace.info() << "(" << nbok << "/" << nb << ") "
         << "Within, bestCapacity : " << Adj::bestCapacity()
         << "== 72 ?" << std::endl;
  trace.endBlock();

  
  
  return nbok == nb;
}


// Standard services - public :

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

  bool res = testMetricAdjacency() && testLocalGraphModel(); // && ... other tests
  trace.emphase() << ( res ? "Passed." : "Error." ) << endl;
  trace.endBlock();
  return res ? 0 : 1;
}
//                                                                           //