testShapesFromPoints.cpp

#include <iostream>
#include <fstream>
#include <sstream>
#include <exception>
#include <boost/program_options.hpp>

#include "DGtal/base/Common.h"
#include "DGtal/kernel/PointVector.h"
#include "DGtal/shapes/fromPoints/CircleFrom3Points.h"
#include "DGtal/shapes/fromPoints/CircleFrom2Points.h"

#include "DGtal/io/boards/Board2D.h"

#include "ConfigTest.h"


using namespace std;
using namespace DGtal;

// testing Functions

template <typename Coordinate>
bool testCircleFrom3Points()
{

  trace.beginBlock("Simple test for CircleFrom3Points"); 
  
  typedef PointVector<2,Coordinate> Point; 
  CircleFrom3Points<Point> c; 
  Point o(0,0); 

  {
    c.init( Point(0,1), Point(150,18), Point(250,-48) ); 
    trace.info() << c << endl;
    trace.info() << o << " is at distance " << c.signedDistance(o) << endl;
    if (c.signedDistance(o) != -4026300) return false; 
  }
  {
    c.init( Point(0,1), Point(150,18), Point(100,48) ); 
    trace.info() << c << endl;
    trace.info() << o << " is at distance " << c.signedDistance(o) << endl;
    if (c.signedDistance(o) != 442200) return false; 
  }
  {
    c.init( Point(8,4), Point(9,3), Point(10,0) ); 
    trace.info() << c << endl;
    trace.info() << o << " is at distance " << c.signedDistance(o) << endl;
    if (c.signedDistance(o) != 0) return false;
    double cx, cy, r; 
    c.getParameters(cx, cy, r); 
    if ( (cx != 5.0)||(cy != 0.0) )
    {
      trace.emphase() << "error in returned center" << endl;
      trace.info() << "(" << cx << "," << cy << ") " << endl;
      trace.info() << "instead of (5,0)" << endl;
      return false; 
    } 
  }
  {
    c.init( Point(0,0), Point(150,20), Point(15,2) ); 
    trace.info() << c << endl;
    trace.info() << o << " is at distance " << c.signedDistance(o) << endl;
    if (c.signedDistance(o) != 0) return false; 
  }

  Board2D board; 
  board << SetMode(o.className(), "Grid") << o; 
  {
    c.init( Point(8,4), Point(9,3), Point(10,0) ); 
    trace.info() << "arc drawing" << endl;
    
    board << SetMode( c.className(), "Arc" );
    Display2DFactory::draw(board, c, Point(5,10), Point(8,4));
    
    board << SetMode( c.className(), "Sector" );
    Display2DFactory::draw(board, c, Point(9,3), Point(10,0) );
    
    board << SetMode( c.className(), "Annulus" );
    Display2DFactory::draw(board, c, Point(5,-10), Point(2,-4) );
    
    board.saveEPS("arcDisplay.eps");
  }


  trace.endBlock(); 
  
  return true; 
}

template <typename Coordinate>
bool testCircleFrom2Points()
{

  trace.beginBlock("Simple test for CircleFrom2Points"); 
  
  typedef PointVector<2,Coordinate> Point; 
  Point o(0,0);
  Point pole(0,1); 
  CircleFrom2Points<Point> c( pole ); 

  {
    c.init( Point(15,2), Point(10,5) ); 
    trace.info() << c << endl;
    trace.info() << o << " is at distance " << c.signedDistance(o) << endl;
    if (c.signedDistance(o) != 470) return false; 

  }

  trace.endBlock(); 
  
  return true; 
}


// Standard services - public :

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


  bool res = testCircleFrom3Points<int>()
  && testCircleFrom3Points<double>()
#ifdef WITH_BIGINTEGER
  && testCircleFrom3Points<BigInteger>()
#endif
  && testCircleFrom2Points<int>()
; 


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