testHalfPlane.cpp

Go to the documentation of this file.

#include <cstdio>
#include <cmath>
#include <fstream>
#include <vector>
#include <iostream>
#include <iterator>
#include "DGtal/base/Common.h"
#include "DGtal/base/Exceptions.h"
#include "DGtal/kernel/SpaceND.h"
#include "DGtal/kernel/domains/DomainPredicate.h"
#include "DGtal/kernel/domains/HyperRectDomain.h"
#include "DGtal/io/boards/Board2D.h"
#include "DGtal/shapes/fromPoints/StraightLineFrom2Points.h"
#include "DGtal/shapes/fromPoints/Point2ShapePredicate.h"

#include "DGtalCatch.h"


using namespace DGtal;
using namespace DGtal::functors;
using namespace LibBoard;

TEST_CASE( "HalfplaneUnit tests" )
{
  typedef int Coordinate;
  typedef PointVector<2,Coordinate> Point;
  typedef StraightLineFrom2Points<Point> StraightLine;

  Point p(0,0);
  Point q(5,2);
  Point r1(2,0);
  Point r2(0,2);
  Point r3(10,4);
  StraightLine line(p,q);
  
  //halfplane upward oriented and closed 
  typedef Point2ShapePredicate<StraightLine,true,true> UpClosedHalfPlane;
 
  //halfplane upward oriented and open 
  typedef Point2ShapePredicate<StraightLine,true,false> UpOpenHalfPlane; 
  //halfplane downward oriented and closed 
  typedef Point2ShapePredicate<StraightLine,false,true> DownClosedHalfPlane; 

  //Location
  trace.beginBlock("Creation of a halfplane and test some points");
  std::string res;
  trace.info() << "created line:" << std::endl;
  trace.info() << line << std::endl;

  //line is an instance of StraightLine
  UpClosedHalfPlane hp(line);
  trace.info() << hp << std::endl;
    trace.info() << "Does " << r1 << " belongs to hp (no)?" << std::endl;
  REQUIRE( !hp(r1) );
  REQUIRE(hp(r2));
  REQUIRE(hp(r3));


  UpOpenHalfPlane hp2(line);
  REQUIRE(! hp2(r1));
  REQUIRE(hp2(r2));
  REQUIRE( ! hp2(r3) );

  DownClosedHalfPlane hp3(line);
  REQUIRE(hp3(r1));
  REQUIRE(!hp3(r2));
  REQUIRE(hp3(r3));
}