testEuclideanShapesDecorator.cpp

Go to the documentation of this file.

 
#include <iostream>
#include "DGtal/base/Common.h"
#include "DGtal/helpers/StdDefs.h"
 
#include "DGtal/shapes/EuclideanShapesDecorator.h"
 
#include "DGtal/shapes/implicit/ImplicitBall.h"
 
 
using namespace DGtal;
 
// Functions for testing class EuclideanShapesDecorator.
 
bool testEuclideanShapesDecorator()
{
  unsigned int nbok = 0;
  unsigned int nb = 0;
 
  trace.beginBlock ( "Testing Unary operation on euclidean shapes ..." );
 
  typedef ImplicitBall< Z2i::Space > ShapeA;
  typedef ShapeA::RealPoint RealPoint;
 
  ShapeA shapeA(Z2i::RealPoint( -2.501, 0.0 ), 2.5013);
  ShapeA shapeB(Z2i::RealPoint( 2, 0.0 ), 2.5013);
  ShapeA shapeC(Z2i::RealPoint( 0.0, 0.0 ), 2.5);
 
  typedef EuclideanShapesCSG< ShapeA, ShapeA > CSG;
  CSG s_union ( shapeA );
  s_union.plus( shapeB );
 
  CSG s_intersec ( shapeA );
  s_intersec.plus( shapeB );
  s_intersec.intersection( shapeC );
 
  CSG s_minus ( shapeA );
  s_minus.minus( shapeC );
 
 
  nbok += (( s_union.orientation( RealPoint( -5.1, 0.0 )) == INSIDE )
          || ( s_union.orientation( RealPoint( -5.1, 0.0 )) == ON  ))? 0 : 1;
  nbok += (( s_union.orientation( RealPoint( 0.0, 0.0 )) == INSIDE )
          || ( s_union.orientation( RealPoint( 0.0, 0.0 )) == ON  ))? 1 : 0;
  nbok += (( s_union.orientation( RealPoint( 4.6, 0.0 )) == INSIDE )
          || ( s_union.orientation( RealPoint( 4.6, 0.0 )) == ON  ))? 0 : 1;
  nbok += (( s_union.orientation( RealPoint( 4.5, 0.0 )) == INSIDE )
          || ( s_union.orientation( RealPoint( 4.5, 0.0 )) == ON  ))? 1 : 0;
  nbok += (( s_union.orientation( RealPoint( -5.0, 0.0 )) == INSIDE )
          || ( s_union.orientation( RealPoint( -5.0, 0.0 )) == ON  ))? 1 : 0;
 
 
  nbok += (( s_intersec.orientation( RealPoint( -2.6, 0.0 )) == INSIDE )
          || ( s_intersec.orientation( RealPoint( -2.6, 0.0 )) == ON  ))? 0 : 1;
  nbok += (( s_intersec.orientation( RealPoint( 2.5, 0.0 )) == INSIDE )
          || ( s_intersec.orientation( RealPoint( 2.5, 0.0 )) == ON  ))? 1 : 0;
  nbok += (( s_intersec.orientation( RealPoint( 2.6, 0.0 )) == INSIDE )
          || ( s_intersec.orientation( RealPoint( 2.6, 0.0 )) == ON  ))? 0 : 1;
  nbok += (( s_intersec.orientation( RealPoint( 0.0, 2.5 )) == INSIDE )
          || ( s_intersec.orientation( RealPoint( 0.0, 2.5 )) == ON  ))? 0 : 1;
 
 
  nbok += (( s_minus.orientation( RealPoint( -5.0, 0.0 )) == INSIDE )
          || ( s_minus.orientation( RealPoint( -5.0, 0.0 )) == ON  ))? 1 : 0;
  nbok += (( s_minus.orientation( RealPoint( -1.0, 0.0 )) == INSIDE )
          || ( s_minus.orientation( RealPoint( -1.0, 0.0 )) == ON  ))? 0 : 1;
  nbok += (( s_minus.orientation( RealPoint( 0.0, 2.0 )) == INSIDE )
          || ( s_minus.orientation( RealPoint( 0.0, 2.0 )) == ON  ))? 0 : 1;
  nbok += (( s_minus.orientation( RealPoint( -2.6, 0.0 )) == INSIDE )
          || ( s_minus.orientation( RealPoint( -2.6, 0.0 )) == ON  ))? 1 : 0;
 
 
  nb = 13;
 
 
  trace.info() << "(" << nbok << "/" << nb << ") "
               << "true == true" << std::endl;
  trace.endBlock();
  return  nbok == nb;
}
 
// Standard services - public :
 
int main( int argc, char** argv )
{
  trace.beginBlock ( "Testing class EuclideanShapesDecorator" );
  trace.info() << "Args:";
  for ( int i = 0; i < argc; ++i )
    trace.info() << " " << argv[ i ];
  trace.info() << std::endl;
 
  bool res = testEuclideanShapesDecorator(); // && ... other tests
  trace.emphase() << ( res ? "Passed." : "Error." ) << std::endl;
  trace.endBlock();
  return res ? 0 : 1;
}
//                                                                           //