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)