StarShaped2D.h

 
#pragma once
 
#if defined(StarShaped2D_RECURSES)
#error Recursive header files inclusion detected in StarShaped2D.h
#else // defined(StarShaped2D_RECURSES)
 
#define StarShaped2D_RECURSES
 
#if !defined StarShaped2D_h
 
#define StarShaped2D_h
 
// Inclusions
#include <iostream>
#include <vector>
#include "DGtal/base/Common.h"
#include "DGtal/kernel/NumberTraits.h"
 
namespace DGtal
{
 
  // template class StarShaped2D
  template <typename TSpace>
  class StarShaped2D
  {
 
  public:
    typedef TSpace Space;
    typedef typename Space::RealPoint RealPoint;
 
    StarShaped2D() = default;
 
    StarShaped2D ( const StarShaped2D & other ) = delete;
 
    StarShaped2D & operator= ( const StarShaped2D & other ) = delete;
 
    virtual ~StarShaped2D() = default;
 
  // ------------------------- Implemented services -------------------------
  public:
    virtual RealPoint interiorPoint() const
    {
      return center();
    }
 
  // ------------------------- Abstract services ----------------------------
  public:
 
    virtual RealPoint getLowerBound() const = 0;
 
    virtual RealPoint getUpperBound() const = 0;
 
 
    virtual RealPoint center() const = 0;
 
    virtual void moveTo( const RealPoint& newCenter ) = 0;
 
    virtual double parameter( const RealPoint & p ) const = 0;
 
    virtual RealPoint x( const double t ) const = 0;
 
    virtual RealPoint xp( const double t ) const = 0;
 
    virtual RealPoint xpp( const double t ) const = 0;
 
 
  // ------------------------- star-shaped services -------------------------
  public:
 
    Orientation orientation( const RealPoint &p) const;
 
    RealPoint tangent( const double t ) const;
 
    RealPoint normal( const double t ) const;
 
    double curvature( const double t ) const;
 
    double arclength( const double t1, double t2, const unsigned int nb ) const;
 
    RealPoint findIntersection( const RealPoint& inner, const RealPoint& outer, const double epsilon ) const;
 
    RealPoint closestPointWithWitnesses( const RealPoint& p, const RealPoint& left, const RealPoint& right, const int step) const;
 
  // ----------------------- Interface --------------------------------------
  public:
 
    void selfDisplay ( std::ostream & out ) const;
 
    bool isValid() const;
 
  // ------------------------- Hidden services ------------------------------
  private:
 
    template <typename T>
    inline
    bool isAlmostEqual( T x, T y ) const
    {
      return std::abs(x - y) <= std::numeric_limits<T>::epsilon();
    }
 
  }; // end of class StarShaped2D
 
 
  template <typename T>
  std::ostream&
  operator<< ( std::ostream & out, const StarShaped2D<T> & object );
 
} // namespace DGtal
 
 
// Includes inline functions.
#include "DGtal/shapes/parametric/StarShaped2D.ih"
 
//                                                                           //
 
#endif // !defined StarShaped2D_h
 
#undef StarShaped2D_RECURSES
#endif // else defined(StarShaped2D_RECURSES)