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)