Ellipse2D.h

 
#pragma once
 
#if defined(Ellipse2D_RECURSES)
#error Recursive header files inclusion detected in Ellipse2D.h
#else // defined(Ellipse2D_RECURSES)
#define Ellipse2D_RECURSES
 
#if !defined Ellipse2D_h
#define Ellipse2D_h
 
// Inclusions
#include <iostream>
#include "DGtal/base/Common.h"
#include "DGtal/shapes/parametric/StarShaped2D.h"
#include <cmath>
 
namespace DGtal
{
 
  // template class Ellipse2D
  template <typename TSpace>
  class Ellipse2D final:  public StarShaped2D<TSpace>
  {
  // ----------------------- Standard services ------------------------------
  public:
 
    typedef TSpace Space;
    typedef typename Space::RealPoint RealPoint;
    typedef typename Space::RealVector RealVector;
 
    Ellipse2D() = delete;
 
    Ellipse2D( const double x0, const double y0,
         const double a0, const double a1, const double theta );
 
    Ellipse2D( const RealPoint& aPoint,
        const double a0, const double a1, const double theta );
 
     Ellipse2D( const Ellipse2D& other );
 
    Ellipse2D& operator=( const Ellipse2D& other ) = delete;
 
    ~Ellipse2D() = default;
 
  // ------------- Implementation of 'StarShaped' services ------------------
  public:
 
    RealPoint getLowerBound() const
    {
      return myCenter - myAxis1;
    }
 
    RealPoint getUpperBound() const
    {
      return myCenter + myAxis1;
    }
 
    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 myAxis1;
 
 
    double myAxis2;
 
    double myTheta;
 
  // ----------------------- 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 Ellipse2D
 
 
  template <typename T>
  std::ostream&
  operator<< ( std::ostream & out, const Ellipse2D<T> & object );
 
} // namespace DGtal
 
 
// Includes inline functions.
#include "DGtal/shapes/parametric/Ellipse2D.ih"
 
//                                                                           //
 
#endif // !defined Ellipse2D_h
 
#undef Ellipse2D_RECURSES
#endif // else defined(Ellipse2D_RECURSES)