Astroid2D.h

 
#pragma once

 
#if defined(Astroid2D_RECURSES)
#error Recursive header files inclusion detected in Astroid2D.h
#else // defined(Astroid2D_RECURSES)
#define Astroid2D_RECURSES
 
#if !defined Astroid2D_h
#define Astroid2D_h

// Inclusions
#include <iostream>
#include "DGtal/base/Common.h"
#include "DGtal/shapes/parametric/StarShaped2D.h"
#include <cmath>
 
namespace DGtal
{

  // template class Astroid2D
  template <typename TSpace>
  class Astroid2D : public StarShaped2D<TSpace>
  {
    // ----------------------- Standard services ------------------------------
    public:
 
    typedef TSpace Space;
    typedef typename Space::RealPoint RealPoint;
    typedef typename Space::RealVector RealVector;

    Astroid2D() = delete;

    Astroid2D( const double x0, const double y0,
              const double a, const double b );

    Astroid2D( const RealPoint& aPoint, const double a, const double b );

    Astroid2D( const Astroid2D& other );

    Astroid2D& operator= ( const Astroid2D& other ) = delete;

    ~Astroid2D() = default;
 
  // ------------- Implementation of 'StarShaped' services -------------------
  public:

    RealPoint getLowerBound() const
    {
      return RealPoint( myCenter[0] - myA, myCenter[1] - myB );
    }

    RealPoint getUpperBound() const
    {
      return RealPoint( myCenter[0] + myA, myCenter[1] + myB );
    }

    RealPoint center() const
    {
      return myCenter;
    }

    inline
    void moveTo( const RealPoint& newCenter )
    {
      myCenter = newCenter;
    }

    double parameter( const RealPoint & p ) const;
 

    RealPoint x( const double t ) const;

    RealVector xp( const double t ) const;

    RealVector xpp( const double t ) const;
 
 
  // ------------------------- data -----------------------------------------
  private:

    RealPoint myCenter;

    double myA;

    double myB;
 
  // ----------------------- 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 Astroid2D
 
 

   template <typename T>
   std::ostream&
   operator<< ( std::ostream & out, const Astroid2D<T> & object );
 
} // namespace DGtal
 

// Includes inline functions.
#include "DGtal/shapes/parametric/Astroid2D.ih"
 
//                                                                           //
 
#endif // !defined Astroid2D_h
 
#undef Astroid2D_RECURSES
#endif // else defined(Astroid2D_RECURSES)